AD-A064  728  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OHIO  SCH— ETC  F/G  9/2 


microprocessor  based  data  acquisition  and  PROCESSING  SYSTEM. (U) 
DEC  78  S IFTEKHAR 

UNCLASSIFIED  AFIT/GE/EE/78-29  NL 


> >-* 
Q- 
O 
CJ> 


a 

t l-L- 


MICROPROCESSOR  BASED  DATA  ACQUISITION 
AND  PROCESSING  SYSTEM 

THESIS 


AFIT/GE/EE/78-29 


L 

1 1 


Saleem  Iftekhar 
Flt/Lt  PAF 


K 


Approved  for  public  release;  distribution  unlimited. 


A i 


I i' 


r i2  . 

^AFIT/GE/EE/78-29 


Ck 


MICROPROCESSOR  BASED  DATA  ACQUISITION 
AND  PROCESSING  SYSTEM, 


A)^4-efV  THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


C Jj?  Sa  1 eem / 1 ftekhar 

F+t/tt PAF 

Graduate  Electrical  Engineering 

SiT',  Decartwr  $*78 


92A 


Approved  for  public  release;  distribution  unlimited. 


r 


Preface 


This  report  is  the  result  of  my  attempt  to  develop  a microprocessor 
based  data  acquisition  and  processing  system  for  the  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio.  The  report 
describes  the  three  phases  of  system  development:  conceptual  phase, 
design  phase,  and  implementation  and  test  phase.  Structured  analysis  and 
design  technique  was  used  in  the  conceptual  and  design  phases.  In  my 
opinion,  this  graphical  technique  makes  it  possible  to  follow  a top-down 
design  approach.  A bench-type  model  was  the  result  of  the  implementation 
and  test  phase. 

I wish  to  express  my  appreciation  to  Dr.  Gary  B.  Lamont,  my  thesis 
advisor,  for  his  many  hours  of  help  and  advise.  I greatly  appreciate  his 
detailed  reading  of  the  report  to  bring  to  light  errors  which  I had  made 
and  his  suggestions  for  improvement.  I would  like  to  thank  the  committee 
members,  Capt  Gregg  L.  Vaughn  and  Dr.  Thomas  C.  Hartrum  for  their  valuable 
guidance  during  the  course  of  the  project. 

Also,  I would  like  to  thank  Dr.  Lynn  C.  Rogers  and  Mr.  Robert  W. 
Gordon  of  the  Air  Force  Flight  Dynamics  Laboratory  who  provided  the 
thesis  topic.  Their  interest  in  my  work  and  assistance  during  the  pro- 
ject were  invaluable. 

I would  like  to  mention  my  appreciation  to  Lieutenant  Peter  D. 
Summers,  III,  and  Mr.  Frank  E.  Beitel  of  the  Air  Force  Materials  Labora- 
tory, Computer  Activities  Division,  Wright-Patterson  Air  Force  Base,  Ohio, 
for  extending  their  help  during  the  software  development. 

The  manuscript  in  its  present  form  would  not  be  possible  without 
the  efforts  of  my  typist.  Miss  Joyce  L.  Wilson. 
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Abstract 

The  Air  Force  Flight  Dynamics  Laboratory  at  Wright-Patterson  Air 
Force  Base,  Ohio,  plans  to  design  a constrained  layer  damping  treatment, 
for  aircraft  components.  This  design  requires  the  knowledge  of  vibra- 
tion and  temperature  variations  that  a component  encounters  in  service.. 

The  development  of  a microprocessor  based  data  acquisition  and  pro- 
cessing system  (MIBDAPS)  to  monitor  the  vibration  and  temperature  varia- 
tions of  a component  under  test  is  discussed  in  this  paper.  The  designed 
system  is  capable  of  acquiring  and  processing  data  in  real  time  and  gen- 
erating cumulative  damage  factor  versus  frequency  and  temperature  con- 
tours. These  contours  are  used  in  the  design  of  the  damping  treatment* 

The  system  uses  an  LSI-1 1 M microcomputer  as  its  central  processor. 

The  support  hardware  for  this  microcomputer  includes  an  analog-to-digital 
converter,  core  memory,  real  time  clock,  and  serial  line  interface  module. 
The  software  developed  for  this  system  handles  the  acquisition  and  pro- 
cessing of  data.  The  real  time  processing  includes  the  computation  of  the 
FFT  and  power  spectrum  of  the  vibration  signal  and  subsequent  calculation 
of  the  cumulative  damage  factor  versus  frequency  and  temperature  contours. 
The  bench  model  of  the  designed  system  was  implemented  and  tested  for 
functional  performance. 


I.  Introduction 


The  Air  Force  Flight  Dynamics  Laboratory  (AFFD),  Wright-Patterson 
Air  Force  Base,  Ohio,  plans  to  develop  a constrained  layer  damping  treat- 
ment for  aircraft  components  such  as  the  airframe,  avionics  hardware,  and 
other  accessories.  The  specific  damping  layer  treatment  design  requires 
knowledge  of  temperature  and  vibration  variations  encountered  by  the  com- 
ponent. To  generate  this  information,  a data  acquisition  and  processing 
system  is  required  to  serve  as  an  in-flight  monitor  for  the  component 
under  test.  This  report  describes  the  design,  development,  and  bench 
model  implementation  of  such  a system. 


Background 

Research  has  established  that  high-cycle  fatigue  damage  due  to 
resonant  vibration  is  one  of  the  major  reasons  for  structural  failure  of 
components  (Ref  13).  A constrained  layer  damping  treatment  would  be  cost 
effective  in  extending  the  life  of  the  component.  In  order  to  determine 
the  correct  type  of  damping  layer,  the  temperatures,  vibration  nodes, 
and  frequencies  at  which  damage  occurs  in  service  must  be  known. 

It  has  been  determined  that  there  are  approximately  ten  vibration 
frequency  bands  of  interest  (Ref  13).  They  lie  in  the  range  from  100  HZ 
to  1000  HZ  and  are  spaced  100  HZ  apart  (Figure  1).  A damage  factor  (DF) 
at  each  frequency  of  interest  can  be  calculated  using  the  following 


relationship  (Ref  13): 


DF  = 


K(A) 


(f) 


23-1 


(1) 


1 


Here,  3 is  approximately  equal  to  3-2,  and  A is  the  root-mean-square 
(RMS)  value  of  acceleration  in  the  frequency  band  of  interest. 


A summation  of  the  DF  over  a predetermined  period  of  time  would  yield  a 
cumulative  damage  factor  (CDF)  and  curves  as  shown  in  Figure  2 could  be 
drawn. 

Existing  in-flight  data  acquisition  techniques  employ  an  analog 
signal  conditioner  and  a multichannel  analog  tape  recorder  (Ref  13).  The 
analog  signals  from  the  vibration  and  temperature  sensors  are  conditioned 
and  recorded  on  a magnetic  tape.  This  system  is  bulky  and  takes  excessive 
setup  time.  Skilled  personnel  are  required  to  install  the  equipment  on 
an  aircraft,  and  calibration  is  difficult  and  time-consuming. 

To  simplify  this  process,  the  design  of  a self-contained  Micropro- 
cessor Based  Data  and  Processing  System  (MIBDAPS)  appeared  as  an  attrac- 
tive alternative.  Such  a system  should  be  easy  to  install  and  remove, 
and  only  minimum  operator  intervention  during  its  operation  should  be 
necessary. 

System  Concept 

The  underlying  concept  of  MIBDAPS  was  a self-contained  system  that 
would  monitor  the  analog  signals  from  the  temperature  and  vibration  sen- 
sors, digitize  these  signals,  process  the  data  in  real  time,  and  store 
the  results  on  an  appropriate  media.  The  real  time  processing  would 
include  computation  of  the  power  spectrum  of  the  vibration  signal  and 
calculation  of  DF  at  each  frequency  of  interest.  The  results  of  process- 
ing would  be  stored  in  a CDF  versus  frequency  and  temperature  array,  respec- 
tively. Therefore,  at  the  termination  of  system  operation,  the  CDF  versus 
frequency  and  temperature  contours  (Figure  2)  should  be  available. 

MIBDAPS  should  be  capable  of  either  15  to  20  hours  of  continuous 
operation  (for  gathering  data  on  transport  aircraft  hardware)  or  time 
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segmented  operation  (for  gathering  data  on  fighter  aircraft  hardware  which 
has  limited  flying  endurance).  During  a time  segmented  operation,  the  sys- 
tem should  be  able  to  operate  for  a total  aggregate  of  15  to  20  hours.  The 
data  processed  during  this  time  interval  would  enable  determination  of 
vibration  nodes  at  which  the  damage  occurs. 

Approach 

The  development  of  MIBDAPS  was  carried  out  in  three  phases:  the  con- 
ceptual phase,  the  design  phase,  and  the  implementation  and  testing  phase. 
During  the  conceptual  phase,  the  detailed  requirements  placed  on  the  system 
were  defined  using  a structured  analysis  and  design  technique  (SADT)  (Ref  1). 
The  design  phase  included  the  software/hardware  partition  of  the  system, 
selection  of  hardware  components,  and  design  of  software  and  hardware 
modules  to  implement  different  functions  of  MIBDAPS.  During  the  implemen- 
tation and  testing  phase,  the  software  and  hardware  modules  were  given  a 
functional  check.  The  system  integration  and  performance  check  of  the 
integrated  system  was  also  completed  during  this  phase. 

The  Development  Environment 

The  LSI-11M*  microcomputer  was  used  as  the  system  processor  because»of 
its  powerful  instruction  set,  speed,  and  availability  of  support  hardware' 
and  software.  The  use  of  the  LSI-11  microcomputer  system  for  MIBDAPS 
determined  the  general  development  requirements  for  both  hardware  and  soft- 
ware modules  of  MIBDAPS.  The  off-the-shelf  hardware  modules  used  where 
analog-to-digital  converter  (ADC),  real  time  clock  (RTC),  and  the  4K  by 
16-bits  core  memory.  All  of  these  modules  are  manufactured  by  the  Digital 
Equipment  Corporation  (DEC). 


♦LSI-1 1M  is  the  militarized  version  that  is  manufactured  by  Norden  Systems,  Inc. 
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Figure  2.  CDF  Versus  Frequency  at  Different  Temperature  Values 
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The  software  development  system  used  was  a PDP-11/03  minicomputer 
which  has  the  LSI-11  as  the  processor.  The  RXV011  floppy  disk  based 
operating  system  was  available  on  this  minicomputer.  The  development 
system  was  located  at  the  Air  Force  Materials  Laboratory  (AFML),  Wright- 
Patterson  Air  Force  Base,  Ohio. 

Scope  of  Thesis 

The  structured  analysis  and  design  technique  was  used  to  develop  the 
activity  model  of  MIBDAPS  (Ref  Chap  II),  and  the  LSI-1 1 M microcomputer 
was  used  as  the  processor  for  MIBDAPS  (Ref  Chap  III).  Hardware  and  soft- 
ware modules  were  designed  and  tested  (Ref  Chaps  IV  and  V).  The  hardware/ 
software  modules  were  then  integrated  into  a bench  model  of  MIBDAPS  which 
was  then  functionally  tested  (Ref  Chap  VI).  The  bench  model  of  MIBDAPS 
performed  the  required  functions  of  data  acquisition  and  processing. 
Testing  of  the  integrated  system  for  computational  accuracy  has  not  been 
completed  at  this  time.  Recommendations  for  further  work  have  been  gener- 
ated in  order  to  complete  the  implementation  of  the  prototype  module  (Ref 
Chap  VII). 
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II.  System  Requirements  for  MIBDAPS 


Introduction 

The  conceptual  idea  is  not  enough  to  form  the  basis  of  design  and 
implementation  of  any  complex  system.  The  design  can  be  facilitated  if 
the  complete  system  is  divided  into  functional  subunits.  These  sub- 
units are  easier  to  analyze  and  design.  The  division  of  a system 
into  functional  subunits  according  to  the  requirements  placed  on  the 
system  can  be  called  the  "requirement  definition"  phase  of  the  design. 

A carefully  prepared  requirement  definition  ensures  that  ambiguities  in 
the  design  and  implementation  phase  are  minimized. 

Analysis  and  Design  Methodology 

Structured  Analysis  and  Design  Techniques  (SADT)  is  a design 
methodology  for  performing  the  functional  analysis  and  design  of  complex 
systems  (Ref  1).  A functional  model  of  the  system  is  created  first 
which  helps  in  understanding  what  activities  the  system  will  perform, 
followed  by  a design  model  to  show  how  the  system  will  be  implemented 
to  perform  the  requisite  functions  (Ref  1).  The  functional  model  also 
helps  to  clarify  the  requirement  placed  on  the  systems 

Appendix  A describes  SADT  methodology  and  its  application  rules. 
Some  changes  in  convention  are  made  to  suit  the  design  of  MIBDAPS;  these 
are  also  included  in  Appendix  A. 

The  SADT  was  chosen  as  a functional  analysis  and  design  tool 
because  it  is  a graphical  method  of  analyzing  the  functional  require- 
ments of  a system.  The  activity  model,  which  shows  different  function* 
which  the  system  is  required  to  perform,  is  independent  of  the  hardware/ 
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software  implementation  of  the  system.  Thus,  any  changes  in  the 
activity  model  can  be  made  during  this  phase  of  the  design  without 
affecting  the  final  shape  of  the  product. 

Once  the  activity  model  is  complete,  it  is  easy  to  review  and 
correct  any  oversights  that  might  have  occurred  during  the  functional 
analysis  of  the  system. 

Activity  Model 

Before  the  activity  model  is  drawn,  the  system  is  divided  into 
different  nodes.  The  node  index  given  in  Table  I shows  the  hierarchy  of 
the  system  from  activity  viewpoint.  This  node  index  is  based  on  the 
initial  problem  definition  by  the  Air  Force  Avionics  Laboratory  (AFAL), 

Wright-Patterson  Air  Force  Base,  Ohio,  and  subsequent  discussions  with 
the  sponsor  of  this  thesis.  The  node  index  also  indicates  how  the 
author  visualized  the  system  was  required  to  function-. 

TABLE  I 


Node  Index 


A-l 

Data  Collection  and  Analysis 

A-0 

Acquire  and  Process  Data 

Acquire  and  Process  Data 

A1 

Control  All  Activity 

A2 

Acquire  Data 

A3 

Process  Data 

M 

Store  Data 

The  explanation  of  each  node  follows. 
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Node  A-1,  Data  Collection  and  Analysis  (Figure  3).  This  node  shows 
the  overall  research  activity  and  the  place  where  MI6DAPS  fits  in  the 
overall  function  of  the  organization.  Block  3,  Plot  Data,  activity  will 
be  carried  out  in  the  laboratories. 

Node  A-0,  Acquire  and  Process  Data  (Figure  4).  Node  A-0  is  the 
model  of  the  MIBDAPS.  The  node  shows  the  system  requirements.  Data 
input  comes  in  the  form  of  analog  signals— analog  vibration  signal  (II) 
and  analog  temperature  signal  (12).  The  output  03  of  the  system  is  a 
three-dimensional  array.  This  array  contains  the  cumulative  damage 
factor  (CDF)  as  a function  of  frequency  and  temperature.  Output  02  is 
the  scale  factor  for  the  array  elements,  and  output  01  is  the  total  opera- 
tional time  for  which  the  system  was  operational. 

The  MIBDAPS  will  be  able  to  continuously  acquire  and  process  data  for 
a total  flying  time  of  15  to  20  hours.  This  is  the  requirement  placed  by 
AFAL  because  the  CDF  versus  temperature  and  frequency  array  for  this  dura- 
tion will  be  a good  estimate  of  the  fatigue  damage  profile  of  the  test 
piece  (Ref  13).  Input  signal  specifications  are: 

Analog  Vibration  Signal  (II)  0-5V 

Analog  Temperature  Signal  (12)  0-5V 

Node  A0,  Acquire  and  Process  Data  (Figure  5).  The  MIBDAPS  is 
expected  to  acquire  data,  process  data,  and  store  the  results.  There- 
fore, node  A0  shows  the  decomposition  of  the  system  into  four  primary 
activities: 

Acquire  Data 
Process  Data 
Store  Data 

Control  All  Processes 
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An  explanation  of  how  a signal  travels  through  each  subnode  is  now 
given.  The  activity  Acquire  Data  continuously  digitizes  the  input  signal 
under  the  control  from  node  A1 . The  digitized  signal  is  queued  and  sent 
to  the  processing  node  (A2). 

The  power  spectrum  of  the  signal  is  computed  by  the  processing  node. 
For  each  frequency  of  interest,  a damage  factor  is  calculated  (Ref  Chap  I) 
and  stored  in  an  array.  The  activity  Store  Data  adds  the  most  recent  set 
of  DF  values  to  the  previously  stored  values.  This  results  in  a storage 
of  cumulative  damage  factors. 

Node  A01 , Control  All  Processes,  keeps  all  events  in  the  system 
sychronized.  It  initiates  appropriate  commands  to  ensure  that  data  is 
processed  and  stored  in  proper  order.  The  input  and  output  signal  speci- 
fications are  given  below: 

- Input  Signal 

a.  Analog  Vibration  Signal  (II):  0-5V. 

b.  Analog  Temperature  Signal  (12):  0-5V. 

c.  AC  Supply  (13):  110V,  60  HZ. 

- Output  Signal 

a.  Operational  Time  (01):  Displayed  on  the  front  panel. 

b.  Scale  Factor  (02). 

c.  CDF  Versus  Temperature  and  Frequency  Array  (03). 

Node  A1 , Control  All  Activity  (Figure  6).  Node  A1  subdivides  the 
overall  control  of  the  system  into  the  following  functions: 

Control  Acquisition  of  Data 
Control  Processing  of  Data 
Control  Storage  of  Data 
Control  Power  to  the  System 
Record  Operational  Time 
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Figure  6.  Node  AT,  Control  All  Activity 


Control  signal  Cl  (switch  ON/OFF)  activates  the  system.  The  power 
to  the  system  is  also  controlled  by  total  operational  time;  therefore, 
at  the  end  of  15  to  20  hours  of  operation,  the  system  halts  automatically 
(Ref  13).  The  output  and  control  signals  are  listed  below  along  with 
their  explanation. 

- Output  Signal 

a.  Operational  Time  (01):  Four  binary  coded  decimal  (BCD)  num- 
bers, each  4-bit  long,  display  the  operational  time  elapsed  in  hours  and 
minutes. 

b.  Data  Acquisition  Command  (0202):  This  control  signal  initiates 
the  acquisition  of  N-data  points.  This  signal  acts  as  a control  signal 

to  node  A2. 

c.  Process  Data  Command  (0302):  This  signal  starts  the  process- 
ing of  either  the  raw  data  or  update  of  CDF  versus  temperature  and  fre- 
quency array. 

d.  Store  Data  Command  (0403):  This  command  signal  initiates 
either  the  storage  of  the  N-processed  data  element  or  storage  of  updated 
CDF  versus  temperature  and  frequency  array. 

- Control  Signal 

a.  Switch  ON/OFF  Control  (Cl):  This  is  a logic  "1"  when  the 
input  switch  is  closed  and  a logic  "0"  when  the  input  switch  is  open. 

b.  Acquisition  Done  Status  (02C2):  This  signal  indicates  node 
A14  whenever  data  acquisition  of  N-data  points  is  complete. 

c.  Process  Done  Status  (02C3);  This  control  signal  indicates  the 
completion  of  processing  by  node  A3.  The  processing  is  initiated  by  node 


d.  Storage  Done  Status  (04C4):  Whenever  the  storage  of  data  is 
required,  node  A16  initiates  this  activity.  Once  the  storage  of  data 


is  complete,  this  signal  indicates  the  status. 

Node  A2,  Acquire  Data  (Figure  7).  This  node  converts  the  analog 
vibration  and  analog  temperature  signals  to  an  equivalent  digital  signal. 
The  acquisition  of  data  is  initialized  by  node  A1 , and  either  the  analog 
vibration  signal  or  the  analog  temperature  signal  is  selected.  The 
selected  signal  is  digitized,  and  the  final  output  is  an  equivalent  digi- 
tal value  of  the  input  signal  along  with  a signal  identifier.  The  input/ 
output  and  control  signals  specifications  are: 

- Input  Signal 

a.  Analog  Vibration  Signal  (II):  0-5V. 

b.  Analog  Temperature  Signal  (12):  0-5V. 

- Output  Signal 

a.  Acquisition  Done  Status  (0101):  This  signal  indicates  the 
status  of  the  data  acquisition  process.  Whenever  N-data  elements  have 
been  acquired,  this  status  signal  informs  node  A1  about  the  completion 
of  data  acquisition. 

b.  Formatted  Data  (0302):  This  is  the  digitized  data  along 
with  an  identifier  which  indicates  whether  the  data  point  corresponds  to 
the  analog  vibration  signal  or  the  analog  temperature  signal. 

- Control  Signal 

Data  Acquisition  Control  (01  Cl ) : This  signal  initiates  the  acqui- 
sition of  N-data  elements.  This  signal  controls  node  A2  in  the  following 
manner: 

(1)  It  selects  either  the  analog  vibration  signal  or  the  analog 
temperature  signal. 
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Figure  7.  Node  A2,  Acquire  Data 


r 1 

(2)  It  initiates  digitizing  of  N-samples  of  an  analog  vibration 

signal . 

(3)  It  initiates  digitizing  of  one  sample  of  an  analog  tempera- 
ture signal  for  every  N-samples  of  an  analog  vibration  signal. 

(4)  It  initiates  formatting  of  each  digitized  data  point. 

Node  A3,  Process  Data  (Figure  8).  The  activities  for  this  node  are 
listed  below: 

a.  N-elements  of  raw  data  array  are  taken  from  the  raw  data  storage 
area  and  processed.  The  processing  involves  Fourier  transform  of  N-raw 
data  points.  The  power  spectrum  of  first  M-frequency  components  is  com- 
puted next  (M  is  less  than  N).  The  DF  is  calculated  for  each  of  the 
M-frequency  components. 

b.  M-elements  of  the  final  data  array  (CDF  versus  temperature  and 
frequency)  are  updated  whenever  a new  set  of  DF  array  is  available  after 
processing  of  N-raw  data  points. 

The  input/output  and  control  signals  specifications  are: 

- Input  Signal 

a.  Raw  Data  From  Storage  (0411):  The  N-element  array  from  the 
Raw  Data  Storage  area  is  used  by  the  Process  Raw  Data  module  to  give 
an  M-element  processed  array. 

b.  New  Data  Array  (0412):  This  is  the  data  array  which  has 
resulted  from  processing  N-element  raw  data  array. 

c.  Old  Data  Array  (0413):  CDF  versus  temperature  and  frequency 
array  which  has  to  be  updated  by  adding  in  new  data  array.  The  Update 
Data  Module  creates  updated  data  array. 
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- Output  Signal 

a.  Scale  Factor  (03):  This  is  the  value  which  tells  by  what 
amount  the  elements  in  the  final  arrays  have  been  scaled. 

b.  Processed  Data  (0404):  This  is  an  M-element  data  array  where 
the  Process  Raw  Data  module  stores  the  results  after  operating  on  raw 
data  array. 

c.  Updated  Data  Array  (0405):  This  is  an  M-element  array.  The 
value  of  CDF  as  a function  of  frequency  and  temperature  is  stored  in  this 
array.  Every  time  a new  set  of  M-processed  data  elements  is  available, 
the  final  array  is  updated  by  adding  these  new  values  to  the  old  array. 

- Control  Signal 

Process  Data  Command  (01  Cl):  This  command  initiates  the  process- 
ing of  N-raw  data  elements.  After  the  processing  is  completed,  the  CDF 
versus  temperature  and  frequency  array  is  updated.  The  update  takes 
place  by  adding  the  newly  processed  data  array  to  the  existing  final 
array  (CDF  versus  temperature  and  frequency). 

Node  A4,  Store  Data  (Figure  9).  This  node  also  performs  multiple 
activities  on  the  data.  The  different  activities  are  enumerated  below. 

a.  The  Formatted  Raw  Data  points  are  stored  sequentially  in  an 
N-point  array.  This  array  is  called  Raw  Data  Array,  and  it  serves  as  a 
temporary  storage  area. 

b.  Whenever  raw  data  has  to  be  processed,  N-elements  of  the  raw 
data  array  are  transferred  to  node  A3  for  processing. 

c.  The  processed  data  is  stored  in  M-elements  processed  data  array. 
This  storage  of  processed  data  points  is  temporary. 

d.  Whenever  the  final  output  arrays  have  to  be  updated,  M-data 
elements  from  the  Process  Data  Array  are  transferred  to  node  A3. 
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e.  The  Update  Final  Array  (0313)  is  sent  from  node  A3,  and  it 
is  stored  as  the  final  CDF  versus  temperature  and  frequency  array  (02). 

The  input/output  and  control  signal  specifications  are: 

- Input  Signal 

a.  Formatted  Data  Input  (0211):  This  is  the  formatted  digitized 
version  of  input  analog  signals.  The  data  points  at  this  stage  are  also 
known  as  Raw  Data  Elements. 

b.  M-Element  Processed  Array  (0312).  This  is  an  M-element  pro- 
cessed array,  which  is  the  output  from  node  A3.  This  array  is  generated 
every  time  N-raw  data  elements  are  processed. 

c.  Updated  Final  Array  (0313):  The  final  output  array  (CDF 
versus  frequency  and  temperature)  has  to  be  updated  whenever  a new  M- 
element  processed  array  is  ready.  This  input  is  the  result  of  updating 
the  existing  final  output  array. 

- Output  Signal 

a.  N-Raw  Data  Elements  (0303):  Node  A3  requires  N-raw  data 
elements  whenever  processing  has  to  be  done.  This  output  is  actually  an 
N-raw  data  element  array,  which  is  sent  to  node  A3  for  processing. 

b.  M-Processed  Data  Elements  (0304):  This  output  is  used  to 
update  the  final  output  array,  so  the  CDF  versus  frequency  and  tempera- 
ture is  continuously  updated. 

c.  Storage  Done  Status  (0105):  This  signal  is  used  to  synchro- 
nize different  storage  operations.  The  Control  All  Process,  node  A1 , 
is  informed  about  the  storage  status  of  node  A4. 

d.  Final  Output  Array  (02):  This  array  stores  CDF  as  a function 
of  frequency  and  temperature  and  is  the  final  output  for  MIBDAPS. 
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Summary 

In  this  chapter,  the  overall  requirements  placed  on  MIBDAPS  are 
defined.  All  possible  activities  which  MIBDAPS  is  expected  to  perform, 
along  with  their  interrelationships,  have  been  consistently  and  completely 
documented.  At  this  stage,  engineering  decisions  about  the  mechanization 
are  not  taken;  however,  *the  next  phase  of  the  design  would  be  the  engi- 
neering decision  regarding  the  mechanization  of  the  system.  Hardware  and 
software  subdivision  would  also  take  place  in  the  next  phase. 


III.  Software/Hardware  Division  of  MIBDAPS 


Introduction 

The  functional  analysis  model  for  MIBDAPS,  which  was  developed  in 
Chapter  II,  does  not  discuss  mechanization  of  the  system.  The  choice 
of  the  hardware  for  mechanization  of  MIBDAPS  will,  to  a large  extent,, 
determine  the  subdivision  of  the  system  functions.  In  this  chapter,  the 
engineering  decisions  that  define  the  realization  of  MIBDAPS  are  pre- 
sented. The  implementation  of  different  modules  of  MIBDAPS  and  the  sub- 
division of  different  modules  into  software  and  hardware  activities  are  also 
presented.  SADT  is  again  employed  in  this  development  of  MIBDAPS* 
activities. 

Choice  of  Microcomputer  System 

One  of  the  objectives  in  the  development  of  MIBDAPS  was  a self- 
contained  microprocessor  based  data  acquisition  and  processing  system. 

Thus,  there  were  two  possibilities  for  the  mechanization  of  MIBDAPS: 

i)  Selection  of  an  off-the-shelf  microprocessor  with  the  requisite 
word  length  and  instruction  set,  followed  by  the  development  of  a useable 
microcomputer  system  with  the  required  amount  of  read/write  memory  (RAM y, 
input/output  (I/O)  interface  circuits,  and  development  of  application 
software. 

ii)  Selection  of  an  industry  standard  microcomputer  board  with 
already  developed  I/O  interface  circuits,  on-board  memory  and  memory 
expansion  capability,  and  software  development  tools. 

At  this  stage,  the  approach  taken  for  implementation  of  MIBDAPS  was, 
in  part,  guided  by  the  time  available  for  the  project.  The  first 
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possibility  where  the  system  could  be  developed  around  a selected  micro- 
processor would  encompass  design  of  I/O  interfaces,  memory  boards,  and 
finally  the  development  of  software.  This  would  have  taken  longer  than 
the  allotted  time.  Therefore,  it  was  decided  that  an  industry  standard 
microcomputer  board  along  with  its  existing  I/O  interface  circuits,  memory 
options,  and  software  development  tools  would  be  used  in  the  implementation 
of  MIBDAPS. 

Selection  of  Microcomputer  Board 

The  main  features  which  were  considered  in  the  selection  of  the  micro- 
computer board  are  given  below: 

Word  Size 

I/O  Capability 

Instruction  Set 

Memory  Capacity 

Software  Development  Tools 

Since  the  MIBDAPS  is  to  process  the  input  signals  in  real  time,  the 
choice  of  the  microcomputer  word  size  and  the  instruction  set  would  deter- 
mine the  processing  speed  of  the  system.  The  word  size  determines  the 
number  of  memory  accesses  required  for  execution  of  each  instruction.  For 
the  same  instruction  set,  a machine  with  a smaller  word  size  requires  more 
memory  access  as  compared  with  a machine  with  a larger  word  size  (Ref  4). 
Thus,  a 16-bit  microcomputer  was  an  obvious  choice.  The  signal  processing 
function  (Ref  Chap  I)  of  MIBDAPS  required  that  multiply  and  divide  instruc- 
tions should  be  a part  of  the  instruction  set.  Although  subroutines  can 
be  written  to  carry  out  multiplication/division  of  binary  numbers,  this 
technique  is  orders  of  magnitude  slower  than  hardwired  multiply/divide 
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instructions.  This  is  because  the  hardwired  instructions  are  micro- 
programmed, which  is  a lower  primitive  than  the  software  implemented 
subroutine.  Table  II  lists  the  candidate  microcomputer  boards.  Texas 
Instruments'  TMS  9900  and  the  Digital  Equipment  Corporation's  LSI-11 
have  multiply/divide  instructions  as  part  of  their  instruction  set. 

Also,  both  of  these  mucrocomputers  have  a 16-bit  word  length.  There- 
fore, the  choice  was  narrowed  down  to  these  two  boards. 

The  data  acquisition  function  of  MIBDAPS  (Ref  Chap  I)  dictated  the 
need  for  an  analog-to-digital  converter  (ADC)  board  which  would  be 
compatible  with  the  selected  microcomputer.  Also,  the  need  for  sychro- 
nization  and  timing  of  real  time  signal  processing  made  it  desirable  to 
have  a hardwired  system  to  perform  this  function.  The  Digital  Equipment 
Corporation  (DEC)  manufactures  both  an  ADC  board  and  a real  time  clock 
(RTC)  board  that  are  compatible  with  the  LSI-11.  This  feature  of  the 
LSI-11  made  it  an  obvious  choice.  Recently,  Norden  Systems  Incorporated, 
which  is  a subsidiary  of  United  Technology  Incorporated,  announced  the 
manufacturer  of  a bit-sliced  based  microcomputer  board  (designated  LSI-1 1M) 
which  is  six  times  faster  than  its  commercial  counterpart  (see  Table  II). 
This  board  also  conforms  to  military  specifications  (mil  specs).  Thus, 
the  mil  spec  LSI - 1 1 M was  selected  as  the  microcomputer  board  to  be  used 
for  MIBDAPS. 

LSI-11  Microcomputer  System 

The  LSI-11  microcomputer  system  is  a microprocessor  based  version  of 
DEC'S  PDP-11  series  minicomputer.  The  PDP-11  development  software 
(assemblers,  debuggers,  loaders,  compilers)  is  useable  on  the  LSI-11 
based  system.  H,., 
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fixed-point  arithmetic  only. 


In  the  following  paragraphs,  a brief  description  of  the  LSI-11  pro- 
cessor, interface/memory  options,  and  software  development  tools  that 
were  used  in  the  implementation  if  MIBDAPS  is  given.  These  options  were 
selected  because  they  support  either  the  hardware  or  the  software  reali- 
zation of  different  modules  of  MIBDAPS.  The  detailed  description  of  each 
of  the  above  modules  and  the  installation  instructions  are  given  in  Ref  11. 

KD11-F  Microcomputer  (Ref  11:87-99).  The  KD11-F  microcomputer  is 
contained  on  a single  8 -5-inch  by  10-inch  printed  circuit  board.  The 
module  also  contains  resident  dynamic  4K  by  16-bit  semiconductor  RAM. 

Figure  10  is  a functional  block  diagram  of  the  microcomputer  board.  The 
salient  features  of  this  microcomputer  are  given  below: 

- This  is  a 16-bit  microcomputer  board  using  four  custom  LSI  inte- 
grated circuit  chips. 

- Direct  addressing  of  32K  16-bit  words  or  64K  8-bit  bytes  (K  = 1024). 

- Efficient  processing  of  8-bit  bytes  without  the  need  to  rotate, 
swap,  or  mask. 

- Hardware  memory  stock  for  handling  structured  data,  subroutines, 
and  interrupts. 

- Direct  memory  access  (DMA)  for  high  data  rate  devices  inherent  in 
the  bus  architecture. 

- Eight  general  purpose  registers  that  are  available  for  data  stor- 
age, pointers,  and  accumulators;  two  are  dedicated— stack  pointer  (SP) 
and  program  counter  (PC). 

- A bus  structure  that  provides  position-dependent  priority  to  the 
peripheral  devices  that  are  connected  to  the  I/O  bus. 

- Vectored  interrupt  facility. 

- Typical  operating  speed  is  400  nsec  based  on  a 10  MHZ  clock. 
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Figure  10.  KD11-F  Microcomputer  Logic  Block  Diagram  (From  Ref  11:38) 
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- The  power  supply  requirements  are: 

+5V  +5%  at  2-4A  maximum  current 
+12V  +5%  at  1*10A  maximum  current 

MMVH  4K  by  16-Bit  Core  Memory  (Ref  11:101-102).  The  MMV11-A  4K  by 
16-bit  core  memory  option  provides  a nonvolatile  read/write  storage  of 
the  user's  program  and  data  (Figure  11).  This  optional  module  was  used 
so  that  system  software  which  would  be  developed  for  MIBDAPS  could  be 
stored  permanently.  This  would  eliminate  the  need  to  reload  the  program 
whenever  the  system  is  switched  ON.  The  processed  data  would  also  be 
stored  in  the  core  memory  for  either  updating  or  retrieval  at  a later 
time.  The  salient  features  of  the  MMV11-A  core  system  are: 

- 4096  by  16-bit  capacity. 

- Typical  access  time  is  equal  to  425  nsec  (475  nsec  maximum);  full 
read/restore  cycle  time  is  equal  to  1*15  nsec. 

- User-selected  bank  address;  this  allows  the  user  to  assign  any  4K 
bank  address. 

- Power  requirements  are: 

+5V  +5%  at  2-0A  maximum 
+12V  +5%  at  0-56A  maximum 

ADV11-A  Analog-to-Digital  Converter  (Ref  11:186-189).  The  ADV11-A 
is  a 12-bit  successive  approximation  analog-to-digital  converter  with  a 
built-in  multiplexer  and  sample-and-hold  for  use  on  the  LSI-11  bus  (Fig- 
ure 12).  The  multiplexer  selection  accommodates  16  single-ended  or 
8 quasi-differential  inputs. 


29 


30 


Figure  11.  NMV11-A  Logic  Block  Diagram  (From  Ref  11:105) 


Figure  12.  ADV11-A  Function  Block  Diagram  (From  Ref  11:191) 


This  A/D  converter  was  adequate  for  MIBDAPS  requirements  because 
only  two  analog  signals  needed  to  be  digitized.  Also,  It  was  required 
that  the  conversion  error  should  not  exceed  +1%  of  the  full  scale  (Ref  13). 
For  analog  signals  of  0-5V,  this  requirement  translates  to  an  error  of 
+5mv.  The  one  least  significant  bit  (LSB)  value  for  the  ADV11-A  is  2*5mv 
which  means  that  the  maximum  error  that  can  be  introduced  is  +2-5mv,  which 
is  within  the  requirement. 

The  main  features  of  ADV11-A  are: 

- Analog  input  full  scale  range  is  10 - 24V  bipolar  (-5-12V  to  +5-12V). 

- Logic  input  voltages  are: 

low  (logic  "0")  = 0-0  to  +0-7V 
high  (logic  "1")  = +2V  to  +5V 

- Resolution  of  ADV11-A  12  12-bit  binary  weighted.  The  format  of 
the  digitized  signal  is  parallel  offset  binary  right  justified. 

- The  timing  requirements  are: 

a.  External  start  should  be  a low-level  pulse,  50  nsec  minimum 
duration  and  10  ysec  maximum  duration.  Conversion  starts  on  the  leading 

edge. 

b.  Conversion  time  is  16  times  the  clock  period.  The  maximum 
clock  frequency  is  100  KHZ. 

- Power  requirements  are: 

+5V  +5%  at  2-0A  maximum  current 
+12V  +5%  at  450mA  maximum  current 

KWV11-A  Programmable  Real  Time  Clock  (Ref  11:211-212).  The  func- 
tional diagram  for  the  real  time  clock  is  shown  in  Figure  13.  This  pro- 
grairmable  clock/counter  provides  a variety  of  means  for  determining  time 
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Intervals  or  counting  levels.  It  can  be  used  to  generate  interrupts  to 
the  LSI-11  processor  at  predetermined  intervals  to  synchronize  the  pro- 
cessor to  external  events  or  to  measure  time  intervals  or  establish  pro- 
grammed ratios  between  input  and  output  events.  It  can  also  be  used  to 
start  the  ADV11A  analog-to-digital  converter  either  by  clock  counter  over- 
flow or  by  the  firing  of  a Schmitt  trigger  circuit. 

The  clock  counter  has  a resolution  of  16  bits  and  can  be  driven  from 
any  five  internal  crystal -controlled  frequencies  (100  HZ  to  1 MHZ).  It 
can  also  be  driven  from  a line  frequency  input  or  from  a Schmitt  trigger 
fired  by  an  external  clock  input.  The  performance  specifications  for  the 
KWV11-A  are: 

- Internal  Clock 

The  internal  clock  has  a frequency  accuracy  of  +0.01%.  The  base 
frequency  is  10  MHZ,  and  it  is  divided  into  five  selectable  rates  (1  MHZ, 
100  KHZ,  10  KHZ,  1 KHZ,  and  100  HZ). 

- External  Clock 

The  external  clock  can  either  be  the  line  frequency  input  or  a 
clock  generator  having  the  desired  frequency. 

- Output  Signals 

a.  A clock  overflow  signal  occurs  whenever  a certain  timing  event 
is  completed.  The  characteristics  of  this  signal  are:  500  nsec  low 
asserted  level  pulse  which  is  TTL  compatible. 

b.  Two  trigger  outputs,  Schmitt  trigger  1 and  Schmitt  trigger  2, 
are  also  available  for  either  being  used  as  a trigger  for  the  ADV11-A 
(A/D  converter  board)  or  for  any  other  external  device.  The  character- 
istics are  the  same  as  that  for  the  clock  overflow  signal. 
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DLV11  Serial  Line  Unit  (Ref  11:147-148).  The  DLV11  is  the  basic 
interface  module  used  for  connecting  asynchronons  serial  line  devices  to 
the  LSI-11  bus.  The  functional  diagram  for  the  DLV11  is  shown  in  Figure 
14.  The  salient  features  of  the  DVL11  unit  are: 

- Either  an  El A RS232  or  a 20mA  current  loop  device  can  be  inter- 
faced to  an  LSI-11  bus. 

- There  are  13  baud  rates  that  can  be  selected:  50,  75,  110,  134*5, 
150,  200,  300,  600,  1200,  1800,  2400,  4500,  and  9600. 

- The  LSI-11  bus  interface  and  control  logic  for  interrupt  processing 
and  vector  generation  are  also  available  on  the  board. 

- The  control /status  register  (CSR)  and  data  buffer  register  (DBR) 
for  programmed  or  interrupt  driven  transfer  of  data  are  available. 

Software  Development  Tools.  One  of  the  reasons  why  the  LSI-11  micro- 
computer system  was  chosen  for  mechanization  of  MIBDA°S  was  because  it  is 
software  compatible  with  the  PDP-11  computer  family.  The  following  sys- 
tems that  were  available  on  Wright-Patterson  Air  Force  Base  could,  there- 
fore, be  used  for  software  development: 

- The  Department  of  Electrical  Engineering  at  the  Air  Force  Institute 
of  Technology  (AFIT),  Wright-Patterson  Air  Force  Base,  Ohio,  has  a PDP-11/20 
along  with  a paper  tape  software  system.  (A  single  floppy  disk  drive  has 
been  added;  therefore,  the  RXV11  operating  system  could  also  be  used.) 

- The  AFAL  has  a cross-assembler  for  PDP-11  assembly  language  on  its 
DEC-10  computer  system.  Programs  could  be  created  on  the  DEC-10  system^ 
using  text  editors.  The  PDP-11  cross-assembler  (MACY11)  can  then  be  used 
to  create  absolute  loader  formatted  paper  tapes.  There  is,  however,  na 
facility  to  simulate  execution  of  these  programs. 
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Figure  14.  DLV11  Logic  Block  Diagram  (From  Ref  11:150) 
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- The  AFML  has  a PDP-11/03  with  dual  floppy  disk  drive  and  an  RXV11 
operating  system.  This  system  uses  an  LSI-11  microcomputer  and  all  inter- 
face boards.  The  RXV11  operating  system  has  a text  editor,  macro  assembler, 
linker,  FORTRAN  compiler,  and  peripheral  interchange  program  (PIP)  avail- 
able on  a floppy  disk.  This  system  was  suitable  not  only  for  developing 
and  debugging  programs,  but  also  for  real  time  program  execution. 

The  PDP-11/03  system  at  AFML  was  eventually  used  for  the  development 
of  software  for  MIBDAPS. 

Analysis  and  Design  of  MIBDAPS 

The  decisions  about  the  implementation  of  different  functional  modules 
of  MIBDAPS,  either  in  software  or  hardware,  were  guided  by  the  following 
factors: 

- Capabilities  and  limitations  of  the  KDF-11  processor  module  along 
with  4K  words  of  core  memory  used  as  the  microcomputer  system. 

- Capabilities  and  limitations  of  support  hardware  modules,  in 
particular  ADV11-A  (ADC)  and  KWV11  (RTC) . 

- Speed  of  software  execution. 

The  KDF-11  processor  along  with  4K  words  of  core  memory  (subsequently 
referred  to  as  "microcomputer  system")  would  be  used  to  implement  the  con- 
trolling and  processing  functions  of  MIBDAPS  in  software.  The  ADC  board 
(ADV11-A)  and  RTC  (KWV11)  along  with  a system  clock  circuit  (which  would 
be  designed  and  fabricated)  were  visualized  as  the  only  hardware  needed 
to  implement  the  data  acquisition,  signal  processing,  and  synchronization 
functions  of  MIBDAPS.  Keeping  these  two  factors  as  guidelines,  another 
in-depth  study  of  the  system  requirement  model  developed  in  Chapter  II 
was  carried  out.  SADT  was  employed  to  carry  out  a functional  analysis 
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of  the  system  and  a new  activity  model,  keeping  implementation  in  view, 
was  developed.  The  node  index  for  this  activity  model  is  given  in 
Table  III.  It  should  be  noted  that  this  time  the  mechanization  (M)  for 
each  node  is  also  defined  as  hardware  or  software,  or  both.  Software 
implementation  implies  programs  which  would  be  executed  on  the  micro- 
computer system.  Hardware  implementation  implies  use  of  ADC,  RTC,  or 
an  external  system  clock. 

TABLE  III 


Node  Index  for  New  Activity  Model 


A-0 

Acquire  and  Process  Data  (same  as  in  Chap  II) 

A0 

Acquire  and  Process  Data 

A1 

Control  All  Activity 

A13 

Switch  Power  to  System 

A14 

Stop  Execution 

A15 

Start  Execution 

A16 

Start  Acquisition  and  Processing  of  Data 

A162 

Initiate  Data  Acquisition 

A163 

Initiate  Signal  Processing 

A2 

Acquire  Data 

A3 

Process  Data 

A32 

Compute  Fourier  Transform 

A33 

Compute  Power  Spectrum 

A34 

Compute  Damage  Factor 

A35 

Update  CDF  Array 
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The  explanation  of  each  node  along  with  data  structures  used  and 
Input/output  signal  specifications  follows. 

Node  A0,  Acquire  and  Process  Data  (Figure  15).  The  difference 
between  node  A0  as  compared  with  node  A0*  (Figure  4)  is  the  absence  of 
module  A4  (Store  Data).  This  difference  has  come  about  because  the 
Process  Data  Module  will  perform  the  task  of  storing  and  updating  result 
that  are  stored.  The  various  activities  taking  place  in  this  node  are 
listed  below. 

- The  Control  All  Activity  module  performs  the  function  of  managing 
the  control  and  execution  of  all  other  activities.  It  also  provides  the 
DC  supply  to  all  hardware  modules. 

- The  Acquire  Data  module  continuously  samples  and  digitizes  the 
input  analog  signals,  analog  vibration  signal  and  analog  temperature 
signal.  Whenever  a set  of  data  values  have  been  digitized,  it  sets  the 
Acquisition  Done  flag. 

- The  Process  Data  module  takes  in  the  input  raw  data  and  processes 
it  to  give  an  array  of  damage  factor  versus  frequency.  This  array  is 
added  on  to  one  of  the  cumulative  data  factor  versus  frequency  arrays. 
There  are  11  CDF  versus  frequency  arrays,  one  for  each  expected  tempera- 
ture intervals  (Ref  Chap  II). 

The  input/output  and  control  signal  specifications  are: 

- Input  Signal 

a.  Analog  vibration  signal  (II):  0-5  volts. 

b.  Analog  temperature  signal  (12):  0-5  volts. 

c.  AC  power  (13):  110  volts,  60  cycles. 
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- Output  Signal 

a.  Total  Operational  Time  (01):  A display  of  operational  time 
In  hours  and  minutes. 

b.  DC  Power  to  System  (02): 

c.  Cumulative  Scale  Factor  (03):.  This  is  a 16-bit  octal  number 
which  gives  the  scale  factor  for  the  output  arrays. 

d.  CDF  Versus  Frequency  and  Temperature  Arrays  (03):  These  are 
13-word  long  arrays  having  the  data  structure  as  shown  in  Figure  16. 


CSTR1 : 

TEMPERATURE  #1 

CDSF1 

CDF ( fQ) 

V ° HZ 

CDF ( f -| ) 

f1  = 100  HZ 

- 

' 

CDF(f10) 

f = 1000  HZ 

Figure  16.  Data  Structure  for  CDF  Vs  Frequency  and  Temperature  Array 


- Control  Signal 

System  ON/OFF  Control  (Cl):  Once  the  system  is  switched  ON  and  the 
AC  power  is  available,  the  system  starts  functioning. 

Node  Ai , Control  All  Activity  (Figure  17).  This  node  performs  the 
activities  necessary  to  control  the  synchronization  of  all  other  activities 
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of  MIBDAPS.  The  mechanization  of  different  activities  in  this  node  are 


1 

either  done  by  hardware  or  software  or  a combination  of  both.  It  is 
assumed  that  the  LSI-11  computer  is  used  to  execute  the  software;  however, 
hardware  mechanization  includes  various  other  modules  which  are  being 
Interfaced  to  the  microcomputer.  The  various  activities  performed  in  this 
node  are  listed  below. 

- The  DC  power  required  by  different  modules  of  the  LSI-11  micro- 
computer system  is  provided  by  this  node.  The  110-volt  AC  supply  is  con- 
verted to  +5V  and  +12V  DC  supplies.  Whenever  the  system's  ON/OFF  switch 
is  selected  to  the  ON  position,  the  DC  supply  is  switched  ON.  This  is 
done  by  a Switch  Power  to  System  module.  This  module  also  generates  an 
interrupt  in  case  the  DC  power  is  within  an  acceptable  limit.  In  case 

of  a power  failure,  an  interrupt  is  also  generated. 

- Total  operational  time  for  the  system  is  also  computed  by  the 
Record  Operational  Time  module.  In  case  the  operational  time  reaches 

• the  selected  time  limit  (15  to  20  hours),  an  interrupt  is  generated. 

- The  Start  Execution  module,  after  receiving  the  Power  OK  inter- 
rupt, initializes  the  system  and  sends  a signal  to  Start  Acquisition 
and  Processing  of  Data. 

- The  Stop  Execution  module  stops  execution  of  the  system  func- 
tion in  case  the  following  interrupts  occur: 

Power  Fail  Interrupt 
Time  Limit  Interrupt 

- The  Start  Acquisition  and  Processing  of  Data  module  initializes 
different  variables  which  are  used  by  the  software.  It  then  signals 
nodes  A2  and  A3  to  start  acquisition  and  processing  of  data. 

Specifications  for  the  input/output  and  control  signals  are: 
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- Input  Signal 

AC  Supply  (13):  110  volts,  60  HZ. 

- Output  Signal 

a.  Total  Operational  Time  (01):  Four  4-bit  BCD  numbers. 

b.  Start  Acquisition  Control  (0202):  A software-generated  com- 
mand. The  control /status  register  of  the  ADV11-A  A/D  converter  (ADC)  is 
loaded  with  the  appropriate  value  to  enable  the  ADC.  This  initiates  the 
data  acquisition  cycle. 

c.  Start  Processing  Control  (0303):  This  is  also  a software- 
generated command,  which  activates  the  processing  modules  for  input  data 
processing. 

- Control  Signal 

System  ON/OFF  Control  (Cl):  A signal  which  initiates  the  opera- 
tion of  the  system.  A logic  "1"  value  is  available  when  the  system  is 
turned  ON  and  logic  "0"  value  when  the  system  is  turned  OFF. 

Node  A13,  Switch  Power  to  the  System  (Figure  18).  This  node  carries 
out  the  following  activities: 

- Converts  the  AC  power  supply  to  +5  volts  and  +12  volts  DC. 

- Generates  an  interrupt  to  the  microcomputer  when  DC  power  is 
available. 

- Detects  a power  failure  and  generates  a Power  Fail  interrupt  to 
the  microcomputer. 

The  input  and  output  signal  specifications  are: 

- Input  Signal 

AC  Supply  (13):  This  is  the  only  input  to  this  node.  The  system 
operates  on  11 OV  AC,  60  HZ. 
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Figure  18.  Node  A13,  Switch  Power  to  the  System 


- Output  Signal 

a.  DC  Power  Supply:  The  +5V  and  +12V  DC  power  supply  is  routed 
to  all  modules  of  MIBDAPS. 

b.  Start  Execution  (1401):  A software-generated  signal  which  is 
used  to  start  execution  of  the  program.  Bit-0  of  variable  STEXEC  is  set 
to  initiate  data  acquisition  and  processing. 

c.  Power  Fail  Interrupt  (1501);  A software-generated  signal 
which  stops  execution  of  the  program.  Bit-0  of  variable  STPXC  is  set  to 
stop  execution  of  the  program. 

Node  A14,  Stop  Execution  (Figure  19).  This  node  receives  the  Power 
Fail  interrupt;  it  then  accesses  the  vector  interrupt  address  for  the 
service  routine.  The  service  routine  is  executed  and  bit-0  of  the  vari- 
able STPXC  is  set  to  a logic  "1."  This  variable  is  sent  as  an  output 
signal  to  stop  execution  of  the  system  execution. 

Node  A15,  Start  Execution  (Figure  20).  This  node  receives  the  Power 
OK  interrupt;  it  then  accesses  the  vector  interrupt  address  for  the  serv- 
ice routine.  The  service  routine  is  executed  and  bit-0  STEXEC  switch  is 
set  to  logic  "1."  This  variable  is  sent  as  to  the  GO  Command  to  start 
system  execution. 

Node  A16,  Start  Acquisition  and  Processing  of  Data  (Figure  21).  Once 
the  system  has  been  turned  ON  and  DC  power  is  made  available  to  the  sys- 
tem, this  node  starts  execution  of  the  data  acquisition  and  data  process- 
ing functions  of  MIBDAPS. 

The  Stop  Execution  Command  (14C2)  and  Start  Execution  Command  (14C1) 
are  used  to  start  and  halt  the  execution  of  data  activity.  The  input/out- 
put signal  specifications  are: 
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Figure  19.  Node  A14,  Stop  Execution 


Figure  20.  Node  A15,  Start  Execution 


Acquisition  and  Processing  of  Data 


- Input  Signal 

a.  Acquisition  Done  Status  (0211):  Whenever  the  acquisition  of 
n-data  points  is  complete,  the  variable  ADFLG  is  set  to  logic  "1."  This 
variable  along  with  ADCNT  (number  of  data  points  to  be  converted)  and 
CDATA  (starting  address  of  data  buffer  where  data  is  to  store)  are  passed 
to  the  Initiate  Data  Acquisition  module.  All  of  these  variables  need  to 
be  initialized  whenever  data  acquisition  has  to  be  reinitiated. 

b.  Processing  Done  Status  (0312):  The  following  variables  are 
passed  to  the  Initiate  Processing  module  for  initialization:  DATA  - Con- 
tains starting  address  of  data  to  be  processed.  SFACTR,  CUMSF1 , and 
CDSFi  (i=l,  2...)  - Scale  factors  used  in  node  A3,  which  are  explained 
later;  however,  each  has  to  be  cleared  before  processing. 

- Output  Signal 

a.  Start  Acquisition  (0201):  This  signal  is  generated  when  the 
control/status  register  (ADCSR)  for  the  A/D  converter  is  loaded  with  the 
appropriate  bit  pattern.  Once  the  ADCSR  is  loaded,  the  acquisition  of 
data  starts, 

b.  Start  Processing  (0302):  This  is  a software-generated  com- 
mand, and  it  initiates  the  processing  of  data. 

- Control  Signal 

a.  Start  Execution  (14C1):  This  is  a software  switch,  STEXEC, 
which  is  set  to  logic  "1"  whenever  DC  power  is  available  and  the  system 
has  been  initialized. 

b.  Stop  Execution  (14C2):  This  is  also  a software  switch,  STPXC, 
which  is  set  to  logic  "1"  whenever  a power  failure  occurs. 

Node  162,  Initiate  Data  Acquisition  (Figure  22).  This  node  shows 
explicitly  the  functional  decomposition  of  the  module  Initiate  Data 
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Acquisition  of  node  A16.  The  STEXEC  switch  (161  Cl)  and  STPXC  switch 
(161C2)  control  the  activity  in  this  node.  The  input  signal  0211— CDATA, 
ADCNT,  and  ADFLG— is  initialized  by  modules  1,  2,  and  3,  respectively. 

The  Generate  Start  Acquisition  Command  loads  the  ADCSR  with  appropriate 
bits  to  initiate  acquisition  of  data  by  ADC,  only  after  all  of  the  above 
mentioned  switches  are  initialized  and  STEXEC  is  set  to  logic  "1." 

Node  163,  Initiate  Signal  Processing  (Figure  23).  This  node  shows 
the  functional  decomposition  of  the  Initiate  Signal  Processing  module  of 
node  A16.  The  activities  in  this  node  are  similar  to  node  A162.  The 
input  signal  Processing  Done  Status  (0321)  and  output  signal  Start  Pro- 
cessing (0301  ) are  different.  The  input  includes  SFACTR,  CUMSF1 , CDSFi 
SPCTRA,  and  DATA  switches.  These  are  initialized  in  this  node.  The 
Start  Processing  Command  generated  in  this  node  starts  the  execution  of 
the  processing  of  acquired  data. 

Node  A2,  Acquire  Data  (Figure  24).  This  node  performs  the  activity 
of  acquiring  data  and  storing  it  in  temporary  storage.  Once  the  acqui- 
sition of  N-data  points  is  complete,  the  acquisition  is  stopped  and  the 
ADFLG  switch  is  set  to  logic  "1."  When  the  analog  temperature  signal  is 
selected  for  acquisition,  the  temperature  interval  flag,  TFLG,  is  loaded 
with  a number  from  0 to  11.  The  temperature  interval  number  that  is 
loaded  in  TFLG  depends  on  the  value  of  the  temperature  signal  at  the 
time  of  conversion.  The  input  and  output  signal  specifications  are: 

- Input  Signal 

a.  Analog  vibration  signal:  0-5V. 

b.  Analog  temperature  signal:  0-5V. 
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Start  Data 
Acquisition 


Figure  24.  Node  A2,  Acquire  Data 


- Output  Signal 


/ 


a.  Temperature  Interval  Flag  (0301):  A decimal  number  between 
0 and  11  is  loaded  in  TFLG.  This  number  depends  on  the  value  of  the 
temperature  signal  at  the  time  of  acquisition. 

b.  N-Input  Data  Points  (0302):  The  data  structure  for  the  data 
buffer  storage  is  given  in  Figure  25.  CDATA  is  the  switch  which  contains 
the  starting  address  of  the  data  buffer  to  be  used.  ISTR1  and  ISTR2  are 
the  starting  addresses  of  the  two  buffers.  Each  buffer  is  N-words  long. 

c.  Acquisition  Done  Flag  (0103):  Bit-0  of  switch  ADFLG  is  set 
to  logic  "1"  whenever  data  acquisition  of  N-points  is  complete. 

Node  A3,  Process  Data  (Figure  26).  The  activities  performed  by  this 
will  be  easier  to  understand  if  the  flow  of  data  through  each  module  is 
explained.  The  Input  Data  (0211)  is  N-points  of  acquired  data  which  are 
stored  in  a buffer  area.  The  starting  address  of  this  buffer  is  available 
in  switch  DATA.  The  fourier  transform  of  these  N-points  is  computed,  and 
the  results  are  stored  back  in  the  same  buffer  area.  Switch  SFACTR  con- 
tains the  scale  factor  for  the  resulting  array.  The  fourier  transformed 
array  serves  as  an  input  to  the  Computer  Power  Spectrum  module.  This 
module  computes  the  power  spectrum  for  first  M- frequency  components 
(including  DC  term).  The  output  for  this  module  is  M-element  power  spec- 
trum array  and  switch  CUMSF1 , which  has  the  scale  factor  for  this  array. 
The  M-element  DF  array  is  added  to  one  of  the  existing  CDF  versus  fre- 
quency arrays.  The  value  of  switch  TFLG  is  used  to  select  the  final 
array,  where  the  DF  array  is  added.  The  data  structure  associated  with 
each  module  of  this  node  is  explained  in  the  subordinate  nodes. 

Node  A32,  Computer  Fourier  Transform  (Figure  27).  This  node  takes 
the  input  data  from  the  buffer  storage,  computes  the  fourier  transform 
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Figure  27.  Node  A32,  Computer  Fourier  Transform 


w ^ 

of  the  N-data  points,  and  stores  the  results  in  the  original  data  buffer 
array.  The  input/output  data  structure  is  shown  in  Figure  28.  Switch 
DATA  contains  the  starting  address  for  the  input/output  data  buffer  area. 

Each  data  buffer  area  is  N-words  long. 

Node  A33,  Computer  Power  Spectrum  (Figure  29).  This  node  takes  the 
N-point  transform  array  and  computes  the  power  spectrum  for  first  M-fre- 
quency  components.  The  resultant  M-point  array  is  stored  in  another 
temporary  data  storage  area  whose  starting  address  is  stored  in  switch 
SPCTRA.  The  input/output  data  structure  is  shown  in  Figure  30. 

Node  A34,  Computer  Damage  Factor  (Figure  31).  The  M-point  power 
spectrum  array  is  input  to  this  node  (3311).  It  takes  the  starting 
address  from  the  SPCTRA  switch  and  computes  the  damage  factor  for  each 
M-data  points  in  the  input  array.  The  resultant  DF  array  is  passed  on 
to  Node  A35. 

Node  A35,  Update  CDF  Array  (Figure  32).  This  node  takes  the  damage 
factors  computed  by  node  A34  and  updates  one  of  the  11  CDF  versus  fre- 
quency arrays.  Switch  TFLG  (02C1)  is  used  to  select  the  starting  address 
of  the  output  array  which  corresponds  to  the  correct  temperature  interval. 

The  output  data  structure  for  CDF  versus  frequency  arrays  is  given  in 
Figure  33. 


Summary 

In  this  chapter,  the  selection  of  a microcomputer  system  to  imple- 
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ment  MIBDAPS  was -matte.'  Based  on  the  selected  microcomputer  system,  the 
partitioning  of  MIBDAPS  into  hardware  or  software  realizable  modules  was 
performed.  A detailed  design  of  MIBDAPS,  including  mechanization  of  dif- 
ferent modules,  input/output  data  structures,  and  breakdown  of  different 
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modules  into  primitives,  was  carried  out  using  SADT. 
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Figure  28.  Input/Output  Data  Storage  for  Node  A32 
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Figure  29.  Node  A33,  Computer  Power  Spectrum 


Figure  31.  Node  A34,  Computer  Damage  Factor 


IV.  Design  and  Implementation  of  Software 


The  general  design  of  MIBDAPS  discussed  in  Chapter  III  presented 
the  foundation  for  system  implementation.  The  division  of  the  overall 
system  into  software  and  hardward  modules  was  also  accomplished.  In  this 
chapter,  the  algorithm  for  each  software  module  is  developed.  Subsequently, 
the  LSI-11  assembly  language  code  is  written  for  the  designed  algorithms.  The 
source  listing  for  each  software  module  is  included  in  Appendices  E,  F,  and  G. 
Before  the  software  algorithms  are  designed,  a mathematical  analysis  of 
the  signal  processing  function  of  MIBDAPS  is  presented.  This  analysis  is 
performed  to  determine  system  parameters  and  their  variation  during  the 
signal  processing  of  data. 

Determination  of  System  Parameters 

The  development  of  any  software  module  requires  a priori  knowledge 
of  the  computation  taking  place  within  the  module  and  the  relationship  of 
input  and  output  parameters.  Therefore,  a mathematical  analysis  of  the 
data  acquisition  and  signal  processing  functions  of  MIBDAPS  is  required. 

This  analysis  should  result  in  the  type  and  amount  of  computation  required 
within  each  software  module  of  MIBDAPS  and  the  different  input/output 
parameters.  Because  the  processing  of  the  acquired  data  is  to  be  done 
in  real  time  (Ref  Chap  I),  timing  constraints  placed  on  the  system  would 
also  be  determined  in  this  analysis. 

Data  Acquisition  Analysis.  The  MIBDAPS  is  required  to  digitize  the 
analog  vibration  signal  and  the  analog  temperature  before  any  signal  pro- 
cessing is  carried  out.  The  acquisition  analysis  for  the  analog  vibration 
signal  is  carried  out  first.  This  would  be  followed  by  the  acquisition 
analysis  for  the  analog  temperature  signal. 
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The  analcg-to-digital  conversion  of  a continuous  signal  is  a two- 
step  progress,  as  shown  in  Figure  34.  The  analog  signal  is  first  sam- 
pled, and  then  followed  by  quantization  of  the  amplitudes  of  the  sampled 
signal  (Ref  2:155). 


Figure  34.  Sampling  and  Quantization  of  an  Analog  (From  Ref  2:155) 

The  limitations  and  applicability  of  these  two  steps  will  first  be 
considered  and,  particularly,  the  consequences  arising  from  the  selection 
of  a given  sampling  rate  and  number  of  quantizing  levels. 

The  first  problem  that  arises  from  sampling  a continuous  signal  is 
that  of  aliasing.  The  nature  of  this  problem  is  illustrated  in  Figure  35 
which  shows  that  the  same  set  of  sampled  data  points  can  describe  a num- 
ber of  time  series  histories,  which  are  indistinguishable  to  the  digital 
computer. 
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Figure  35.  Aliasing  (From  Ref  2:157) 
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To  overcome  this  problem,  the  continuous  signal  which  has  a finite 
band  width  up  to  and  including  B rad/sec  should  be  sampled  at  ir/B  sec 
intervals  (Ref  2:160).  Therefore,  the  sampling  interval  AT  is  given  as: 


AT  = tt/B  (2) 

B = 27rfm  (3) 

where 

fm  = maximum  frequency  content  of  the  signal 
From  Eq  (2)  and  Eq  (3)  we  get 


AT 


'2fm 


(4) 


The  sampling  frequency  F,  which  is  the  reciprocal  of  AT,  is  given  by 


F = 


(5) 


Equation  (5)  defines  the  sampling  frequency  required  to  reduce  the  effect 
of  aliasing.  This  rate  is  also  known  as  the  Nyquist  sampling  rate  (Ref  6 
225).  The  analog  vibration  signal  input  to  MIBDAPS  is  expected  to  have 
a finite  band  width  extending  from  OHZ  to  1000  HZ  (Ref  Chap  I);  there- 
fore, Eq  (5)  specifies  that  the  sampling  frequency  for  this  signal 
should  be  equal  to  or  greater  than  2000  HZ. 

The  second  step  in  data  acquisition  is  quantization.  "Quantization" 
Is  defined  as  the  representation  of  a variable  amplitude  series  of  Dis- 
crete sample  values  as  an  equivalent  series  of  discrete  numbers  represent 
ing  their  amplitude  values..  This  process  can  only  be  an  approximation 
because  the  number  of  bits  in  a digital  representation  is  limited,  even 
though  the  analog  signal  can  assume  an  infinite  number  of  values.  The 
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numerical  values  of  the  quantized  variable  may  be  represented  by  some 
form  of  binary  code  to  permit  entry  into  the  digital  computer  (Ref  2:162). 
Figure  36  shows  the  transfer  characteristics  of  the  quantizer  which  gives 
a 3-bit  code  for  an  input  analog  signal.  The  analog  values  are  quantized 
by  partitioning  the  continuum  into  eight  discrete  ranges.  All  analog 
values  within  a given  range  are  represented  by  the  same  digital  code,  which 
corresponds  to  the  nominal  mid-range  value.  There  is,  therefore,  an 
inherent  quantization  uncertainty  of  +}$  LSB.  The  only  sure  way  to  reduce 
this  quantization  uncertainty  is  to  increase  the  number  of  bits  for  the 
representation  of  output  and  have  more  quantization  levels  (Ref  8:69). 


Figure  36.  Quantizer  Transfer  Characteristics  (From  Ref  8:69) 
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The  ADC  used  for  data  acquisition  in  MIBDAPS  is  a 12-bit  successive 
approximation  converter  which  has  one  LSB  value  equal  to  25mv  (Ref  Chap 
III).  Thus  the  quantization  error,  which  is  +%  LSB,  is  equal  to  +12-5mv 
The  analog  vibration  and  temperature  signal  will  have  a full-scale  value 
of  5 volts;  therefore,  the  quantization  error  will  be  +0.25  percent  of 
the  full-scale  value. 

Signal  Processing  Analysis.  The  signal  processing  function  of 
MIBDAPS  has  been  broken  down  into  three  sub-functions  (Ref  Chap  III, 

Node  A0).  The  mathematical  analysis  of  the  signal  processing  in  each 
subunit  will  now  be  carried  out. 

1 . Processing  an  Analog  Vibration  Signal.  The  power  spectral 
density  of  the  digitized  vibration  signal  needs  to  be  computed  first 
before  the  damage  factor  for  each  frequency  of  interest  is  computed 
(Ref  Chap  III).  This  can  be  accomplished  by  first  determining  the 
frequency  domain  spectrum  of  the  input  time  domain  signal.  The  Fourier 
transform  is  a method  that  converts  a signal  from  time  to  frequency 
domain.  Since  the  input  signal  is  a sampled  version  of  the  continuous 
time  domain  signal,  the  discrete  fourier  transform  (DFT)  method  was 
used  for  the  transformation.  Equation  (6)  gives  the  expression  for  the 
conversion  of  N sampled  data  points  from  time  domain  to  frequency  domain 

N-l 

F(n)  = Z x(k)exp(--j^k))  n = 0,  1 , 2 ...  N-l  (6) 

k=0  N 

where  F(n)  is  the  nth  point  in  the  frequency  domain  and  n(k)  the  kth 
point  in  the  time  domain.  The  x(k)  may  be  complex  and  the  F(n)'s  are 
always  complex.  We  can  rewrite  Eq  (6)  as 


70 


(7) 


F(n)  = 


N-l  ■, 

I x(k)Wnk,  n = 0,  1,  2 
k = 0 


N-l 


where 


W = exp(^21) 


(8) 


Figure  37  shows  the  required  frequency  domain  spectrum  (Ref  Chap  I) 
from  0 HZ  to  1000  HZ.  The  frequency  interval  Af,  which  is  100  HZ,  will 
be  used  to  determine  the  number  of  time  domain  samples  required  for  the 
Fourier  transformation.  Equation  (9)  gives  the  relationship  between  Af 
and  the  sampling  interval  T (Ref  3:87). 


> 

II 

-4^ 

(9) 

or 

T = V 

(10) 

'Af 

also 

T = MxAT 

(ID 

therefore 

(12) 

Equation  (12)  gives  the  number  of  sampled  points  M required  to  get  a fre- 
quency spectrum  shown  in  Figure  37.  The  relationship  between  the  signal 
in  the  frequency  domain  and  time  domain  is  shown  in  Figures  37  and  38, 
respectively.  From  Eq  (4)  AT  is  computed  equal  to  0.5  msec,  and  Eq  (10) 
gives  T equal  to  10  msec.  Therefore,  the  number  of  sampled  points  M 
required  to  get  the  desired  frequency  spectrum  is  20,  which  is  calculated 
using  Eq  (12).  Once  the  frequency  spectrum  of  the  time  domain  signal  is 
computed,  the  next  step  in  signal  processing  is  to  compute  the  power  spec- 
tral density  for  each  frequency  component.  After  the  Fourier  transforma- 
tion, each  frequency  component  F(n)  in  Eq  (6)  is  a complex  number  as  given 
in  Eq  (13): 
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F(")  = xn  + jy„ 

the  power  spectral  density  can  be  calculated  using  Eq  (14): 

A(n)  = (xn)2  + (yn)2  (14) 

However,  the  nature  of  the  input  signal  is  random  (Ref  13)  and  to  get 
a good  representation  of  the  power  spectral  density  a mean  value  for 
the  P set  of  spectral  density  points  is  taken  (Ref  2:195).  This  is 
expressed  mathematically  in  Eq  (15): 

A(n)  = l s"1  A(n) , n = 0,  1,  ...  M-l  (15) 

p k=0  K 

where  subscript  k denotes  a unique  set  of  M spectral  density  values,  and 
A is  the  mean  power  spectrual  density. 

Once  the  mean  power  spectral  density  is  computed,  the  damage  factor 
at  each  frequency  component  can  be  evaluated  using  Eq  (16)  (Ref  13): 

DF(n)  = , n = 0,  1 ...  M-l  (16) 

(n) 

where 

J[n)  = / A(n)  (17) 

Here  fn  is  the  frequency  at  which  the  damage  factor  is  being  calculated. 
From  Eqs  (16) and  (17)  the  following  relationship  is  derived: 


(18) 


(19) 


1 
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Equation  (19)  eliminates  the  need  to  compute  the  root  mean  square  (RMS) 
of  the  average  power  spectral  density.  The  CDF,  which  is  a function  of 
frequency  and  temperature,  is  obtained  by  continuous  summation  of  DF 
values  over  a predetermined  interval.  Equation  (20)  gives  the  required 
relationship: 

R 

CDF(t. , m)  = ? DF.(m),  m = 0,  1 ...  M-l  (20) 

1 J—  l J 

Here,  t.  represents  the  ith  temperature  interval  for  which  the  CDF  is  com- 
puted, and  subscript  j indicates  the  number  of  DF  value  sets  that  have 
been  summed  together. 

2.  Analog  Temperature  Signal.  The  temperature  signal  is  sampled 
every  time  a set  of  20  data  points  of  the  vibration  signal  is  processed 
and  the  power  spectral  density  has  been  computed.  Since  P set  of  power 
spectral  density  values  are  summed  [Eq  (15)Q,  there  are  P temperature 
values  that  can  be  summed  up  to  give  an  average  temperature  value  T.  The 
average  temperature  value  would  be  available  when  the  DF  is  computed  and 
the  final  output  array  (CDF  versus  frequency  and  temperature)  is  to  be 
updated.  The  following  equation  gives  the  mathematical  relationship  for 
T (average  temperature): 

1 P"1 

T = / E T(k)  (21) 

P k=0 

This  concludes  the  mathematical  analysis  of  the  data  acquisition  and 
signal  processing  functions  of  MIBDAPS.  The  software  modules  to  imple- 
ment different  functional  activities  will  now  be  developed. 
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The  software  modules  which  will  be  used  to  perform  different 
functional  activities  are  developed  in  subsequent  paragraphs. 


Fourier  Transform  Module.  The  Fourier  transform  is  a method  to  con- 
vert the  signal  spectrum  from  the  time  domain  to  a spectrum  in  the  fre- 
quency domain.  Since  MIBDAPS  suses  digitized  data  for  processing,  the  DFT 
algorithm  can  be  used  to  transform  the  digitized  vibration  signal  from 
the  time  domain  to  the  frequency  domain.  Equation  (7)  gives  the  mathe- 
matical relationship  to  convert  a digitized  signal  to  a frequency  domain. 

A computer  program  to  perform  the  summation  given  in  Eq  (7)  can  be  written, 
but  it  turns  out  that  this  can  be  a very  slow  and  time-consuming  process 
if  the  equation  is  implemented  as  given  (Ref  6:267).  The  reason  for  this 
inefficiency  is  that  the  transformation  of  N-data  points  requires  com- 


plex multiplications  and  N(N-l)  complex  additions  (Ref  3:151).  To  speed 
up  the  DFT  implementation  on  the  computer,  Cooley  and  Tukey  (Ref  3:151) 
developed  a fast  Fourier  transform  (FFT)  algorithm,  which  for  N = 2X 
points  reduces  the  complex  multiplication  to  Nx/2  and  the  complex  addi- 
tion to  Nx  (Ref  3:151 ) . 

Appendix  C gives  the  theoretical  development  of  the  base  2 FFT  algo- 
rithms. It  was  determined  earlier  that  20  points  of  DFT  are  required  to 
get  the  required  frequency  spectrum  of  the  vibration  signal.  N = 20  does 
not  satisfy  the  relationship  N = 2X  where  x is  an  integer;  therefore,  a 
base  2 FFT  algorithm  would  not  be  used.  FFT  algorithms  for  arbitrary 
factors  can  also  be  written  (Ref  3:184).  An  FFT  algorithm  for  20-point 
DFT  was  written,  but  it  was  found  that  125  complex  multiplications  and 
108  complex  additions  are  required.  On  the  other  hand,  a 32-point  FFT 
algorithm  requires  80  complex  multiplications  and  160  complex  additions. 
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so  It  was  decided  to  use  a 32-point  FFT  algorithm.  This  would  give  a 
frequency  spectrum  from  0 HZ  to  1500  HZ,  and  by  picking  up  the  first  11 
frequency  components,  the  desired  spectrum  for  MIBDAPS  could  be  obtained. 
Figure  39  shows  the  relationship  between  the  time  domain  signal  and  the 
resultant  frequency  domain  signal.  The  sampling  frequency  F should  now  be 
3200  HZ  [Ref  Eq  (5)],  and  the  sampling  period  AT  becomes  equal  to  0*3125  msec. 


Figure  39.  Time  and  Frequency  Domain  Signals  for  32  Samples 


The  FFT  algorithm  developed  in  Appendix  C can  handle  both  real  and 
complex  input  data.  However,  the  input  to  MIBDAPS  is  real  data,  and  com- 
putation time  would  be  halved  if  2N  points  of  real  data  are  treated  as  N 
points  of  complex  data.  This  would  enable  the  use  of  an  N-point  FFT  algo- 
rithm. The  theoretical  background  for  transforming  2N  points  of  real  data 
using  an  N-point  FFT  is  given  in  Appendix  D;  this  algorithm  was  eventually 
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used  for  implementation  of  the  Fourier  transform  module.  The  descrip- 
tion of  the  algorithm  is  given  in  the  following  paragraphs. 

1.  Algorithm  to  Implement  Fourier  Transform  Module.  Figure  40 
shows  the  flow  chart  of  the  algorithm  that  was  implemented.  The  input 
data  array  is  assumed  to  have  32  real  data  points  in  contiguous  memory 
locations.  Since  the  FFT  algorithm  described  in  Appendix  D would  be 
implemented,  it  is  assumed  that  the  input  array  has  16  complex  data 
points  (all  odd  points  are  treated  as  real  points  and  all  even  points  are 
treated  as  complex  points).  The  bit-inversion  subroutine  is  called  first 
so  that  the  data  elements  are  bit- reversed  before  being  processed  by  the 
FFT  subroutine.  This  is  done  to  ensure  that  the  output  of  the  FFT  sub- 
routine would  be  in  the  correct  order.  The  FFT  computation  is  performed 
on  the  bit-inversed  array;  the  resultant  array  would  be  16  complex  Fourier 
transformed  points.  To  get  the  correct  transformed  array  for  real  data 
points,  the  post-processing  computation  (Ref  Appendix  D)  is  performed; 
this  would  give  the  16  complex  values  corresponding  to  frequencies  from 

0 HZ  to  1600  HZ.  Although  a 32-point  real  data  array  should  result  in 
32  complex  point  Fourier  transformed  arrays,  the  algorithm  being  imple- 
mented only  gives  16  complex  Fourier  transformed  points.  This  discrepancy 
is  because  the  other  16  complex  transformed  points  are  the  complex  con- 
jugates of  these  points  and  are  obtainable  easily.  However,  these  points 
are  not  required  for  further  processing  in  MIBDAPS  and,  therefore,  are 
not  computed. 

2.  Description  of  Different  Subroutines. 

a.  Bit-Inversion  Subroutine.  The  flow  diagram  for  the  bit- 
inversion  (BIV)  subroutine  is  given  in  Figure  41.  The  starting,  address 
of  the  input  data  a^ray  is  stored  in  pointer  DATA.  The  five  LSB  of  the 
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address  for  any  data  elements  are  used  to  compute  the  bit- inverted 
address.  These  two  memory  location  contents  are  then  swapped.  This 
process  is  repeated  for  all  data  elements,  and  the  resulting  array  is  in 
bit-reversed  order. 

b.  FFT  Subroutine.  The  transformation  of  input  data  is  carried 
out  in  four  steps  (or  passes).  The  manipulation  of  data  within  this  sub- 
routine follows  the  scheme  given  in  the  signal  flow  graph  shown  in  Figure 
C-2,  Appendix  C.  During  the  first  pass,  there  are  no  complex  multiplica- 
tions (Ref  Table  C-I,  Appendix  C);  therefore,  the  following  equations 
need  to  be  implemented: 


R(m)'  = R(m)  + R(n) 

(22) 

I(m)'  = I (m)  + I (n ) 

(23) 

R(n)' = R(m)  - R(n) 

(24) 

I (n ) ' = I(m)  - I(n) 

(25) 

The  letter  within  the  parenthesis  signifies  a pair  of 

points  which  are 

used  in  the  calculation. 

The  second,  third,  and  fourth  pass  require  complex  multiplications 

n k 

because  of  the  factor.  The  equations  implemented 

in  these  passes  are 

X(m)  = X(m)  + X(n)Wy 

(26) 

X(n)  = X(m)  + X(n)Wz 

(27) 

Here 

X(m)  = R(m)  + jl(m) 

X(n)  = R(n)  + j I ( n ) 

Wy  = exp(— -j-'-H)  = WZ 
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Also,  the  exponents  of  W,  y and  z,  are  N/2  apart,  such  that  the  follow- 
ing relationship  holds: 


wy  = -WN/2  +y  = -WZ  (28) 

Equations  (26)  and  (27)  can  be  rewritten  using  Eq  (28). 

X(m)  = X(m)  + X(n)wy  (29) 

X(n)  = X(m)  - X(n)wy  (30) 

or 

R(m)  + jl(m)  = R(m)  + jl(m)  + [R(n)  + jI(n)](Cosy  + jSiny)  (31) 

. R(n)  + jl(n)  = R(m)  + jl(m)  - [R(n)  + jI(n)](Cosy  + jSiny)  (32) 


Collecting  terms  and  rearranging  Eqs  (31)  and  (32),  the  following 


relationship  evolves: 

R(m)'  = R(n)Cosy  + I(n)Siny  + R(m)  (33) 
R(n)'  = -R(n)Cosy  - I(n)Siny  + R(m)  (34) 
I(m)'  = -R(n)Siny  + I(n)Cosy  + I(m)  (35) 
I(n)'  = R(n)Siny  - I(n)Cosy  + I(m)  (36) 


Equations  (33),  (34),  (35),  and  (36)  are  finally  used  in  computation 
during  the  second,  third,  and  fourth  pass.  The  flow  diagram  for  the  FFT 
subroutine  is  given  in  Figure  42.  The  flow  diagram  shows  a call  to  sub- 
routine SCALE;  this  subroutine  is  used  to  ensure  that  the  data  array  is 
scaled  to  avoid  any  arithmetic  overflow  during  computations  in  any  pass. 
The  flow  diagram  for  the  subroutine  SCALE  is  shown  in  Figure  43.  Before 
the  subroutine  is  called,  it  is  assumed  that  the  starting  address  is 
loaded  in  Register  RO,  and  the  number  of  points  in  the  array  is  stored  in 
Register  R1 . 
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c.  Post-Processing  Subroutine.  After  the  FFT  of  input  data  is 


completed,  the  post-processing  of  the  resulting  array  is  carried  out  to 
get  the  correct  output  array  (Ref  Appendix  D).  Equation  (D-8),  Appendix  D, 
is  rewritten  below: 

X(n)  = Js[V(n)  + V*(N-n)]  + VjCV(n)  - V*(N-n)]W2Nn  (37) 

where 

V(n)  = Y(n)  + jZ(n)  (38) 


Substituting  Eq  (38)  in  Eq  (37)  and  separating  the  real  and  imaginary 
parts,  the  following  relationships  are  derived. 


Re[X(n)]  = Js[Y(N-n)]  + Cos(^-)[Z(n)  + Z(N-n)] 

- Sin(™)[Y(n)  - Y(N-n)]  (39) 

Im[X(n)]  = %[Z(n)  - Z(N-n)]  - Cos(^)[Y(n)-Y(N-n)] 

- Sin(^-)[Z(n)  + Z(N-n)]  (40) 

If  Rp,  Ip,  Rn,  and  Im  are  defined  as  follows 


Rp  = Y(n)  + Y(N-n)  (41) 
Ip  = Z(n)  + Z(N-n)  (42) 
Rm  = Y(n)  - Y(N-n)  (43) 
Im  = Z(n)  - Z(N-n)  (44) 


and  using  the  fact  that  -'  Sin = Sin(^)  and  Cos = -Cos(^), 
Eqs  (39)  and  (40)  can  be  rewritten: 


Re[X(n)]  = Rp  + Ip  Cos(M)  - Rm  Sin(M)  (45) 

Im[X(n)]  = Im  - Ip  Sin(HJI)  - Rm  Cos(^)  (46) 
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(47) 

(48) 


Re[X(N-n)]  = Rp  - Ip  Cos(^)  - Rm  Sin(^) 

Im[X(N-n)]  = Im  - Ip  Sin(^)  + Rm  Cos(^p) 

N N 

Equations  (45),  (46),  (47),  and  (48)  were  eventually  used  in  the 
implementation  of  the  post-processing  subroutine.  The  flow  diagram  for 
the  post-processing  subroutine  is  given  in  Figure  44. 

3.  Assembly  Language  Code  for  FFT.  The  FFT  algorithm  was  coded  in 
LSI-11  assembly  language.  A deviation  from  the  flow  diagram  of  Figure  42 
was  made  and  all  other  subroutines  (SCALE,  BIT  INVERSION,  POST-PROCESSING) 
were  included  in  the  FFT  subroutine.  Appendix  E gives  a source  listing  of 
FFT  subroutines.  The  complete  analysis  of  the  subroutine  execution  time 
is  also  given  in  Appendix  E. 

Power  Spectrum  Module.  The  Fourier  transformed  array,  which  is  the 
output  from  the  Fourier  transform  module,  contains  16  complex  frequency 
components.  The  power  spectrum  for  the  first  11  frequency  components 
(from  0 HZ  to  1000  HZ)  needs  to  be  computed.  This  is  accomplished  using 
a subroutine  called  PRSPEC  (Power  Spectrum).  The  flow  diagram  from  the 
subroutine  PRSPEC  is  shown  in  Figure  45,  which  shows  that  the  first 
11  frequency  components  of  the  power  spectrum  are  calculated  by 
squaring  the  real  and  imaginary  parts  and  summing  them  together.  The 
output  is  added  to  the  previous  elements  of  an  array  whose  starting 
address  is  stored  at  pointer  SPCTRA.  This  feature  helps  in  computing 
the  average  power  spectrum  which  would  then  be  used  for  computing  the 
damage  factor  at  each  frequency  of  interest  [Ref  Eq  (15)].  The  assembly 
language  code  for  PRSPEC  is  given  in  Appendix  F;  this  appendix  also  con- 
tains the  analysis  of  execution  time  for  this  subroutine. 
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Figure  44.  Flow  Diagram  for  Post-Processing  Subroutine 


Damage  Factor  Module.  The  subroutine  to  calculate  the  DF  at  each 


frequency  of  interest  implements  Equation  (19),  which  is  rewritten: 


DF(n)  = n = 0,  1,  ...  10 

f(n)5’4 


DF(n)  = K'CAtn)]1'6 


where 


‘ ■ V„)5-4 

I 

Since  the  factor  K is  constant  for  each  frequency,  the  function  imple 
mented  in  this  subroutine  is  [A(n)]^  To  implement  this  exponential 
factor,  the  expected  values  of  A(n),  which  are  going  to  be  less  than 
equal  to  1*0,  were  divided  into  10  intervals  as  shown  in  Figure  46. 


(A)1*6 


■1  -2  -3  -4  -5  -6  *7  -8  .9  1*0 


Figure  46.  Damaie  Factor  Versus  A Curve  Divided  Into  10  Intervals 


A polynomial  curve  fit  program  PLSCF,  available  as  an  AFIT  Library  sub- 
routine on  the  CDC  6600  computer  (Ref  1 5: A-5) , was  used  to  fit  poly- 
nomial curves  to  each  of  the  10  intervals  shown  in  Figure  46.  It  was 
experimentally  determined  that  if  the  ordinate  (DF)  was  divided  into  8 
intervals,  shown  in  Figure  47,  a second-order  polynomial  gave  the  least 
RMS  error  for  these  intervals  as  compared  with  RMS  errors  for  other  com- 
binations of  intervals.  Appendix  K gives  a description  of  the  curve 
fitting  program.  The  equation  for  a second-order  polynomial  and  the 
coefficients  for  the  8 intervals  (Figure  47)  are  also  given  in  Appendix  K. 


Figure  47.  Damage  Factor  Versus  A Curve  with  the  Selected  Intervals 
for  Curve  Fitting  Program  (8  Intervals) 

The  flow  diagram  for  the  subroutine  to  compute  the  damage  factor  is 
given  in  Figure  48;  this  subroutine  is  called  POLYN.  Note  that  the 
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resulting  DF  values  are  added  to  the  existing  elements  of  the  selected 
output  array  (Ref  Chap  III).  This  is  done  to  get  the  CDF  as  a function 
of  temperature  arid  frequency.  The  assembly  language  code  for  subroutine 
POLYN  is  given  in  Appendix  G.  This  appendix  also  gives  an  analysis  of 
execution  time  for  this  subroutine. 

Control  Executive  (EXEC).  The  Control  Executive  is  the  systems  pro- 
gram that  handles  the  execution  of  all  the  modules  in  MIBDAPS.  The  con- 
trol and  synchronization  of  different  activities  is  also  handled  by  this 
program. 

The  state  transition  diagram  that  shows  the  initial  switching  on 
sequence  is  given  in  Figure  49.  Once  the  power  is  switched  ON  to  the 
system,  the  microcomputer  executes  its  own  interval  microcode  and  jumps 
to  location  (17300)g  where  if  finds  a WAIT  instruction.  Once  the  POWER  OK 
or  POWER  FAIL  interrupt  occurs,  the  service  routine  for  either  of  these 
interrupts  is  executed.  In  case  the  POWER  ON  sequence  is  successful,  the 
execution  of  the  system  starts. 

Figure  50  shows  the  flow  diagram  for  the  continuous  execution  of  the 
system.  The  acquisition  of  data  is  done  using  the  interrupt  feature  of 
ADV11-A.  Therefore,  the  data  acquisition  occurs  while  the  processing  of 
previously  acquired  data  is  going  on.  The  real  time  clock  would  be  used 
for  synchronization  for  the  entire  data  acquisition  and  processing  cycle, 
but  this  feature  has  not  yet  been  included  in  the  flow  diagram.  The 
"dashed"  block  in  Figure  50  (page  94)  would  accomplish  the  synchroniza- 
tion of  different  processing  cycles. 


The  assembly  language  code  for  EXEC  is  given  in  Appendix  H.  The 
execution  time  analysis  is  also  included  in  Appendix  H. 
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Figure  49i  State  Transition  Diagram  for  Initial  Power  ON  Sequence 
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ADCLOP:  Counter  that  tells 

the  system  which  data 
buffer  to  use  for 
acquisition  of  data 
(Ref  Chap  III). 

PRCLOP:  This  counter  allows 
summation  of  10  sets 
of  power  spectrum 
computations  to  obtain 
the  average  spectrum 
[Ref  Eq  (18)]. 
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Figure  50.  Flow  Diagram  for  Continuous  Execution 
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Sumnary 

In  this  chapter,  the  mathematical  analysis  of  data  acquisition 
and  processing  functions  of  MIBDAPS  was  presented.  The  theoretical 
background  necessary  for  the  development  of  each  software  module  was 
discussed.  This  was  followed  by  the  generation  of  a flow  diagram  and 
assembly  language  code  for  each  software  module.  The  design  of  hard- 
ware realizable  modules  of  MIBDAPS  is  undertaken  in  the  next  chapter. 
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V.  Design  and  Implementation  of  Hardware 


Introduction 

The  design  and  implementation  of  software  realizable  modules  of 
MIBDAPS  was  presented  in  Chapter  IV.  The  next  step  in  the  development 
of  MIBDAPS  was  the  design  and  implementation  of  hardware  realizable 
modules  (Ref  Chap  III).  In  this  chapter,  the  general  design  approach 
is  discussed  and  the  implementation  of  the  designed  circuits  using  stand- 
ard off-the-shelf  integrated  circuits  is  explained.  At  the  end  of  this 
chapter,  the  DC  power  supply  requirements  are  discussed  and  the  use  of 
standard  DC  power  supply  modules  is  suggested.  The  author's  conception 
of  the  prototype  module  of  MIBDAPS  along  with  its  physical  dimensions  is 
also  presented. 

Design  of  Hardware  Circuits. 

The  design  of  different  hardware  realizable  modules  is  described  in 
the  following  paragraphs.  The  design  is  based  on  the  functional  require- 
ments spelled  out  in  Chapter  III. 

System  Clock  Generation  Circuit.  The  data  acquisition  function  of 
MIBDAPS  required  that  the  vibration  analog  signal  be  sampled  at  a fre- 
quency of  3200  HZ  (Ref  Chap  IV).  Therefore,  to  gather  data  at  this  sam- 
pling rate,  it  was  necessary  to  trigger  the  ADC  (ADV11-A)  by  a 3200  HZ 
external  clock  (Ref  11:335-338).  The  same  clock  was  also  required  to 
synchronize  different  processing  loops  during  system  execution  (Ref  Chap 
IV).  Figure  51  is  the  schematic  that  shows  the  hardware  configuration 
where  the  system  clock  was  used.  It  should  be  noted  that  the  output  of 
the  clock  is  fed  to  the  RTC  before  being  applied  to  the  ADC.  This  has 
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been  done  to  ensure  that  the  clock  is  conditioned  by  the  Schmitt  trigger 
circuit  on  the  RTC  board  (Ref  11:211-212)  before  being  applied  to  the 
ADC. 


Figure  51.  Hardware  Configuration  Showing  System  Clock  Generator 

To  design  a stable  3200  HZ  clock,  SN74LS424  clock  generator  IC  and  a 
64K  HZ  crystal  were  used  as  the  reference  oscillator  (the  choice  of  crystal 
frequency  was  arbitrary).  This  reference  frequency  was  divided  down  to 
get  the  required  3200  HZ  clock.  Two  SN7490  ICs  were  configured  to  accom- 
plish the  appropriate  frequency  division.  Figure  52  shows  the  circuit 
diagram  of  the  designed  clock  generator. 

POWER  ON  Interrupt  Circuit.  The  LSI-11  bus  foundation  module 
(DRV11-P),  which  is  manufactured  by  DEC,  can  be  used  for  user  defined 
interfaces  (Ref  11:169).  This  module  is  supplied  with  the  logic  necessary 
for  interfacing  to  the  LSI-11  bus.  This  logic  includes  bus  transceivers, 
a device  address  comparator,  protocol  logic,  interrupt  logic,  vector 
address  comparator,  and  bus  receivers  and  invertors  (Ref  11:169).  The 


device  interrupt  generator  and  vectors  address  generation  feature  of 
this  module  were  used  in  the  design  of  the  POWER  ON  interrupt  circuit. 

The  POWER  OK  interrupt  occurs  whenever  the  DC  power  to  the  system 
has  been  switched  ON  and  +5V,  +12V  DC  supplies  are  within  acceptable 
limits.  The  POWER  ON  transition  and  a signal  from  the  power  sense  cir- 
cuit would  be  used  to  generate  the  requisite  sequence  of  signals  neces- 
sary to  activate  the  interrupt  and  vector  address  generator  feature  on 
the  DRV11-P.  (The  reader  is  referred  to  Ref  11:177  for  an  understanding 
of  the  logic  requirement  placed  on  the  requesting  device.) 

The  block  diagram  of  the  designed  circuit  is  given  in  Figure  53. 


Figure  53.  Block  Diagram  of  PWROK  Interrupt  Generator 


This  circuit  assumes  that  the  +5V  supply  is  available  to  itself  and  the 
DRV11-P  module  before  the  POWER  ON  switch  (SW1)  makes  the  OFF-to-ON 
transition.  This  assumption  would  create  a problem  if  there  is  only  one 
power  ON/OFF  switch  in  the  system.  The  solution  to  overcome  this  problem 
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Is  to  have  a RUN/STANDBY  switch,  which  is  put  to  the  STANDBY  position 
before  the  POWER  ON/ORR  switch  is  switched  ON.  The  RUN/STANDBY  switch 
Is  then  switched  to  the  RUN  position,  which  would  generate  the  appropriate 
signal  sequence  for  the  interrupt  and  interrupt  vector  generation  logic 
on  the  DRV11-P.  Figure  54  shows  the  designed  circuit  diagram  along  with 
the  timing  diagram  of  the  logic  signals  generated. 

POWER  FAIL  Interrupt  Generator.  This  circuit  would  be  used  to  gen- 
erate an  interrupt  whenever  a power  failure  occurs  or  when  the  RUN/STANDBY 
switch  makes  a RUN-to-STANDBY  transition.  Figure  55  shows  the  block  dia- 
gram of  the  hardware  configuration  required.  This  circuit  was  not  designed 
due  to  the  shortage  of  time.  The  DRV11-P  could  be  used  to  generate  the 
POWER  FAIL  interrupt;  however,  a logic  circuit  to  initiate  the  interrupt 
request  and  interrupt  vector  cycles  would  have  to  be  designed. 


Figure  55.  POWER  FAIL  Interrupt  Generation  Logic 

Operational  Time  Recorder.  The  block  diagram  for  the  operational 
time  recorder  is  given  in  Figure  56.  The  system  clock  is  fed  to  this 
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circuit,  and  it  keeps  track  of  the  time  for  which  the  system  remains 
operational.  Whenever  the  operational  time  reaches  the  15-hour  limit 
an  interrupt  is  generated  and  the  system  executives  a "halt"  cycle. 
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The  operational  time  is  also  displayed  on  the  front  panel.  In  case 
either  a power  failure  occurs  or  the  system  is  switched  to  STANDBY  mode, 
the  current  operational  time  would  be  stored  in  the  core  memory.  This  is 
done  to  obtain  the  total  operational  time  for  which  the  system  monitors  a 
test  piece  in  case  the  system  was  operated  in  different  time  segments. 
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Figure  56.  Operational  Time  Recorder 


Power  Supply  Requirement 

The  power  requirement  for  each  hardware  module  used  in  MIBDAPS  is 
tabulated  in  Table  IV.  The  +5V  DC  supply  is  required  to  source  a maxi- 
mum current  of  16- 63A,  and  the  +12V  is  required  to  source  2 - 21 A maximum. 

A comparative  analysis  of  the  +5V  and  +12V  DC  supplies  available  from  dif- 
ferent manufacturers  was  done.  The  following  characteristics  were  compared 
during  this  analysis:  maximum  output  current,  load  regulation,  ripple  and 
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noise,  case  size,  and  weight.  Based  on  this  analysis,  model  numbers 
HE  252  and  HE  212,  the  +5V  DC  and  +12V  DC  supplies,  manufactured  by 
Computer  Products  Incorporated,  are  the  recommended  DC  power  supplies 
for  MIBDAPS.  Table  V lists  the  salient  features  of  these  recommended 
supplies;  it  should  be  noted  that  the  maximum  rated  current  for  the  sup- 
plies is  more  than  the  required  maximum  current  for  MIBDAPS  (see  Table  IV). 
This  would  enable  future  additions  of  the  required  hardware  modules. 


Mechanical  Description  of  MIBDAPS 

The  author's  conception  of  the  final  shape  of  MIBDAPS  (called 
MIBDAPS-I)  is  given  in  Figure  57.  The  detailed  diagram  showing  the  layout 
of  different  subunits  and  their  dimensions  is  given  in  Figure  58.  It 
should  be  noted  that  the  rear  view  shows  an  AC  inlet  socket  and  an  RS-23Z 
outlet.  The  RS-232  would  be  provided  to  connect  the  system  to  a TTY, 
which  would  be  used  as  an  I/O  device  to  get  a printout  of  CDF  versus 
frequency  and  temperature  arrays. 


Surranan 


In  this  chapter,  the  design  and  implementation  of  two  hardware  modules 
was  presented  along  with  the  block  diagram  descriptions  of  the  other  twa 
hardware  modules  which  need  to  be  designed.  The  power  supply  requirements 
was  determined  and  a recommendation  for  off-the-shelf  DC  power  supplies 
was  made. 

The  physical  dimensions  for  the  mainframe  and  the  location  of  sub- 
units  within  the  mainframe  were  made.  Diagrams  depicting  the  physical 
chape  and  dimensions  of  MIBDAPS  have  also  been  drawn.  Integration  of  the 
system  was  accomplished  next.  The  results  of  testing  the  integrated  sys- 
tem are  discussed  in  the  following  chapter. 
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TABLE  IV 


Power  Requirement  for  MIBDAPS 


Description 
of  Modules 

Current  Required 

From  +5V  Supply 
(nominal /maxi  mum) 

Current  Required 

From  +12V  Suoply 
(nominal /maximum) 

LSI -1 1 M 

Processor 

1-8A/2-4A 

0-8A/1 • 1A 

DLV11  Serial 

Line  Unit 

1 *6A/1 *6A 

0-18A/0-25A 

ADV11-A  A/D 
Convertor 

2-0A/2-0A 

0-45A/0-45A 

MMV11-A  4K 

Core  Memory 

3-0A/7-0A 

0-6A/0-6A 

REV11-C  Boot- 
strap ROM  and 
Terminator 

1-0A/1-88A 

— 

KWV11-A  Real 

Time  Clock 

1-75A/1-75A 

0-01  A/0 -01 A 

Total  current  requirement  for  +5V  supply  = 
Total  current  requirement  for  +12V  supply 

11  - 1 5A/1 6 * 63A. 

= 2-04A/2-21A. 

TABLE  V 

Recommended  Power  Supply  Characteristics 


Description 

Output 

Current 

(maximum) 

Load  Reg. 
Maximum 
(NL-FL) 

Ripple 

and 

Noise 

Case  Size 
(in.) 

L x W x H 

Weight 

(lbs.) 

+5V 

HE  252 
(1 10V  AC 
output) 

20 -0A 

+0-1% 

50MV  P-P 
(13  MV 
RMS) 

6*5x4*5x 

3-25 

+12V 

HE  212 
(110V  AC 
output) 

3-OA 

+0-1% 

20MV  P-P 
(2MV  RMS) 

6-5x4-5x 

1 -8 

1-7 
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VI . Testing  and  Experimental  Results 

Introduction 

The  testing  of  the  software  during  the  development  stage,  the  inte- 
gration of  hardware  and  software,  and  the  testing  of  the  integrated  sys- 
tem are  discussed  in  this  chapter.  The  experimental  setup  used  in  both 
the  software  development  stage  and  the  system  integration  stage  is 
explained,  and  the  results  of  the  experimentation  and  software  developed 
for  testing  are  presented  at  the  end  of  the  chapter. 

Software  Debugging  and  Testing 

The  software  algorithms  that  were  designed  and  coded  in  Chapter  IV 
were  debugged  and  tested  on  a PDP-11 /03  minicomputer  system.  The  config- 
uration of  the  development  system  used  is  given  in  Figure  59.  A three- 
step  approach  was  taken  to  develop  software.  First,  the  assembly  language 
code  was  assembled  and  syntax  errors  were  removed.  Second,  the  logical 
errors  detected  during  program  execution  were  debugged.  And,  lastly,  the 
testing  of  software  for  computational  accuracy  was  done.  The  flow  diagram 
given  in  Figure  60  shows  the  complete  software  development  cycle. 


Figure  59.  Configuration  of  the  Software  Development  System 
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Figure  60.  Software  Development  Cycle 
109 


Testing  of  Software  Modules.  After  the  syntax/semantics  and  logical 
error  had  been  removed,  the  testing  of  each  software  module  was  under- 
taken. The  approach  followed  was  to  write  FORTRAN  simulated  programs 
which  duplicated  the  computations  performed  by  the  software  module  being 
tested.  These  simulation  programs  were  executed  on  the  CDC  6600.  The 
input  data  to  the  simulation  programs  was  similar  in  format  and  value  as 
that  used  during  the  execution  of  the  actual  software  module.  The 
results  obtained  from  the  simulation  were  compared  with  earlier  results 
from  the  system  software  and,  in  case  of  discrepancies,  the  system  soft- 
ware code  was  checked  for  computational  errors.  Appendix  J gives  the 
simulation  program  used  to  check  the  subroutines  FOURIE,  PRSPEC,  and 
POLYN  (Ref  Chap  IV).  The  test  data  used  and  the  results  obtained  from 
the  simulation  programs  and  system  programs  are  also  tabulated  for  com- 
parison (Ref  Appendix  J) . These  results  are  discussed  at  the  end  of 
this  chapter. 

Integration  of  Software  and  Hardware 

The  system  integration  was  done  once  the  software  modules  had  been 
developed  and  tested.  The  system  integration  phase  required  setting  up 

a bench  model  for  MIBDAPS  and  testing  the  control  execution.  The  timing 

« 

constraints  for  both  data  acquisition  and  processing  were  also  determined 
at  this  stage. 

Bench  Model  for  MIBDAPS.  The  configuration  of  the  bench  model  for 
MIBDAPS  is  given  in  Figure  61.  The  function  generator  was  used  to  simu- 
late the  vibration  signal,  and  the  square  wave  generator  was  used  as  the 
system  clock.  The  bench  model  was  tested  by  running  the  control  execu- 
tive program  which  exercised  the  data  acquisition  and  processing  function 
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of  MIBDAPS.  To  check  the  results  at  different  stages  of  acquisition  and 
processing,  a decimal  print-out  subroutine  called  DECPR  was  included;  this 
subroutine  was  called  by  a macro  called  DECMAC.  The  source  listing  for 
subroutines  DECPR  and  macro  DECMAC  is  included  in  Appendix  I.  The  final 
source  listing  for  the  control  executive,  which  is  given  in  Appendix  H 
does  not  contain  this  print-out  feature  because  during  the  actual  opera- 
tion of  MIBDAPS  the  intermediate  results  are  not  required  to  be  printed. 

The  results  of  the  control  executive  operation  are  discussed  later  in  the 
chapter. 

Results  of  the  Experimentation 

The  result  obtained  from  testing  individual  software  modules  and  the 
operation  of  the  control  executive  are  presented  in  the  following  sections. 

Software  Test  Results.  The  results  of  software  testing  are  presented 
as  the  root  means  square  (RMS)  error  between  the  output  from  the  simula- 
tion program  and  the  output  from  the  actual  software  module.  These 
results  are  given  in  the  tables  in  Appendix  J.  The  RMS  error  criterion 
was  used  because  an  estimate  of  the  average  error  for  each  software 
module  was  desired.  The  mean  of  average  RMS  error  is  .0387  (see  Appen- 
dix J).  This  error  is  acceptable  because  it  is  much  less  than  the 
specified  tolerance  of  +10%  (Ref  13). 

Execution  Results  for  Control  Executive.  The  two  aspects  that  are 
important  from  the  functional  standpoint  of  MIBDAPS  are  the  execution 
time  for  the  acquisition  and  processing  loop  and  the  accuracy  of  time 
computations. 

The  execution  times  for  each  software  module  were  computed  sepa- 
rately (Ref  Chap  IV),  and  the  total  execution  time  for  the  control 
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executive  was  computed  using  the  flow  diagram  given  in  Figure  62.  The 
loop  that  computes  the  power  spectrum  for  a new  set  of  data  points  takes 
29-391  msec,  and  the  time  taken  to  compute  the  damage  factor  once  the 
average  power  spectrum  has  been  computed  (Ref  Chap  IV)  is  equal  to  4-596 
msec.  Figure  63  shows  the  timing  relationship  between  the  actual  acqui- 
sition and  processing  time  and  the  desired  time  as  defined  in  the  initial 
requirement  (Ref  Chap  I).  The  designed  system,  therefore,  would  update 
the  CDF  versus  frequency  and  temperature  array  every  298-48  msec.  For  a 
new  set  of  32  data  points,  the  power  spectrum  computation  takes  29-36  msec; 
the  average  over  10  such  intervals  would  yield  a good  estimate  of  the 
power  spectral  density  of  the  input  vibration  signal.  This  estimate  would 
not  be  in  error  as  compared  with  the  average  power  spectral  density  obtained 
if  the  power  spectrum  computed  had  taken  less  than  10  msec  (as  originally 
stipulated).  This  statement  is  true  because  the  input  vibration  signal 
can  be  assumed  to  be  wide-sense  stationary  (Ref  13),  and  the  estimate  of 
the  average  power  spectral  density  of  such  a signal  depends  on  the  number 
of  power  spectrum  samples  that  are  averaged  and  not  on  the  rate  at  which 
these  samples  are  available  (Ref  6). 

Interpreting  the  Output  Array.  The  CDF  versus  frequency  and  tempera- 
ture array  would  contain  fractional  binary  numbers  that  have  to  be  corre- 
lated  with  actual  physical  data.  An  analysis  of  the  scaling  introduced 
during  data  acquisition  and  processing  of  data  is  now  discussed. 

Assume  that  an  accelerometer  is  used  as  the  vibration  signal  sensor 

I 

which  gives  out  a voltage  of  K1  for  every  32  ft/sec^  of  acceleration. 

(This  acceleration  is  referred  to  as  one  "g.")  Therefore,  the  input 
voltage  from  the  vibration  signal  is  given  by  the  equation: 
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Figure  63.  Timing  Relationship  Between  Actual  and  Desired  Execution  Times 


e(k)  = K^k) 


(51) 


Here,  G(k)  is  the  acceleration  encountered  by  the  test  piece. 

The  analog-to-digital  convertor,  ADV11-A,  gives  a 12-bit  binary 
representation  of  the  input  voltage.  Since  the  computation  carried  out 
In  MIBDAPS  assumes  fractional  numbers  (Ref  Chap  IV),  an  input  voltage  of 
one  volt  would  be  represented  by  the  fraction  K2  which  is  defined  below. 

(007777).  x 2'15) 

u o / co\ 


Therefore,  each  data  point  x(n)  is  given  by  the  following  equation: 

x(k)  = K2e(k)  (53) 

The  Fourier  transformation  of  these  data  points  yields  the  equiva- 
lent value  in  the  frequency  domain.  Using  Eq  (6  ),  we  get 


F(n)  = Z x(k)W  , n = 0,  1,  ...  N-l  (54 

K=0 

since  the  Fourier  transformation  module  (subroutine  FOURIE)  scales  the 
data  array  to  avoid  memory  overflow  (Ref  Chap  IV).  Therefore,  the 
resultant  Fourier  transformed  array  is  given  by  the  following: 


Z(n)  = (SFACTR)F(n) 


Here,  SFACTR  is  the  value  by  which  the  Fourier  transformed  array  has  been 
scaled.  (This  value  in  subroutine  FOURIE  has  been  fixed  equal  to  3.) 

The  power  spectrum  subroutine  PRSPEC  introduces  another  scale 
factor  called  CUMSF1  (Ref  Chap  IV).  Also,  since  an  average  power  spec- 
trum density  is  computed,  Eq  (15)  can  be  modified  to  give: 
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A(n)  = K3A(n) 

(56) 

where 

K3  = CUMSFl/p 

(57) 

and 

A(n)  = |Z(n)  |2 

(58) 

The  cumulative  damage  factor  computation  introduced  a scale  factor 
called  CDSFi  (here  i = 1,  2,  ...11  ),  and  Eq  (19)  can  be  modified  and 
rewritten  as: 


CDF(n)  = 

K4[a']1,6 

(59) 

where 

K4  = 

CDSFi 

(60) 

Lf (n)]5'4 

The  equations  derived  here  can  be  used  to  interpret  the  output  array  CDF 
versus  frequency  and  temperature. 

Summary 

In  this  chapter,  the  testing  procedures  used  to  check  the  software 
modules  as  well  as  the  integrated  system  were  discussed.  The  experi- 
mental results  obtained  for  individual  software  modules  and  also  for- 
the  overall  system  were  presented.  The  scale  factors  introduced  by  each 
software  module  and  the  equations  necessary  to  interpret  the  output 
array  were  also  derived- 


VII.  Recommendations  and  Conclusions 
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Introduction 

At  present,  a bench  model  of  MIBDAPS  exists;  the  development  of  a 
prototype  model  will  require  additional  work.  This  effort  would  include 
developing  software  to  test  the  system's  accuracy,  configuring  the  hard- 
ware in  a mainframe,  and  incorporating  the  system  error  detection  and 
recovery  features.  This  chapter  presents  recommendations  for  future 
development  and  conclusions  of  the  investigation. 

Recommendations 

The  recommendations  for  future  work  are  listed  in  the  following 
paragraphs: 

- Software  to  test  the  performance  and  accuracy  of  results  should 
be  developed.  This  software  should  include  provisions  to  generate  input 
test  signals  which  duplicate  real  world  signals.  The  software  should 
gather  memory  overflow  statistics  for  varying  operating  intervals. 

These  statistics  would  enable  the  determination  of  the  fail-safe  operat- 
ing period. 

- The  software  to  output  the  CDF  versus  frequency  and  temperature 
arrays,  at  the  conclusion  of  the  system  operation,  needs  to  be  developed. 
This  software  should  either  print  the  output  arrays  in  a tabular  form  or 
plot  CDF  versus  frequency  and  temperature  contours.  (The  availability 
of  a plotter  is  essential  for  the  latter.) 

- At  present,  the  control  executive  includes  a latent  code  whlcfr 
needs  to  be  activated  in  the  future.  This  latent  code  has  been  written 
to  handle  the  acquisition  of  a temperature  signal.  It  would  also 
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compute  the  average  temperature  during  the  interval  when  the  power 
spectral  density  is  being  computed. 

- The  ADC  used  in  the  bench  model  can  handle  16  single-ended  in- 
puts. At  present,  only  two  inputs  are  being  used  for  acquiring  the  analog 
temperature  and  analog  vibration  signals.  Therefore,  similar  data  from 
seven  additional  test  pieces  can  be  acquired  and  processed.  The  possi- 
bility of  increasing  the  number  of  test  pieces  monitored  by  MIBDAPS  should 
be  explored.  Changes  in  existing  control  executive  would  be  required  for 
handling  the  additional  data. 

- The  need  to  minimize  the  existing  MIBDAPS  configuration  is  vital 
to  those  applications  where  size  and  weight  of  MIBDAPS  is  critical  (Ref 
13).  In  case  this  venue  is  explored  for  developing  the  prototype,  then 
the  following  changes  are  suggested  to  the  existing  system: 

a.  Use  of  static  RAM  in  place  of  core  memory  would  decrease  the 
weight  and  power  requirement. 

b.  Use  of  custom-designed  dual  channel  ADC  would  eliminate  the 
need  for  the  DEC  manufactured  ADC  board. 

c.  Design  of  a dedicated  FFT  hardwired  module  would  decrease  the- 
memory  requirements  for  the  software.  The  FFT  integrated  circuit  chip 
R5601*,  manufactured  by  Reticon,  would  be  a candidate  device.  If 
designed,  the  FFT  module  should  be  bus-compatible  with  the  LSI-11W 
microcomputer  and  act  as  a DMA  device. 

- The  DC  power  sense  circuit  which  would  sense  DC  supply  variations 
needs  to  be  designed.  At  switch  ON,  if  the  DC  supply  is  within  limits, 
this  circuit  would  send  a signal  to  the  POWER  OK  interrupt  generator. 

♦Ref  Preliminary  Data  Sheet:  Quad  Chipped  Transversal  Fitter, 

27  March  1978. 
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In  the  event  of  power  failure,  this  module  would  signal  the  POWER 
FAIL  interrupt  generator  circuit  which,  in  turn,  would  interrupt  the 
processor  (Ref  Chap  IV). 

- An  operational  time  recorder  is  required  to  keep  track  of  the 
time  for  which  the  system  remains  in  operation.  There  should  be 
provisions  to  move  the  current  value  in  the  buffer  register  (Figure  56) 
to  memory.  This  would  enable  the  permanent  record  of  total  operational 
time.  There  should  also  be  provisions  to  move  the  previously  stored 
operational  time  back  to  the  buffer  register  and  recording  the  addi- 
tional operational  time.  This  would  enable  the  system  to  record  the 
aggregate  operational  time  even  if  the  system  is  required  to  monitor 
the  same  test  piece  at  different  intervals. 

- The  DC  power  backup  for  use  in  case  of  a main  power  failure 
should  also  be  designed  for  the  prototype.  The  backup  power  should 
maintain  the  system  for  the  duration  it  takes  to  execute  the  power 
fail  service  routine. 

Conclusions 

The  bench  model  of  MIBDAPS  as  implemented  is  the  first  step  towards 
developing  a prototype  system.  The  initial  design  and  implementation  of 
different  functions  of  MIBDAPS  were  successfully  carried  out.  The 
choice  of  the  LSI-11M  microcomputer  as  the  processor  for  MIBDAPS  and 
the  development  of  software  to  acquire  and  process  data  fulfilled  the 
initial  requirements  placed  on  MIBDAPS  (Ref  Chap  II). 

The  time  constraints  placed  on  the  processing  function  of  MIBDAPS 
were  exceeded  because  of  the  software  implementation  of  the  FFT  module. 
The  individual  software  modules  have  been  tested  and  found  accurate 
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within  the  bounds  of  a 16-bit  machine.  Therefore,  the  testing  of  the 
integrated  system  would  be  a fairly  easy  task.  The  incorporation  of 
the  recommendation  discussed  earlier  in  this  chapter  would  result  in 
a prototype  model.  This  prototype  system  would  then  accomplish  the 
acquisition  and  processing  of  requisite  data  for  damping  the  layer 
breakout  design  as  envisioned  by  AFFDL. 

The  computation  accuracy  of  the  data  acquisition  and  processing 
functions  of  the  integrated  system  (MIBDAPS  bench  model)  could  not  be 
determined  because  of  the  shortage  of  time.  But  because  the  performance 
of  the  individual  software  modules  was  acceptable,  it  is  a fair  assump- 
tion that  the  integrated  system  would  perform  accurately,  too. 

The  designed  system  in  its  present  form  shows  the  practicability 
of  using  the  state-of-the-art  technology  to  accomplish  the  task  of  data 
acquisition  and  processing. 

The  application  of  such  a system  is  not  limited  to  monitoring  air- 
craft hardware  for  vibration  and  temperature  data  only.  It  can  be  con- 
figured to  act  as  an  environmental  recorder  where  additional  analog 
signals  need  to  be  monitored  and  processed.  This  could  be  accomplished 
by  making  appropriate  changes  in  the  software.  In  its  present  form, 
the  system  is  capable  of  computing  the  FFT  and  power  spectrum  of  a 
signal;  therefore,  this  feature  can  be  used  for  spectrum  analysis  func- 
tions. This  investigation  has  successfully  demonstrated  the  feasibility 
of  developing  a microcomputer  based  system  from  the  conceptual  phase  to 
a working  bench-type  model. 
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APPENDIX  A 


Structured  Analysis  and  Design  Technique 

This  appendix  (Ref  9)  gives  a short  description  of  how  structured 
analysis  (SA)  models  are  constructed  and  explains  the  SA  diagram  con- 
ventions used  in  this  paper.  It  must  be  noted  that  the  format  used  to 
present  the  models  in  this  paper  is  not  standard  according  to  the  rules 
developed  by  SofTech.  The  changes  were  made  to  present  the  models  in  a 
manner  which  is  more  familiar  to  readers  who  have  no  experience  with  SA 
models.  Although  the  format  is  not  that  used  by  SofTech,  the  diagrams 
of  the  models  are  organized  and  related  according  to  SofTech  procedures, 
and  the  conventions  used  to  construct  individual  diagrams  are  standard. 

The  structured  analysis  design  technique  (SADT)  is  a general  purpose 
top-down,  modular  technique  for  modeling  functions.  The  functions  may  be 
as  varied  as  farming  or  manufacturing,  but  SA  was  developed  primarily  as 
a software  requirements  definition  and  design  tool.  Although  a complete 
SA  model  actually  consists  of  two  models,  one  for  activities  and  one  for 
data,  this  paper  employs  only  activity  models  so  the  conventions  described 
here  are  those  which  apply  to  activity  models. 

An  SA  activity  model  consists  of  a series  of  diagrams  which  present 
in  progressively  more  detail  the  activities  necessary  to  perform  some 
function.  Each  diagram  represents  a self-contained  activity  which  is 
part  of  the  overall  function.  A diagram  shows  how  its  activity  is  decom- 
posed into  subactivities  and  how  the  subactivities  are  related  to  each 
other.  The  subactivities  in  each  diagram  may  then  be  decomposed  on  sep- 
arate diagrams  which  leads  to  a tree  structure  of  several  levels.  At  the 
top  is  one  diagram  which  represents  the  whole  function,  and  at  the  bottom 
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Figure  A-1  shows  how  an  SA  model  would  appear  if  all  the  diagrams 
were  on  one  page.  Of  course,  in  real  SA  diagrams  only  one  level  of  decom- 
position is  shown,  but  the  figure  demonstrates  the  top-down  nature  of  SA 


and  the  way  activities  are  grouped  into  modules.  In  the  figure,  as  in 
real  models,  one  large  box  represents  the  whole  function,  and  that  is 
decomposed  into  successive  levels  of  related  activities.  The  decomposi- 
tion process  continues  until  the  desired  amount  of  detail  has  been  devel- 
oped, which  may  require  more  levels  than  shown  in  Figure  A-1.  Another 
thing  to  note  is  that  while  the  figure  shows  only  three  subactivities  in 
each  decomposition,  any  number  from  three  to  six  is  acceptable. 

From  Figure  A-1,  it  should  be  apparent  that  SA  diagrams  are  con- 
structed with  boxes  and  arrows.  In  an  activity  model,  each  box  represents 
an  activity  and  is  called  a "node."  Arrows  represent  data  where  the  word 
"data"  is  used  in  a very  general  sense  to  include  anything  that  is  not  an 
activity.  Figure  A-2  shows  the  different  meanings  given  to  arrows  depend- 
ing on  which  side  of  a box  they  enter  or  leave.  An  input  is  data  that  is 
modified  by  the  activity  to  produce  an  output.  A control  is  data  which 
may  or  may  not  be  converted  into  output,  but  which  in  some  way  restricts 
the  activity  (starts  or  stops  it,  for  example).  Every  box  must  have  at 
least  one  control  arrow.  A mechanism  is  a person  or  thing  which  acts  as  a 
processor.  Mechanism  arrows  are  often  omitted  when  the  processor  is  the 
same  for  all  nodes.  No  limit  is  placed  on  the  number  of  arrows  which  may 
interface  with  a side  of  a box,  but  it  is  common  practice  to  group  related 
types  of  data. 

Between  boxes,  arrows  may  split  and  join.  In  general,  all  branches 
of  an  arrow  contain  the  same  data  unless  a branch  is  given  a separate 
label.  This  convention  is  summarized  in  Figure  A-3,  which  also  gives  two 
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forms  of  OR  branches.  The  OR  branches  are  used  to  show  that  data 
follows  one  path  or  the  other,  but  not  both. 

When  two  nodes  are  related  so  that  the  output  of  each  is  a control 
for  the  other,  a special  two-way  arrow  may  be  used.  Figure  A-4  shows  a 
mutual  control  situation  with  a two-way  arrow  and  the  equivalent  form 
with  normal  allows.  An  arrow  showing  mutual  control  has  two  labels  sepa- 
rated by  a slash;  the  first  label  identifies  data  going  forward,  and  the 
second  is  the  feedback  data. 

A special  numbering  system  is  used  to  distinguish  between  nodes  at 
different  levels  and  between  nodes  at  the  same  level.  In  an  activity 
model,  node  numbers  are  prefixed  with  the  letter  "A."  For  preliminary 
nodes,  A is  follows  by  a dash  and  a number.  Node  A-l  may  be  used  to  show 
the  model  in  relationship  to  other  functions  (Figure  3).  Node  A-0  serves 
as  a cover  sheet  for  the  model;  the  node  is  simply  a box  showing  inputs, 
outputs,  controls,  and  mechanisms  for  the  function  which  the  model  is  to 
describe  (Figure  4).  Decomposition  begins  in  node  A0.  Note  in  Figure  5 
that  each  box  of  the  decomposition  is  numbered;  the  boxes  on  all  decom- 
position diagrams  are  numbered,  and  this  number  is  used  to  form  the  node 
number.  For  the  activities  subordinated  to  node  A0,  the  node  number  is 
simply  the  box  number  on  A0,  Process  Data  in  Figure  15,  for  example, 
becomes  node  A3.  From  this  level  on,  the  node  number  is  a combination  of 
the  node  number  of  the  parent  diagram  and  the  box  number  of  the  subordi- 
nate. As  an  example,  the  decomposition  of  Process  Data  is  given  in  Fig- 
ure 26.  Box  2,  Compute  Fourier  Transform,  is  assigned  the  node  number  A32. 
Subordinates  of  A32,  if  diagrammed,  would  have  the  numbers  A321,  A322,  and 
so  on  through  the  last  box  number. 
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Figure  A-l . Top-down  View  of  an  SA  Model 


(Ref  9) 


Figure  A-4.  Arrows  Showing  Mutual  Control 


(Ref  9) 


A special  code  called  an  I COM  code  (input,  control,  output,  mechan- 
ism) is  used  to  identify  arrows.  The  code  contains  a number,  a letter. 


i 


and  another  number  (Figure  A-5). 


• • -nnbXc 

^7\ 

number 

letter  number 

•••nn  = node  identification 

c = number  to  differen- 

X = type  of  signal 

tiate  between  arrows 
of  same  type 

b = box  identifier 

■ 

Figure  A-5.  ICOM  Code 


The  first  number  has  two  parts:  the  digit  immediately  preceeding  the 
letter  is  that  of  source  or  destination  box,  and  the  digits  preceeding 
this  number  identify  the  node  on  which  the  box  is  located.  The  letter 
refers  to  the  type  of  arrow:  I for  input,  C for  control,  0 for  output, 
and  M for  mechanism.  The  last  number  distinguishes  between  arrows  of  the 
same  type  on  a box.  One  exception  to  the  rule  is  node  A0  where  the  ICOM 
code  has  no  number  preceeding  the  letter;  this  has  been  done  to  differen- 
tiate between  signals  within  the  system  and  signals  which  are  input,  out- 
put, or  control  from  the  outside  world.  If  any  signal  that  appears  on 
node  A0  is  also  present  on  a subordinate  node,  the  identifying  ICOM  code 
remains  unchanged. 

A few  important  points  about  the  text  describing  each  SA  diagram 
must  be  included  here.  The  text  is  intended  to  point  out  the  highlights 
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of  a diagram  and  not  repeat  all  the  details.  As  an  aid  to  following 
the  discussion,  the  ICOM  code  is  included  in  parenthesis  following  any 
reference  to  specific  diagram  features.  Finally,  the  text  describing 
the  diagrams  in  this  paper  includes  I/O  specifications,  which  are  not  a 
standard  SA  practice. 


APPENDIX  B 


Instruction  Timing 

LSI-11M  Instruction  Execution  Time  (Ref  10) 

The  execution  time  for  an  instruction  depends  on  the  instruction 
itself,  the  modes  of  addressing  used,  and  the  type  of  memory  referenced. 
In  most  cases,  the  instruction  execution  time  is  the  sum  of  a basic  time, 
a source  address  (SRC)  time,  and  a destination  address  (DST)  time. 

Instr  Time  = Basic  Time  + SRC  Time  + DST  Time 
(Basic  Time  = fetch  Time  + Decode  Time  + Execute  Time) 

Some  of  the  instructions  require  only  some  of  these  times.  All  timing 
information  is  in  microseconds,  unless  otherwise  noted.  Times  are  typi- 
cal; processing  timing  can  vary  +20%. 

The  following  talbes  give  the  basic  time  for  the  LSI-11M  instruc- 
tion set. 
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TABLE  B-I 


Source  and  Destination  Times 


Mode 

SRC  Time 
(Word) 

SRC  Time 
(Byte) 

DST  Time 
(Word) 

DST  Time 
(Byte) 

0 

0 

0 

0 

0 

1 

1 .05ps 

1 . 05yS 

1 . 30ps 

1 .30ns 

2 

1.55 

1.80 

1.55 

1.80 

3 

2.55 

2.55 

2.55 

2.55 

4 

1.80 

1.80* 

1.80 

1.80* 

5 

3.35 

3.35 

3.35 

3.35 

6 

2.85 

2.85 

2.85 

2.85 

7 

4.35 

4.35 

4.35 

4.35 

*If  R6  or 

R7  is  used  wit 

h Mode  4 and  B 

yte  operation. 

add  0.25jis  to 

SRC  time  and  0.25ps  to  DST  time. 

I 1 
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TABLE  B-II 


Basic  Time  (DOP) 


Double 

Operand 

Instruction 

SM0 

DM0 

SMI -7 

DM0 

SM0 

DM1 -7 

SMI -7 

DM1  -7 

MOV,  ADD 

1.75ys 

2.55ys 

4.55ys 

4.55ys 

MOVB 

2.25^) 

4. 05^) 

4.55^) 

5.55(2) 

SUB 

2.80 

2.80 

4.80 

4.80 

CMP,  BIT 

3.05 

3.05 

3.05 

3.05 

CMPB,  BITB 

2.55 

3.55 

3.80 

2.80 

BIC 

2.80 

2.80 

2.80 

2.80 

BIS 

3.05 

3.05 

4.30 

4.30 

BICB,  BISB 

2.80 

3.80 

5.55(2) 

6.55(2) 

XOR 

2.30 

- 

4.05 

- 

iiii 

(1)  Add  0.25ps  if  effective  SRC  byte  operand  is  negative. 

(2)  Add  0.25ys  if  DST  byte  address  is  odd .| 
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TABLE  B-III 
Basic  Time  (SOP) 


CLR,  INC,  DEC,  SWAB  2.05ys  4.05ys 

ADC  2.30.  3.05 

SBC  2.30V5)  3.05 (’) 

COM,  NEG,  ROL,  ASL,  ROR,  ASR,  CLRB,  SXT  2.30  4.30 

TST  2.30, 2.30,  _ 

INCB  2.55 jj)  5.05(3) 

ADCB  2.30{5  (9)  3.80(7|(1(!I) 

DECB  2.55  2 ....  5.053) 

SBCB  2 . 30^ 5 ) 0 ) ) 3.80(8H1OJ 

COMB  2.30  5.30‘;3) 

NEGB  2.80(2)  5.30(3) 

ROLB,  ASLB,  RORB,  ASRB  3.55  5.55(4) 

TSTB  3.05  3.05. 

MFPS  (1067DD)  3.05  4.30(3) 

MTPS  (1064SS)*  3.55  4.05 

JMP  10.44  2.50 

JRS  10.54  5.00 

*For  MTPS,  use  byte  DST  times,  not  SRC  times. 

(1)  Add  0.25ys  if  register  byte  operand  = 377  octal. 

(2)  Add  0.25ys  if  register  byte  operand  = 0. 

(3)  Add  0.25ys  if  DST  byte  address  is  odd. 

(4)  Add  0.25ys  if  DST  byte  address  is  even. 

(5)  Add  0.25ys  if  C-bit  is  set. 

(6)  Add  1 .50ys  if  C-bit  is  set. 

(7)  Add  1.75ys  if  C-bit  is  set. 

(8)  Add  2.00ys  if  C-bit  is  set. 

(9)  Add  0.25 ps  if  C-bit  is  set  and  register  byte  operand  = 377g. 

10)  Add  0.25 us  if  C-bit  is  set  and  DST  byte  address  is  odd. 

11)  Add  0.25ys  if  C-bit  is  set  and  register  byte  operand  = 0. 


(1)  Add  0.25ys 

(2)  Add  0.25ys 

(3)  Add  0. 25ys 

(4)  Add  0.25ys 

(5)  Add  0.25ys 

(6)  Add  1 .50ys 

(7)  Add  1.75ys 

(8)  Add  2.00us 

(9)  Add  0.25ys 

(10)  Add  0.25 ps 

(11)  Add  0.25ys 


TABLE  B-IV 


Miscellaneous  Instruction 


Instruction 

Basic  Time 

SOB  (BRANCH) 

3.05ys 

SOB  (NO  BRANCH) 

2.55 

SET  CC 

2.05 

CLEAR  CC 

2.30 

NOP 

2.30 

RTS 

3.30 

MARK 

4.10 

RTI 

5. 60^) 

RTT 

5.60(2) 

TRAP,  EMT,  IOT,  BPT 

10.15 

WAIT 

2. 75^) 

RESET 

3.25+10usINIT+90ys 

(1)  Add  0.50ys  if  NEW  PS  has  bit  4 (T-bit)  set. 

(2)  Add  0.25ys  if  NEW  PS  has  bit  4 (T-bit)  set. 

- 


TABLE  B-V 


Basic  Time  (Branch) 


Branch 

Instruction 

Branch'^  ) 

No 

Branch 

BR  „ 

2.80 

- 

BNE,  BEQ,  BPL,  BMI , BVC,  BVS,  BCC/BLO,  BCS/BHQS 

3.05 

2.00 

BGE 

3.30(2) 

2.25(2) 

BLT 

3.30(3) 

2.25(3) 

BGT 

3.80(2) 

2.00(4) 

BLE 

3.05(5) 

2.75(3) 

BHI 

3.55 

2.00(6) 

BLOS 

3.05(6) 

2.50 

(1)  Add  0.25ms  if  offset  is  negative. 

(2)  Add  0.25ms  if  N-bit  is  set. 

(3)  Add  0.25ms  if  V-bit  is  clear. 

(4)  If  Z-bit  clear:  add  0.75ms  if  N-bit  is  clear;  add  1.00ms  if 

N-bit  is  set. 

(5)  If  Z-bit  clear:  add  0.75ms  if  V-bit  is  set; 

add  1 .00ms  if 

V-bit  is  clear. 

(6)  Add  0.05ms  if  Z-bit  is  clear. 
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TABLE  B-VI 
Basic  Time  (EIS) 


EIS  Instruction* 

Basic  Time 

MUL 

16.70ps  worst  case 

DIV 

23.40ys  worst  case 

ASH  (RIGHT) 

5.80ms  + 0.25ms  per  shift 

ASH  (LEFT) 

5.55ms  + 0.75ms  per  shift 

ASHC  (RIGHT) 

6.60ms  + 0.25ms  per  shift 

ASHC  (LEFT) 

6.35ms  + 0.75ms  per  shift 

*Use  word  DST  times,  not  SRC  times. 

1 

TABLE  B-VI I 
Basic  Time  (FIS)* 


FIS  Instruction 

Basic  Time 

FADD 

20.15  typical 

22.80  worst  case 

FSUB 

20.40  typical 

23.00  worst  case 

FMUL 

50.25  typical 

69.40  worst  case 

FDIV 

69.40  typical 

70.10  worst  case 

*Inst  Time  (FIS)  = Basic  Time 

+ Shift  Time  for  Binary  Points  + 

Shift  Time  for  normalization. 

Binary  Point  Alignment:  0.25ps  per  shift. 

Normalization:  0.75ps  per  shift. 

^ 1 
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APPENDIX  C 

THEORETICAL  DEVELOPMENT  OF  A BASE-2  FFT  ALGORITHM 


APPENDIX  C 


Theoretical  Development  of  a Base-2  FFT  Algorithm 


Development  of  the  Cooley-Tukey  FFT  Algorithm  (Ref  2:257-265) 
The  DFT  of  N-data  can  be  written  as: 


N~1 


,nk 


X(n)  = l x(k)W  , n = 0,  1,  ....  N-l 
k=0 


(C-l) 


where 


Wnk  = exp  ~ ~1nk , n = 0,  1,  ...  N-l,  k = 0,  1,  ...  N-l 


The  direct  evaluation  of  Eq  (C-l)  requires  the  solution  of  the  following 
set  of  equations. 


X(O)1  = W°x(0)  + W°x(l ) + W°x(2)  . 
X(l)  = W°x(0)  + W1x(l)  + W2x(w)  . 
X (2)  = W°x(0)  + W2x(l)  + W4x(2)  . 
X(3)  = W°x(0)  + W3x(l ) + W6x(2)  . 


. + W°x(N-l ) 

. + WN_1  x(N-l) 
. + W2N'2x(N-1) 
. + W3N~3x(N-1) 


X(N-l)  = W°x(0)  + WN-1x(l ) + W2N"2x(2)  . . . vjCN"1  )N-(PI-1  )X(N-1 ) 


The  above  set  of  equations  can  also  be  written  in  the  matrix  form  given 
below: 


x(0) 

W°  W°  W°  . . 

X(1) 

W°  W1  W2  . . 

X(2) 

w°  w2  w4  . . 

X(3) 

• 

w°  w3  w6  . . 

• 

X(N-1]_ 

w w^w2"'2 

. VT 
. W1 


. w 


N-l 

,2N-2 

3N-3 


.(N-l)N-(N-l) 


X(0) 

X(l) 

X(2) 

X(3) 

X(N-l) 
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Therefore,  in  order  to  compute  a DFT  of  N-data  points,  we  have  to  solve 
the  matrix  given  by  Eq  (C-2),  all  terms  of  this  matrix  being  complex: 


[X(n) ] = [Wnk]Cx(k)]  (C-2) 

2 

This  direct  evaluation  requires  N complex  multiplications  and  N(N-l) 
complex  additions  (Ref* 3:151).  The  FFT  algorithm  developed  by  Cooley 
and  Tukey  (Ref  3:151)  reduces  the  number  of  complex  multiplications 
and  additions,  and  computer  execution  of  this  algorithm  results  in  faster 
computation  time  as  compared  with  the  evaluation  of  Eq  (C-l)  (Ref  3:151). 

In  case  N is  limited  to  a power  of  2 (i.e.,  N = 2P),  it  is  possible 

to  express  k and  n in  terms  of  the  index  p,  as  a binary  weighted  series: 

k = k + kp_22p"2  + ...  k^1  + kQ  (C-3) 

n = ry^-1  + np_22p-2  + ...  n^1  + nQ  (C-4) 


Here  kp_-|,  kp_2,  ...  kQ  and  np_-j,  np_2>  ...  nQ  can  take  a value. of  either 
0 or  1 . In  this  way,  we  can  express  all  the  N possible  values  of  the 
indices  in  terms  of  a binary  number  and,  hence,  facilitate  the  consider- 
ation of  storage  in  the  digital  computer.  Using  this  convention,  Eq  (C-l) 
can  be  rewritten  in  the  following  form: 


x(np-r  np-2* 


"o>  * E 


kQ  0,  k.j=0, 
wn(kp-i2f>  1 


+ k 


P-2 


P-1 
2P-2 


=0  x(kP-r  kP-2’ 


k0)  X 


k0) 


(C-5) 


If  the  first  term  of  the  exponential  for  Eq  (C-3)  is  expanded,  Eq  (C-6) 
is  the  resulting  equation: 


w"Vi2P  1 - wnp-i2P  1kp-i2P  ’ + W p'2~  'P 


1,  2P-1 


,nk 


un  k ^P'1 
w o p-1 


-1..  oP-1  + wnp_22p-2kp_12p-1  + ...  Wnl2  kp-l 

(C-6) 
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This  can  be  simplified  to  give  Eq  (C-7)  because  all  other  terms  go  to 
unity.  This  enables 


unk  _ unnk  ^P’1 

w p-1  = W 0 p-1 


(C-7) 


the  innermost  sume  of  Eq  (C-5)  over  k - to  be  written  as  a shorter 


Fourier  transform: 


Vn0*  kp-2  •••  V ~ L ]=0Xl(kp-r  kp-2’  •••  VW 


nokP-i2f> 


(C-8) 


Unlike  the  complete  Fourier  transform,  this  sum  consists  of  a set  of  N 
numbers  only,  each  calculated  from  two  of  the  original  data  points. 
Subsequent  sums  proceeding  outwards  in  Eq  (C-5)  can  be  calculated  using 
a generalized  recursive  expression  for  the  exponential  term: 


,nWp'q  - „(v.2<h  + vz2"'2  + - "(A-'1 


k_  _2p-q 


(C-9) 


q = 1,  2,  3 ...  p 


The  successive  sums  are  evaluated  according  to  the  equation 


Xq(r»o»  ni  •••  Hq_i , kp_q_-|,  ^p_q_2  ••• 


E xq-l ^n0*  ^1  **"  nq-2*  kp-q  ***  k0^  x 
kp-q=0 

uK-i2"'’  + V22"'2  + :::  "cA-'’2'””' 


(c-10) 


q = 1,  2,  3 ...  p 


To  apply  this  recursive  formulae,  the  initial  set  of  data,  x(k),  is 


first  made  equal  to  xQ(k)  so  that  q = 1 , thus 

x ( k ) - x(kp_i , kp_2  •••  ^q)  = X0^p-1  ’ ^p-2  k0^  ( C— 11) 

This  leads  to  the  derivation  of  succeeding  arrays  in  so  that  the 
final  array  will  be  in  X(n).  Equation  (C-12)  given  below  shows  that 
the  elements  of  X(n)  will  be  calculated  in  incorrect  order,  so  some  sort 
of  shuffling  is  required  to  get  the  result  in  correct  order.  (This 
point  is  clarified  when  the  algorithm  for  a 16-point  DFT  using  this  FFT 
method  is  worked.) 

X(np-1’  np-2  •"  n0}  = xp(n0’  *•*  np-2’  np-l^  (C-12) 

A 16- Point  FFT  Algorithm 

In  case  the  number  of  data  points  to  be  transferred  is  16,  the  value 
of  P = 4.  Therefore,  k,n  can  be  expressed  as  4-bit  binary  numbers  as 
given  in  the  following  equations: 


k = 8k3  + 4kg  + 2k^  + kg 
n = 8n3  + 4n2  + 2n^  + ng 


Using  Eqs  (C-13)  and  (C-14),  Eq  (C-5)  can  be  written  as 

X(n3,  ng,  , ng)  = E E E E x(k3>  kg,  k^ , kg) 

k0=0*  kl=0’  k2=0*  k3=0 

j^n(8k3  + 4kg  + 2 k-j  + kg) 


(C-13) 

(C-14) 


(C-15) 


Using  Eqs  (C-7)  and  (C-8),  the  innermost  summation  can  be  written  as 

8k.nn 

x-j  (nQ,  kg,  kj,  kQ)  = E xn(k,,  k,,  k, , kn)W  (C-16) 


3=0 


'O'  3*  2*  1*  O' 
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Subsequent  summations  can  be  written  using  Eqs  (C-9)  and  (C-10): 


x2^0*  ^1*  kl  ’ k0^  ~ ^ VV  k2*  kl  ’ k0)W 

k2=0 


4k2(2ni+  ng) 


(C-17) 


x3(ng,  n1 


k0) 


x9(n 


2'  O’  "1 


k,,k0)W 


2ki  (4n?  + 2n-|  + nn) 


1 


(C-18) 


1=0 


x4^n0’  nl * 


n3} 


x3^n0’ 


V «'2* 


k0)W 


k0(Sn3  + 


4n2  + 


2n, 


nQ) 


(C-19) 


0=0 


X(n3>  n2»  i nQ)  — x^(0g,  n^,  ^2*  (C-20) 

The  set  of  equations  that  result  when  summations  of  Eqs  (C-16),  (C-17), 
(C-18),  (C-19),  and  (C-20)  are  expanded  is  given  in  tables  C-I,  C-II, 

C- III,  C-IV,  and  C-V,  respectively. 

Table  C-V  shows  that  the  final  array  for  X(n)  is  in  bit- reversed 
order;  therefore,  to  get  the  result  in  order,  a post-shuffling  needs  to 
be  done.  Another  way  to  overcome  this  problem  is  to  bit  invert  the  input 
array  [xg(k)],  and  the  output  will  automatically  be  in  order  (Ref  3:178). 
The  signal  flow  graph  for  the  16-point  FFT  without  bit  inverting  the 
data  is  shown  in  Figure  C-l,  and  Figure  C-2  shows  the  signal  flow  graph 
for  data  which  has  been  bit  inverted  before  being  transformed. 
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TABLE  C-I 

Expansion  Results  of  Eq  (C-16) 


Xj(0,  0,  0,  0) 

Xj(0,  0,  0,  1) 

x-j  (0,  0,  1,  0) 

x-|(0,  0,  1,  1) 

x-j (0,  1,  0,  0) 
XjCO,  1,  0,  1) 
Xj(0,  1,  1,  0) 
x1  (0,  1,  1,  1) 
x-|(l,  0,  0,  0) 

Xj(l,  0,  0,  1) 

x^l,  0,  1,  0) 
x^l,  0,  1,  1) 

Xjd,  1,  0,  0) 
Xj(l,  1,  0,  1) 

x1(l,  1,  1,  0) 
x^l,  1,  1,  1) 


xQ(0,  0, 
Xg(0,  0, 
Xg(0,  0, 

Xg(0,  0 , 

Xg ( 0 , 1 , 
Xg (0  » 1 » 

Xg(0,  1 » 
Xg(0,  1 , 
Xg(0,  0, 
Xg  ( 0 , 0 , 
Xg  ( 0 , 0 , 
Xg(0,  0 , 
Xg (0 , 1 > 
Xg(0,  1, 
Xg(0,  1 » 
Xg(0,  1 » 


0,  0)lT 

0,  1)WC 

1,  0)WC 
1,  1)WC 
0,  0)WC 

0,  1)WC 

1,  0)w° 
1,  1)WC 
0,  0)WC 

0,  1)W° 

1,  0)wc 

1,  1)WC 
0,  0)WC 
0,  1)WC 


+ x0^  ’ 

+ x0^  ’ 

+ x0^  ’ 

+ x0  ^ ’ 

+ Xg ( 1 » 

+ Xg(l  , 

+ Xg(l, 

+ Xg(l, 

+ Xg(l , 

+ Xg(l, 

+ Xg  ( 1 9 

+ Xg (If 
+ Xg(l , 
+ Xg(l, 


i,  o)wu  + xnd. 


1,  1)WL 


+ Xg(l, 


o,  o,  o)wu 

0,  0,  1)W° 
0,  1,  0)W° 

0,  1,  1)W° 

1,  0,  0)W° 
1,  0,  1)W° 
1,  1,  0)W° 
1,  1,  1)W° 
0,  0,  0)W8 
0,  0,  1JW8 
0V  1,  0)W8 

0,  i,  nw8 

1,  0,  0)W8 

i,  o,  DW8 

1,  1,  0)W8 
1,  1,  1)W8 


TABLE  C-II 


Expansion  Results  of  Eq  (C-17) 

x2 (0,  0,  0,  0) 
x2(0,  0,  0,  1) 
x2(0,  0,  1,  0) 
x2(0,  0,  1,  1) 
x2(0,  1,  0,  0) 
x2(0 , 1 » 0,  1 ) 
x2(0,  1,1,0) 
x2(0,  1,  1,  1) 
x2(l,  0,  0,  0) 
x2(l,  0,  0,  1) 
x2(l,  0,  1,  0) 
x2Wl,  0,  1,  1) 
x2(l,  1,  0,  0) 
x2(l,  1,  0,  1) 
x2(l,  1,  1,  0) 
x2(l,  1,  1,  1) 
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= Xj(0,  0,  0,  0)W°  + 
= x-j (0,  0,  0,  1 )W°  + 

= Xl (0,  0,  1,  o)w°  + 

= Xl(0,  0,1,1 )W°  + 
= x^O,  0,  0,  0)W°  + 
= x] (0,  0,  0,  1)W°  + 
= x^O,  0,  1,  0)W°  + 
= x-j  (0,  0,  1,  1 )W°  + 
= Xj(l,  0,  0,  0)W°  + 
= Xl(l,  0,  0,  1 )W°  + 
= xT(l,  0,  1,  0)W°  + 
= Xjd,  0,  1,  1)W°  + 
= (1 , 0,  0,  0)W°  + 

- Xld.  0,  0,  1)W°  + 
= Xjd,  0,  1,  0)W°  + 

= Xjd,  0,  1,  1 )W°  + 


Xj (0,  1,  0,  0 )W° 

x-j  (0,  1,  0,  1)W° 
x-j (0,  1,  1,  0)W° 
x-j  (0,  1,  1,  1)W° 
x-j  (0,  1,  0,  OJW8 
x-j  (0,  1,  0,  UW8 
x-j  (0,  1,  1,  0)W8 
Xjfo,  i,  l,  DW8 
x1  (1 , 1,  0,  0)W4 
x^l,  1,  0,  1)W4 

X]d,  i,  i,  o)w4 

Xl(ls  1,  1,  1)W4 
Xl(l,  1,  0,  0)W12 

Xld,  i,  o,  i)w12 

X-j  ( 1 , 1,  1,  0)w12 
x^l,  1,  1,  1)W12 


TABLE  C-III 

Expansion  Result  of  Eq  (C-18) 


x3(0,  0,  0,  0) 
x3(0,  0,  0,  1) 
x3(0,  0,  1,  0) 
x3(0,  0,  1,  1) 
x3(0,  1,  0,  0) 
x3(0,  1,  0,  1) 
x3(0,  1,1,0) 
x3(0,  1,  1,  1) 

x3(l,  0,  0,  0) 
x3(l,  0,  0,  1) 
x3(l,  0,  1,  0) 
x3(l,  0,  1,  1) 
x3(l,  1,  0,  0) 
x3(l,  1,  0,  1) 
x3(l,  1,  1,  0) 
x3(l,  1,  1,  1) 


= x2(0,  0,  0,  0)W°  + x2(0,  0,  1,  Ojli3 

= x2(0,  0,  0,  1)W°  + x2(0,  0,  1,  1 )uP 

= x2(0,  0,  0,  0)W°  + x^O,  0,  1,  0)/ 

= x2(0,  0,  0,  1)W°  + x2(0,  0,  1,  1)1^ 

= x2(0,  1,  0,  0)W°  + x2(0,  1,1,  0)W4 
= x2(0,  1,  0,  1)W°  + x2(0,  1,  1,  1)W4 

= x2(0,  1,  0,  0)W°  + x2(0,  1,  1,  0)W12 

= x2(0,  1,  0,  1)W°  + x2(0,  1,  1,  1)W12 

= x2(l,  0,  0,  0)W°  + x2(l,  0,  1,  OJW2 

= x2(l,  0,  0,  1)W°  + x2(l,  0,  1,  DW2 

= x2(l,  0,  0,  0)W°  + x2(l,  0,  1,  0)W1C 

= x2(l,  0,  0,  1)W°  + x2(l,  0,  1,  1)W1C 

= x2(l,  1,  0,  0)W°  + x2(l,  1,  1,  0)V^ 

= x2(l,  1,  0,  1 )W°  + x2(l,  1,  1,  1)^ 

= x2(l,  1,  0,  0)W°  + X2(l,  1,  1,  0)W14 

= x2(l,  1,  0,  1)W°  + x2(l,  1,1,  1)W™ 


TABLE  C-IV 


Expansion  Result  of  Eq  (C-19) 


x4(0,  0,  0,  0) 
x4(0,  0,  0,  1) 
x4(0,  0,  1,  0) 
x4(0,  0,  1,  1) 
x4(0,  1,  0,  0) 
x4(0,  1,  0,  1) 
x4(0,  1,  1,  0) 
x4(0,  1,  1,  1) 
x4(l,  0,  0,  0) 
x4(l,  0,  0,  1) 
x4(l,  0,  1,  0) 
x4(l,  0,  1,  1) 
x4(l,  1,  0,  0) 
x4(l,  1,  0,  1) 
x4(l,  1,  1,  0) 
x4(l,  1,  1,  1) 


x3(0,  0,  0, 
x3(0,  0,  0, 
x3(0,  0,  1, 
x3(0,  0,  1, 
x3(0,  1,  0, 
x3(0,  1,  0, 
x3(0»  1 > 1 * 
x3(0,  1,  1, 
x3( 1 » 0,  0, 
x3(l,  0,  0, 
x3(T,  0,  1, 
X3O » 0,  1, 
x3(l»  1»  0, 
x3( 1 » 1 » 0 1 

X3 ( 1 » 1 » 1 t 

x3( 1 t 1 » 1 » 


0)WU  + x3(0,  0,  0,  1)W 
0)W°  + x3(0,  0,  0,  1)WS 

n\i  |0  1 t r\  r\  1 i \ii^ 


0,  0,  1)WC 
0)WU  + x3(0,  0,  1,  1)W4 
0)W°  + x3(0,  0,  1,  l)w] 
0)W°  + x3(0,  1,  0,  1)W2 


0)WU  + x3(0,  1,  0,  l)wa 
0)W°  + x3(0,  1,  1,  1)W6 
0)W°  + x,(0,  1,  1,  l)w] 
0)W°  + x3(l,  0,  0,  1JW1 


0)WU  + x3(l, 
0)w°  + x3(l; 
0)W°  + x3( 1 . 
0)W°  + x3(l, 
0)W°  + x3(l: 
0)W°  + x3(l; 


0,  0,  1)WS 
0,  1,  1)W5 

0.  1,  nw1 

1,  0,  1)W3 
1,  1,  1)W] 
i,  i,  nw7 


0)wu  + x3(l,  1,  1,  l)w' 
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TABLE  C-V 


Expansion  Result  of  Eq  (C-20) 


x(0. 

0, 

0, 

0) 

= 

* 

o 

X 

0, 

0, 

0) 

x(0. 

0, 

0, 

1) 

= 

x4(l. 

0, 

0, 

0) 

x(0. 

0, 

1, 

0) 

= 

0k 

0 

X 

1, 

0, 

0) 

x(0. 

0, 

1, 

1) 

= 

X40 » 

1, 

0, 

0) 

x(0. 

1, 

0, 

0) 

= 

*4(0 » 

0, 

1, 

0) 

x(0. 

1, 

o,  1) 

= 

x4(l. 

0, 

1, 

0) 

x(0. 

1, 

1, 

0) 

= 

x4(0. 

1, 

1, 

0) 

x(0. 

1, 

1, 

1) 

= 

x4(l. 

1, 

1, 

0) 

x(l. 

0, 

0, 

0) 

= 

x4(0. 

0, 

0, 

1) 

x(l  , 

0, 

0, 

1) 

= 

X4  ( 1 » 

0, 

0, 

1) 

x(l. 

0, 

1, 

0) 

= 

X4  (0 , 

1, 

0, 

1) 

x(l  , 

0, 

1, 

1) 

= 

x4(l. 

1, 

0, 

1) 

x(l. 

1, 

0, 

0) 

= 

O 

X 

0, 

1, 

1) 

x(l. 

1, 

0, 

1) 

= 

x4(l. 

0, 

1, 

1) 

x(l. 

1, 

1, 

0) 

= 

x4(0. 

1, 

1, 

1) 

x(l. 

1, 

1, 

1) 

s 

x4(l. 

1, 

1, 

1) 
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APPENDIX  D 

CALCULATION  OF  DFT  OF  REAL  DATA 
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APPENDIX  D 


Calculation  of  DFT  of  Real  Data 


Consider  a real  2N  point  sequence  x(k);  the  DFT  of  this  sequence  can 
be  carried  out  using  Eq  ( D- 1 ) : 

2N-1  nk 

X(n)  = E x(k)W™  n = 0,  1,  ...  2N-1  (D-l) 

k=0 

nk  9 .. 

W = exp  - , ...  n = 0,  1 ...  2N-1  (D-2) 

2N  . 


Equation  (D-l)  can  be  rewritten  in  the  following  form: 

N-l  2rn  N-l  (2r  + l)n 

X(n)  = E x(2r)W  + E x(2r  + 1)W  (D-3) 

r=0  2N  r=0  2N 

N-.l  Dn 

X(n)  = z x(2r)  WKn  + w" 
r=0  n H 

Equation  (D-4)  has  resulted  by  splitting  x(k)  in  (D-l)  into  even  and 
odd  sequences  as  shown  in  Equations  (D-5)  and  D-6)  respectively. 


y(r)  = x(2r)  even  sequence  (D-5) 

z(r)  = x(2r+l)  odd  sequence  (D-6) 


A complex  sequence  v(  r)  is  now  defined  such  that  the  even  points  of  x(k) 
[Equation  (D-5)]  becomes  the  real  part  of  v(r)  and  the  odd  point  of  x(k) 
[Equation  (D-6)]  becomes  the  imaginary  part  of  v(r). 


v(r)  = y(r)  + jz(r)  (D-7) 


N-l 

E 

r-0 


x(2r+l)  W, 


rn 


N 


(D-4) 


154 

A 


An  N-point  DFT  of  V(r)  yields  V(n),  which  is  given  in  Eq  (D-8): 


N-l 

V(n)  = E v(r)W™ 
r=0  n 

(D-8) 

or 

N-l  ™ N-l 

V(n)  = E y(r)W™  + jE  z(r)wjn 
r=0  n r=0  n 

(D-9) 

or 

V(n)  = Y(n)  +jZ(n) 

(D— 10) 

The  DFT  of  real  and  imaginary  parts  can  be  written 

in  terms  of  V(n)  as 

(Ref  5:319). 

Y(n)  = l/2[V(n)  + V*(N-n)] 

(D-ll) 

and 

Z(n)  = l/2j[V(n)  - V*(N-n)] 

(D-12) 

here 

V*(N-n)  = Y(N-n)  - jZ(N-n) 

(D-13) 

The  result  obtained  in  Eqs  (D-ll)  and  (D-12)  can  be  used  to  rewrite 
Eq  (D-4)  as: 

X(n)  = Y(n)  + wJNZ(n)  (D-14) 

or 

X(n)  = l/2[V(n)  + V*(N-n)]  + l/2j[V(n)  - V*(N-n)]  (D-15) 

This  equation  is  true  for  O^n^N.  [Note  that  V(N)  = V(0)  ]. 

(Ref  5:320-321) 
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APPENDIX  E 

LSI-11  ASSEMBLY  LANGUAGE  CODE  FOR 
FFT  SUBROUTINE  (FOURIE)  AND 
EXECUTION  TIME  COMPUTATION 


APPENDIX  E 


LSI-11  Assembl;,  Language  Code  for  FFT  Subroutine 
(Fourie)  and  Execution  Time  Computation 

(Ref  6:288-293) 


<r  *-  o <x 

Q <1  r- 

n mu 


u x:  « 

-J  3 

U.  C i/>  i/> 

O t—  O') 

o <z  z 

• o *-« 

Ui  o 

OiZ  X CL 

X C- 

t-  X 

</>  C>  UJ 

(Ac  *-*  _J 

*-•  Z CL 

-*  CVJ  E 

Z O Z O 

— Cu  — O 


►-  x a>  >-  •->  z ac 

Z r-  * r-  — O 

UJ  * OC  O'  O U. 

z o)  rs,  iuaa.(A 

O UJ  it  t—  »—  Z 

u.  CM  Z CO  -J  <X 

X u.  Z C-.  cu 

UJ  CO  O lb  Ui  c 

-*  a.  X QL 

U!  t-  CL  C-  UJ 

x <r  to  p-  go 

»-  0)  x *-•  z 3 

CO  c*  u.  O 

<o  x z a « z p- 

#-«  o o ►—  c-  c* 

*-•  u.  « z UJ 

00  X ►*-  OOu.CZ 

— 3 I - O O <X 


UJ  I X ■ 
£ I K» 

<X  I UJ 
QC  I UJ  Z 
C I (KH 

a.  i 4c 


* UJ  UJ  o 

: x •-«  o • 
: »-  u 

3-OW 

CZ  O -*  UI 
JUJU.  3 

: l-  cz  _j 

J U.  <CUJ  c 

: <r  t > 


© I © CO 
*-*  I mJ 
»-  I J »- 
H • ©IL 


u.  I UJ  UJ 
UJ  | X X 
Q I J-J- 
I 

• I X 
N I 


ono<t 

o -•  o 

OhON 

o o o o 
o a o o 
o a oo 
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BE  ACQUIRED  AND  STORED  IN  THE  2ND  BUFFER  WHILE  PRBC 
-EBSINC  IS  BEINC  DONE  ON  THE  1ST  BUFFER. 


I o oc 

I I—  <X  • IU 

: m u.  u>  *- 

: ai  »-  ►-  2 

► • Z O 

. 2 (A  0.  *-»  O 

► m CO  UJ  O O 

’ d Oi  UJ  U. 

: r-  ct  ^ o 

: 2 O <X  2 

■ O Q O »—  •-* 

U<If  «t- 

I O U. 

i _i  a ^ m 

J iJ  UI  u.  X 

► *-«>-  r-  © CO 

■ 3 OC  2 

aj  3 » i- 

: M > c r-« 

■ i2  2 O U.  CO 

~ M — o *• 


m . X 
J 0>  CO  CM 
lOUJ-Ct 

> a.  h-  r •- 
i O O O 2 
ioo;  i h 


<r  o > 
ui  o uj 

H-  m Ui  0.  OC 

2 a.  0)  «X  O. 

mOCi;3 
O U UltA  O (L 

a.  > r-  <x 

ui  UJ  CC  3 
W j it  O uj  W 

»-*  <X  a.  IX 

X H-  4 UJ 

H-  m m >*  u.  _J 

u£  Qj  U.  E GJ 

I-  O O => 

<X  a.  2 O O O 

— O M — — O 


Z I HH 

O I © u. 

h i a. 

(A  I UI 
OC  I « X 
UJ  I HH 

> i <r 

2 I O H-  • 

Ml  <X  OL 
I H I UJ 

»-  I 3 H-  O 

m i a.  oc 

co  I 2 o o 

I m <A 
Ui  I H- 

X I UJ  — ■ O 

H-  I X O UJ 

■ H UJ  oc 

• I -J  oc 


A 

oc  «e 
- oc 

CM  - N 
Z CO  OC 


> > OC  t) 
O O -I  -i 

z zoo 


> a.  h- 

© z _J 

z o co 


J JJOOOCO.K 
iA0)(AQQ(AZ<l 

ca«x<c<->ua) 


««  M O 

hhA  A 
A A N N 

N*  P»  M M 

n r«*  o o 

M M o o 


f»  O m M -»N*< 

•O  O O O O O *r 

O O O r*  NON 

a a m -o  -e  m o 

O O O m — o o 

o o o o a a o 

^ O CM  o *r  *o 

m cm  cm  cm  nnn 

CM  CM  CM  CM  CM  CM  CM 

a o o o o o o 

o o o o o o o 

o o o o o o o 

iNH  <T  NCO-O  *<N' 


mm  trt 
O O M 

o cm  *r 

O O CM 

M cm  O 


MNNNMNNN 

OOOOO-OCM* 

nnrtNNNHN 

COCm3*COCMO 

□ooooaao 


<ON^  O-COI 
A 4 4 4 N N O ( 
N N N N N N DC 

oooooooc 

oooaoooc 

oooooooc 


Nn>rA4Ntto-OMNn4in4N«9‘OMNn<rA4N«^ 

NNNNNNNNnr»r)r)r)nnrtnn'C>r^«ff4^^^< 


FOURIE  — FFT  SUBROUTINE  RT-U  MACRO  VM02-12  30-0CT-76  23142104  PACE 


CO  I z o 
- I C I Hlk 

in  i u.  i a 

O > 00  >—  I I 3UI 

Z O O O I • I Q X 

*-•  X CO  Z I N I H* 


- I 
Ul  I 

> > cjc  > > q.  > >oco>>oflD>Q.a.ia  t-  i 

o O (ft  0 0X0  00300330XX0  O I 

X X X X C X X C 0/  X C <1  O)  X u o V>  Z I 


0*4  N 

o o n 

NO  n 

to  o m 

o *4  m 

o o o 


O -4  m O *-* 

o o-o  o o 

r»  mm  mo 

-C  <C  <2-  -o  o 

-I  o -< 

o o o o o 


-4  in  nonhn 
MO  O M 3 M O 

— m o — cm  o 

CM  (M  O — 

CM  — -0  >C  — 

O O O O O O 


o M -4  -4  *4  rt  M 
CM  3 M M M -4  -4 

— cm  o o iO  *r 
-4  o m m r*  o 
-O  «C  -•  M M M O 
O *4  O O O O O 


-0  O 

O *4 

n n 
o o 
o o 
o o 


O <r  O *r  «0  O vf'OON.f-OON^-OON 

-4  mm  ri  n n >r  ^^iiTuTiOuT-o-o-O-onn 

n r*.  n n r,  n nnnmnr'imn.-onnrt 

O 3 3 3000  000.3  00000000 

O OO  0 3 0 0 OOOCOOOOOOOO 

O CO  oaoo  030000000000 


HNn^n4N40'04Nn<rh,)4M>^0>«NI<)-rin>0\C>04Nn4^4N«(h04Nl,)>rin4N«^04N 

v4a4V4«^«4V4«4«4«4«4MMMMMMMMM(xr*)r)r)nr)ror<)r*)nn'r<^'r'^*r<r<c>e'r>ronA 


160 


o • 


FOURIE  — FFT  SUBROUTINE  RT-U  MACRO  VM02-12  30-OCT-7B  23H2104  PACE 


(A  Z V)  Z 
O »“*  O ►>» 
uwu  w 


^ I 

* c* 

m s. 

u.  z 


a > 
ui  o 
x u 
3 ^ 


a.  o> 
<x 
x a. 

«r  o 
z 

O CM 


>-  z a _j 


« o-)  iu  uj  <z  a 

<c  a u a.  - 

o uj  - - - o 

tt  J*  -F  (M 

n ai  • » m m 

X 


> > > > 
o o o o 

X X X X 


o »-• 

• o 

x a.  jk 


x - 

- u 

Z ft. 


> X > > 

O J)  O O X 

X ->  XXX 


o 

X 


o 

z 

X 


4N4N 
■f  -r  rtrt 
NNNN 
NNNN 
N N N N 


FFNO 
0 0 0-4' 

o o o o 

0 3 0 0 
0 3 0 0 
O O O O 


-4  O 
f*  N 
r«»  ^ 
n 


O CM 

m r* 
n r* 
r-*  o 
r*  a 

•4  O 


NNNN 
«0  -O  -O  -C 
r**  r-  rs. 
N CM  CM  M 


-4  r* 
O -O 

r*  r*. 
<o  >r 
-4  o 
o o 


4NDFII)4N«>0 


^NDFIflONO  QHNn 


04<FN 
NNPJ-F 
-T  «r  *r  ■# 

o o o o 

3 3 0 0 
O O O O 


«e  o 

m <o 

<r 

o o 
o o 
o o 


•o  o 

FFKJ 

o o o 
o o o 
o o o 


X 

§ 


0HNnFu14N«»04Nr)Fn<IN<»04NI1<f|||4 

NNNNNNNNNNnnnnnr)nnnnFFF-F<F«e-r 


162 


CONTINUED  ON  NEXT  PACE 


in  ui  u|  • 
OC  CC  CC  UJ 
UI  Ui  UJ  X 
X X Z Ui 


CD  CD  CD  <0 
iO(A(A(A 


>->■>>■ 

Z CD  CD  CD 
*-.0  0 0 
cd  o o o 


o 

Ui 

X 


UJ  B 

z iD 

o 

W U N 


z 

>- 

OC 


I—  u.  ul 

o a 

CD  CD  cl  • 

•—  UI  -4  — 

X o x _j  > 

f-  _J  r-  <1 


o • 

<x  z 

U G>  • • 

Ll  *-•  — CM 
>-  >-  * 
oc  a ~ 

O CJ  CD  I— 
U.  *-•  O Ct 


u.  cm  — 
*-•  ►-  o:  _j  >• 


LL  o'  W>  3 >-  t 


> Z CD  . 

> UJ  CD 
4 o')  Ui 
1 UI  Cl 

Ct  Q 
• Cl.  a 
«x 


r—  O h-  * 
UJ  ►—  UJ  CM 

O CD  O CD 


U _J  O I—  * 


CD 


> — 
>- 
CD  — 
*-•  Z 

CD  CD 
■ *-*  * 

I X N 


CD  OO 

<x 

CD  r-» 


<X  UJ 
CJ 

z <x 

*-<  OC 
X o 
I—  I— 
CD 


— <t 
ui  ^ 
OO- 
0.  3 
3X0 
UJ  UJ  -J 
Z *-  UJ 

0 

X UJ 
O X O 
C O UI 
Ui  CD  I 


X 

o 

OC 


Z Z CD  CD 

r-  -4  O O 

0/  CD  CJ  CJ 
NNNN 

t-4  OC  »-4  Z 


OC 

CD 

: x X 


f--.r-.fSIZ  ^ 


> cc  > > o > oc 
o _*  o o o o cd 
x o x x <r  x -> 


o 
— « 

o 


CM  CM 
CC  v 

UJ  - O CM 

Z 3 .J  CD  X CM 

Z O - X CM 

<D  O — - X 

- a - * -4  cd  - 

Cm  CJ  CM  OC  O -4 

CC  m ll.  it  — U • 


> 2 CC  > > _J  I 
O O CD  O O 3 CD 
X <=  3 X X X «X 


CD  Z 

O CM  CM  *-«  CM  CM 

o cc  cc  cd  oc.  x 

CM  - *■  CM  CM  »r  ■*■  * - CM 

X ~ UJ  X X X ■ Ui  oc 

- Z * * - -4-M-Z  - 

cc'~i/>mcccc~~<f>m 


>>_|X>>»->_JX 
O O 3 CD  O O CD  O 3 CD 

cccczrci-rcc 


>0 

o 

o 

XX 

OM 


Z X 
*-•  <x 

>-  Cl 


I I 
I I 

Ui  UJ 
H W 

X X 

S3 


cm  cm  «r 

•O  N Ui 

4 <0  4 
r* 

N M N 


*f  O O O *D 
*r  m *0  *r  o *r  n 

-0-0  0 0 0 0-0 

M N « fM  O -4  *M 

MNO  fM  O fM 


o -4 

CM  O 

o o 

fM  O 
PM  O 


-*-4  0 >4  0 0-4  -4  0 


o -4  cm 
o 000 

fM  OOO 

fM  fM  O P- 

rM  jm  a fM 

-4  -4  0-4 


o -4 

in  o 

ui  o 

fM  O 
fM  O 

-4  O 


OOOO 

OOOO 

OOOO 

OOOO 

OOOO 

OOOO 


V o o 

0 0 0-4 

it1)  *1  in  in 
000c 
— : o o 
> o o 


o . 

O I 


fM  rM  O -4  — . CM  fM 

0 0 0 0 3 0 0 

fM  O fM  O S fM  fM 

o in  O O O o -r 

-*0-4-4-0-4  3 

0000000 


CM  O -r  O CM  O CM 
-.[MNnnn-f 
in  ui  in  in  ui  in  in 

0 0 0 0 3 0 0 
0 0 0 0 0 0 3 

0000000 


rM  cm  pm  rM  cm  fM  rv 
-fl  O -O  -C  O < N 
CM  fM  fM  Psj  4 CM  f*J 

o cm  «r  o -4  o n 

4<C4«>N 

OOOOOOO 

-d  cm  -o  cm  *c  o «r 
-r  D D -O  ■<.  fM  N 

ui  ui  ui  ui  in  in  ui 
coooooo 
0000000 

OOOOOOO 


rMCMfMNrM-r^rCMfMfM 
-CO-CCM-COCMO-CCM 
N-NfM^  — N-CNN 
o -4  a n o o -i  -4  o n 

.0— «f^rM 

OOOOOOOOOO 

D-r-ON-ON-r-OO-/ 

occ-4-4CMCMrsPifi 

-o-c-c-o-o-c-o-o-o-o 

OOOOOOOOOO 

0000300000 

OOOOOOOOOO 


-4wr>-#tn-opM«o-a-4 


O-4CMn-rui-DfMc0OO-4fMn-ru1*0r^O£>O  — fMOvrcl-orM<io-0 
H44H-4-4-4NNNNNNNNNNnnnnnnnnnn-f 


163 


Cl  00G6C0  010267  177636  MOV  R2.I2S1NY  I J 2S1NY- I 2«SIN ( Y > 

C?  0006CC  011C02  . MOV  (RO.R2 

C3  0006C6  070267  177SC2  MUL  C0S.R2  II2»C0S(Y>. 

CC  000632  073227  000001  ASHC  11. R2 

CS  000636  010267  177626  MOV  R2.I2C08Y  I J2C06Y-12«C0«<YI 


164 


f .UK IE  — FFT  SUBROUTINE  RT-II  MACRO  VM02-12  30-0CT-7B  23UZIOA  PACE  10 

POST  PROCESSING  ALGORITHM 


Z CL  Oft 
O C •« 
Q O C ft. 


O 

z 

<x 

Q 

z 

<X  £ 


Oft 

o 

U 


CA 

t-  »-  Z O 
O 

Z O)  W • 

p-«  o * x 
oiuto; 


■ — u z 


Q ♦ 
z — 
<z  r 


z u.  o> 

P-*  3 O Z - 

oft  a.  03 

<6  Z I-  Pi  O 
M r-1  3 r- 
C->  O.  C OJ 
O M P-  3 02 
X X X 0 

ihOUJZ 

Ul 

P-  UJ  lil  Ui  > 
OUT  I I « 
O 3 h-  »-  O 


~ ~ ~ Z 
»—  P—  *— « 

o to 

</)  z o * 

Om  * ft. 

UO)u.H 
» * r-  I <0  ft£ 

ft.  Z I C ->  n 

u.  *-•  u a. 

♦ l U I Z ftC 

■ft-  X U It 

a;  * • 

II  II  X X z — — 

— — i i n r r 

x x z z XI  l 

w * z z 


v>  - 

UJ  3 
-J  o 


II  X u X Ul 
X Z X z o 

CC<I«H 


— — Oft 


X c 

LU  — ' 

U.  »-•  QtC 


x o.  x 

3 OC  U. 


at  u.  ooooooooo 


: « - ..  ~ ..  X Z X z X 

ft.u.Z£  jCiHHX 

I 0> 

: <r 

> a. 


ooooooooo 

ooooooooo 

ooooooooo 

ooooooooo 

ooooooooo 

ooooooooo 


-r-oaM'T'OOt**# 

■^^cs.tNirMtsip'inrt 

ooooooooo 


ooooooooo 

ooooooooo 


-<CMn^i0-flN«»0^Nn-rin’0Nft)*‘O-<Nn'r«T-0N«)fra-*Nn^irt<N«(>'O«MNn^ 

^^WH»««^«^«NNNNNCMNNNNnnnr»r)r»rtnnn«p<^## 


166 


CONTINUED  ON  NEXT  PACE 


FOURIE  — FFT  SUBROUTINE  RT-11  HACRO  VH02-12  3O-0CT-7B  23H.2>0<.  PACE  11 

POST  PROCESSING  ALGORITHM 


N <0  N 4 NO 

«r  in  1/1  ■«  n 
O O O O O O 


in4i04N«>OHNn4gi 


167 


CONTINUED  ON  NEXT  PACE 


168 


169 


• r-  o e-  oc  < 

I KHf  Ui 
lUJSuiZU 
l b)  h (A  >•  < 

O _J  • 

■r  ui  n Z it  H 

-•  _i  <t  Z 

> I <Z  I H H 

, f-UhHIUO 
: mMhuiia^ 
03  * CO  O * 

03  - Oi  H- 

KAUMAOtUX 


- • a >- 

) * »—  u.  *-* 

• uiUOd 

: -j  <x 

i <t  u.  m ui 
low  z 

> V>  • ui  o 

- X 

lUUJr> 

: u,  j at 

03  Li. 

<r  o »- 

* v >-*  x 
« W>  a *£  O 
: <r  o 

* >ca 

: ac 

'OLUI-O 
Z X Ui 

> <X  ►-  03  ►- 

< z a.  u. 


- O M O 

♦ t>i  o ui  o n «r 

■**  -j  a o o j oc 

a <z  o 3 o «z  - 

C*  O (V.  < Q N O h 

'■'(/>  u m u m 03  m 


o >-  ui  ►-  c*  » 
UJHZHIUO 
Z 00  (fi  CO  CO  x 


tvj  NN^NN-e 

p o o o -*  o a o 

oo  *rrvOrs.^f*h 

fNiO  u~>  c*  c*  •+  p* 

-«  o onono** 

O OOOOOO 


■O  O M 

m •«  -o 

*c  «r 


«r  •<»  0N«0^4 

*o  -o  rvr^r>ooo 

>r  >r  <r  <r  <#  in  M)  if) 


• tNin^iO'Or^«(hO«^<^r)«#in-or^*oo*o^c»jr)'Tin«op^«o(ho^fMn«Tin«orv«thO-N*^n^'ift^r**« 

M^*iHHMMM>4MNNNNNNNNNNnr)r)r)nnnnr)r>>r<«<r<<<<f<r<r« 


171 


c*  ‘ >— 


. *03 
>NUI«I 
QC  QC 
> 3 O 

• Q Wt- 

: z z 
; <z  ui  uj 


(AXtf 
►-•A3 
o o 


<c 

O >-  UJ  o 

hfflzo  • 

H JN 
O O f-  QC 
Ui  ►-  3 >- 

A 0X0 
3QXOZ 
UJ  J3S« 

WK3UJ 
M Z (O  t 
*-«  QC 
UJ  O CO  u. 

ZQ.HOZ 
X •-» 


• </> » 


A 


CO  Q 
O O CO  UJ  UJ 
X-ZXQ 
CU  I-  Q <x 
X <Z  J o o 

v V U J d J 

o <c 

s/i  -J  U Ui  Ui 
~ X QC 

x >-  uj  ►-  « 

►-  QC  QC 

O O *—  Q 
CU.CUJ 

uj  uj  x a. 

x cc.  •—  » 


rt  QC  — 
QC—  0* 

- - QC 

^ ISJ  » 

QC  QC  rt  O 

— — oca. 


> > > (• 
o o o »- 
cccac 


<r 

3 

CO 


3 

O 

QC 


a 


I 

C*4 

g 


O 

QC 

o 


4. 

3 


QC* 

Aft. 

3« 

AS 


O 
QC  QC 

©3 


a. 

a 

3 

A 


n^NN 

o o 

HNrt  N 

— • — • a o 

H H H o 

o o o o 


<OON# 
m 4 4 4 
QQQQ 


JNn^w4N4o-o-«Nrt'rin4N«^o^Nn^in4r'«>o. 


172 


FOURIE  --FFT  SUBROUTINE  RT-U  MACRO  VH02-12  30-0CT-7B  23H2I06  PACE  17 


I -O 

J _1  O -I 

Z <X  ~ u UJ 

i > <0  I (Q 

» s 
ZCU1Z 

J UJ  * 0‘>  UJ 

C UJ  <-!«  > 

- »-  CL  O *-» 

X ^ (J 

5 *-*  2 C/> 

-«  CO  *-«•-«  </) 
05  I F*» 

<0  p- 

■*  <Z  uC  Ui 
£ X O 2 _l 

> U.HQ1 

J « . <x 
:nuj  • f- 

0 CD  Ct 

<£C«U1 
H >>  2 

iZUZH 

r o>  uj  * (/> 

C * 3 O 
JClU 

- uj  <r  — i 

- _i  > co  UJ 

D Qj  0 2 

> <r  h-  o *-• 
f-  X 05 

5 o a 

1 • Pi  o Ui 

z ui  ui  u.  x 


m 

m o o o 

m ~ -«  >r 

o r>  i*5  m 


O N M M 

ox  m at 

o o fSl  m . 

o -j  x m o 
•r  i/i  o - - 

• U _l  • » 


»-  UJ  05  O 
hzq:hq 

CD  00  CD  f&  « 


N N CM  m 

o o or* 

MM  M \n 

mm  m cm 

mm  m cm 

o o o «« 


mo  o o 

n •«  «r 

o rt  n n 

r>  M MM* 

o o o *4 


OHNPi  m m o -c 

O O M M M M O M 

OhND  MfOO-O 

O r*5  M O MOO'*' 

OMMn  Mnoo 

OOOO  O O O -« 


M -r  M M 

o o o o o 

M O ^ M M 
(SI  O CM  M 

n o o *r  o 
o o o o o 


MM-o*r  <o<r  <o<f 

«0  M M O 0-4-4M 

m m m >o  <>o<-o 

OOOO  OOOO 
OOOO  OOOO 


■ON^MN 
Mrinn-r 
o o «o  -o  -c 


ir)'rm-0M«0'O-NMn^m'0M«0'  o m r>  ^/KNooHNn-nn-os 

•*•*•*-*•*•*"*•«•*»*  m mm  m MMMMMMr>mr>rinrjr>n 


173 


36  001646  006302  L0K2I  ASL  R2 

39  001650  016202  001566*  MOV  6NTAB2 ( R2 ) . RZ 

60  001656  000207  RIB  PC 


based  on  the  instruction  execu- 
tion times  of  LSI-11M  (Ref  Appen- 
dix B). 


Figure  E-l . Execution  Time  for  Subroutine  FOURIE 
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APPENDIX  F 


LSI-11  ASSEMBLY  LANGUAGE  CODE  FOR 
POWER  SPECTRUM  SUBROUTINE  (PRSPEC) 
AND  EXECUTION  TIME  COMPUTATION 


appendix  f 


LSI-11  Assembly  Language  Code  for  Power  Spectrum 
Subroutine  (PRSPEC)  ancMlxecution  Time  Computation 
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37  I THE  ARRAY • 
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NOTE:  These  execution  times  are  based 
on  the  instruction  execution  times  of 
LSI-1 1M  (Ref  Appendix  B). 


Total  Execution  Time  = 3*964  msec. 


Figure  F-l . Execution  Time  for  Subroutine  PRSPEC 
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THEREFORE  HE  DOUBLE  XT 


NOTE:  These  execution  times  are 
based  on  the  instruction  execution 
times  of  LSI-1 1 K (Ref  Appendix  B). 


Total  Execution  Time  = 1*3655  msec- 
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APPENDIX  J 


Simulation  Program  for  Testing  Computational 
Accuracy  of  Software  Modules  and  Test  Results 


Source  Code  for  Test  Proaram  (TTFT) 


C * THIS  “R03RAW  IS  USED  to  TEST  THE  COMPUTATI-  • 

C * ONAl  ACCURACY  OE  DIFFFRENT  MOOULES  OF  THE  • 

C » MICROPROCESSOR  BASED  fata  ACQUISITION  ano  * 

C » PROCESSING  SYSTEM  (MI»OAPSt.  THE  PROGRAM  * 

C • ACCENTS  32  DECIMAL  OAT  A VALUES  AN 5 PERFOR-  * 

C » NS  THE  cAME  °PQC  ESSInO  Oc  OATA  AS  NITOAPS.  * 

C » SINNCE  THIS  PROGRAM  WAS  EXECUTED  ON  A CDC-  * 

c • s6ao, which  is  a co-pit  machine  therefore  * 

C • THE  COMPUTATIONAL  RESULTS  WOULD  BE  MORE  * 

C * ACCURATE  AS  COMPARED  WITH  THE  RESULTS  ORT-  • 

C * AINEO  ON  POP-11/03,  WI'HICH  IS  A 15-8IT  MA-  * 

c * chine.* 

C * THE  RESULTS  PRINTEr  OUT  9Y  THIS  PROGPA-  * 

C * H ARE  IN  (TECIMAL  AS  WELL  AS  IN  OCTAL  FORM-  * 

C * AT.  THIS  HAS  9cEN  DONE  FOR  EASY  COMPARISON  • 

C * WITH  THE  RESULT  FROM  ASSEMBLY  L AN 3 J AGE  •- 

C * SOFTWARE  WRITTEN  FOR  MIBDAPS,  * 

c 

c 

c 

PROGRAM  TFFT ( INPUT , OUTPUT , TAPE5=0UT »U T,  T AFE6=INPUTI 
DIMENSION  IPOATA (32) ,PRS(12> 

DIMENSION  F(3?) ,PP<32>,CCF(32> 

OIHENSI ON  C ON F (3?) 

DIMENSION  ERR1(32> ,ERR2<10) ,ERR3<10» 

DIMENSION  LD#TA(3?» ,CF<12> 

DIMENSION  X<32> 

COMPLEX  ROA T A ( 32) 

WRITE  < 5,150) 

150  FORMAT (1H1) 

DO  444  1 = 1,29,4 
PRINT*, “ " 

PRINT*,"  INPUT  FOUR  CATA  VALUES “ 

PRINT*,"  “ 

READ*,  X(I>  , X (!♦ 1)  ,X(1  + 2),X(I*3I 

X(I)=X(I>». 99999  S X(I«-1)=X(I*1>».  993999  S X (1  + 2)  =X  ( I *2)  * .999999 
X(I*3)=X(I»3)* .999999 
444  CONTINUE 
Ms  32 

WRITE( 5,222) 

PRINT*,"  “ 

WRITE (5,232) 

222  FORMATllOX, 30MTA9LE  9ELCW  SHOWS  32  REAL  DATA) 

23?  FORHATdOX,  32HP0INTS  WHICH  ARE  FEO  IV  AS  INPUT X 
PRINT*,"  " 

PRINT*,"  “ 

WRITE( 5 ,121) 

WRITE (5, 1311 

121  FORMAT  OX, 11HCATA  POINTS , 19X , 11H0 AT  A POINTS) 

131  FORMAT  OX, 10HIN  DECIMAL , 20X , 8HIN  OCTAL) 


THIS  PROGRAM  IS  USEO  to  TEST  THE  COMPUTATI- 
ONAL ACCURACY  OF  DIFFFRENT  MOOULES  OF  THE 
MICROPROCESSOR  9ASE0  FATA  ACQUISITION  ANO 
9 ROC E SSI NG  SYSTEM  (MIROAPSt.  THE  PROSRAH 
ACCEPTS  32  DECIMAL  DATA  VALUES  ANO  PERFOR- 
MS THE  CAME  OPOCESSINF  0C  DATA  AS  M 19  OA  PS « 
SINNCE  THIS  PROGRAM  WfS  EXECUTED  ON  A COC- 
5600, which  IS  6 CO-PIT  MACHINE  THEREFORE 
THE  COMPUTATIONAL  RESULTS  WOULD  9E  MORE 
ACCURATE  AS  COMPARED  WITH  THE  RESULTS  OBT- 
AINED ON  POP-1 1/  0 3,  WI'HICH  IS  A 15-811  MA- 
CHINE.* 

THE  RESULTS  PRINTED  OUT  9Y  THIS  PROGPA- 
M ARE  IN  nECIKAL  AC  WELL  AS  IN  OCTAL  FORM- 
AT. THIS  HAS  9 CEII  DOME  FOP  EASY  COMPARISON 
WITH  THE  RESULT  FROM  ASSEMBLY  LANGUAGE 
SOFT WA  R c WRITTEN  FOR  miboaps. 
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v ^ 


WRITE  (>,2421 

242  FORMAT  (5  X,3HOC0,RX, 4HEV EN ,9X , 3HO0 0,5*, 4HEVEN) 

CALL  OORINT (X,H) 

00  10  1=1,32 

10  RDATA ( I) =OMPLX(X(I) ,0,0) 

PRINT*,"  “ 

PRINT*,"  - 
WRITE  1 5,  252) 

WRI TE ( 5, 262 ) * 

252  FORMAT  OX, 3LHTASLE  pfLOV  SHOWS  THE  COMPLEX  DATA) 

262  FORMAT (ex, 36HARRAY  WHICH  IS  FEO  TO  FFT  SUBROUTINE) 

PRINT*,"  “ 

WR:TE(5,121)  1 WRITE(5, 131) 

WRITE(5, 272) 

272  FORMAT (7 X , k UREAL ,9X ,4HIHAG,9X,4HREAL,  5X,  4 HIM AG) 
CALL  OPRINT (RCATA.64) 

N=5  $ XI=1. 

CALL  FFT(RDATA,N,XI) 

00  111  J=l, 32 
REAO* , F( J) 

111  CONTINUE 

00  276  1=1,16 

CON  F (2*1-1) =REAL (ROATA(I) ) 

CONF (2*1) = AIM AStROATA (I) ) 

275  CONTINUE 

no  277  J=l, 32 
ERR1 ( J)  = ( CONF  (J) -F<J))**2 
277  CONTINUE 

PRINT*,"  “ 

PRINT*,"  ~ 

CALL  0PRINT1(R0ATA,32) 

PRINT*,"  “ 

CALL  OPRIHT 1 ( F, 32) 

PRINT*,"  " 

PRINT*,"  “ 

CALL  0PRINT1(ERP1,32) 

CALL  pRSPEC  (RCATA,P(?S,32,ERR2) 

CALL  P0LY(PRS,CF,11,ERR3) 

DO  777  1=1,32 
777  SUM1=SUM1*ERR1(I) 

CONTINUE 
PRINT*,"  “ 

PRINT*,-  “ 

PRINT*,"  RMS  ERROR  FOR  FFT ="  , (S'JMl  732.  ) •*  .5 

00  885  1=1, 10 

SUM2=SU-2»ERR2(I) 

SUM3=SU1?*ERR3(I) 

888  CONTINUE 

PRINT*,"  “ 

PRINT*, “ “ 

PRINT*, "RMS  ERROR  FOR  PWRSPEC=  ",  (S JN2/10 . ) ••. J 
PRINT*,"  " 

PRINT*, “ “ 

PRINT*,"  PM S ERROR  FOR  COF=  ",  ( SUM3/10. ) •• .5 

STOP 

END 
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Source  Code  for  Power  Spoctrun  Subroutine 


c 

SUBROUTINE  PRS»EC (ROATA ,PRS .LEN.ERR2) 

DIMENSION  PRS(ll) 

DIMENSION  PF(  10)  ,ERR2U0> 

COMPLEX  ROATACLEN) 

DO  10  K*l,ll 

PRS<<> =<REAL<ROATACK) >) * * 2. * < AIMAGI RO ATI ( «)  > ) * 
10  CONTINUE 

DO  373  J=l,10 
333  REAO* i “FIJI 
CONTINUE 
DO  4*4  J = l,10. 

441,  ERR2< J> = <PRS< J> -PFI J)»**2 
CONTINUE 
PRINT*,”  " 

PRINT* ,“  " 

CALL  0 PRINT  1 ( FRS, 10) 

PRINT*, “ “ 

PRINT*,-  “ 

CALL  0°RTNT 1 ( PF, 10) 

PRINT*,-  “ 

PRINT*,"  “ 

CALL  Ot>RTNTl(ERR2«10l 
RETURN  t END 
SUBROUTINE  OPRINT1 (X,NP) 

DIMENSION  X CNP) 

PRINT*,-  “ 

PRINT*,"  - 
PRINT*,"  - 
00  10  1=1,  MP 
I1=X(I)»2.**15 
WRITE<5, 100)  X(I),I1 
10  CONTINUE 

100  FORMAT <4X,G13.4,4X,06> 

RETURN 

END 


Source  Code  for  CDF  Subroutine 


SUBROUTINE  POLT (pps,CF,NP,ERR3) 
DIMENSION  PRS  1N°) , CF ( NP) ,CCF(10> 
DIMENSION  E RR3< 10> 

00  10  K= 1 ,NP 
CF«)  = PRStK)**1.6 

10  CONTINUE 

00  555  J=l,10 

555  READ*,  CCF(J) 

CONTINUE 
DO  666  J=l, 10 

665  ERR3  < J > = ( CF (J)-CCF(J))**2 

CONTINUE 
PRINT*,"  " 

PRINT*, " “ 

CALL  0 PRINT  1 (CF, 10) 

PRINT*, “ “ 

PRINT*,"  " 

CALL  OPRINT 1 (CCF,10) 

PPINT*,"  “ 

PRINT*,-  “ 

CALL  OPRINT 1 ( ERRS, 10) 

RETURN  S ENO 


Source  Code  for  FFT  Subroutine 


5 

5 

7 


10 

so 

30 

c 

c 


SUBROUTINE  FFT(X,M,XI>  T COMPLEX  X(1),'J,W,T  1 N=2»*i  t NV2=N/2 
NM1=N-1  t J - 1 t PI  = 3.1M5926535b98  £ 90  7 1=1, NM1 
IF(I.GE.J)  GO  TO  5 I T=X(J)  $ X(J)=X(I)  $ X(I)  »T 
K * NV2 

IF(K.GE.J)  GO  TO  7 1 J=J-K  J K=K/2  £ GO  TO  6 

J=J*<  t 00  20  L=1,M  I LF=2**L  $ LEl=-E/2  * U=(l.iO.) 

M=CEXP(;MPLX  ( 0.  ,-Xl’°I/LEl)  ) $ PO  2C  J* 1 , LEI  ! DO  10  I>J,N,lE 

I P=I  *L  El  J T=X(IP)*U  S X(IP)=X(I)-T 

X(X)  = XU)*  *•  T 

U=U*W  $ IFIXI.GT.O.)  RETURN  $ 00  30  t=t,N 
X(I)=X(I)/N  S RETURN  I ENO 


Source  Code  for  Printing  Out  Test  Results. 


SUBROUTINE  OPRINT (X,NP) 

DIMENSION  X (NP) 

PRINT*  ,“  “ 

PRINT*,-  - 
DO  10  I = 1 , N P , 2 
I1=X(I>*2.**15 

I2=X(I ♦l)*2.**15  . < 

WRITE  (5,100)  X(I)  ,XU»1>  ,11,12 
10  CONTINUE 

10  0 FORMAT  (3X,2G1 1. 1 X , 2 < 3X  , 06)  » 

RETURN  t ENO 


I 

SUBPOUTINF  OPRINT1 (X,NP) 

DIMENSION  X (NP) 

PRINT*,"  “ 

PRINT*,"  " 

PRINT* ,~  - 
DO  10  1*1, NP 
I1  = X(I)*2.*M? 

WRITE (5, 100)  X»I>,I1 
10  CONTINUE 

100  FOPMAT  ('4  X , G 1 3 .4  , MX  , 06) 

I RETURN 

EM) 


210 


o o c>  o n o o 
Or  on  o o o 
i>o  oo  no  p 
o O o o oo  o 
oo  oo  o a o 


OJlir.  r-KSKCOOK  CmcT'K  OK  K tOH  roN  okkkk.sk 

H O o C O O O OO  O O OO  O O O O O O **  o H o oo  oo  ooo 

••*••••••••  •••••••••  at*  iiii  lit 

Uili  U:  U i UJ  U'  UI  U'lii  li  I UI  Uili*  It  u i U i U i It'  Ui  Ui  It*  UJ  UJ  Ul  Li  bJ  U'  U)  UJ  UJ 

KjK-u ) ^ U)  O'  K O UN  -*bC  vPO  iD  rj  n_  cr  (T-  (V  f\i  K- 

O'  o J o'  vp  f'  v.  (Vi  vO  K <0  w J-  r ^ a)  »t  1/v  r »KHNmMi(\KK« 

ciintK^vOO'NK^ortia  o ,£  c K ^ J r\if\jMif'K'C’fyooe 

• o<y  o*k  h ^ hk  (pnK'^(r\Htp^u>,o«o^tn(riM\jKi 


u 

tO  o *(  ||\K  OOOvfiKWJT^-t  ^tOIAiOHON  OJ  K vD  H O in  r(  «c  N 4 2 

■4  c=>  JM  rt  K W iO  (NJ  o ^ HK.  *OHiX)UN^((Cf>.  ^QC  J5  O >P  »-»  vP  *0  O K J '** 

If.  c*  j u\m  y.  k a)  jt  K c »<  x.  if  Mf  in  (*•  o n 4 ins  r->  *-<  ip  »M  .4  ip  r \£  K*  u 

K o ivj  ua.  i»  k o ^ 4 j (/,  j 1-  i .£  « o.  ^ ui  r it  r.  ci  k n j a'  IT  (Shf  O 

(ilonjwK  ^KKuNOMMf.K  -4  (V  (\J  J Ti  O H K 4 Ifi  4 M K M H IIMOK  >- 

oo  o oa  K K KK  OKK  KO  kkokoooookk  k o ook  OK 


III  • • 

UJ  UI  UJ  UJ  ui 

t\j  cd  O'  vO  h in  r)  ro  o r<  (O  -4  <vj  nooo'.tifliAo^fji«ift4(PW*KK 
<\J  o IP  IP  cdWfyff'O'-fK'f'JP'  K tO  K Ul  (»  CH  H |/\  (Ji  a)  K.  Min  vDO>  O 0>(b 

O'  O'  K>K.  HO^iOifiKCTOK  t(J  K-  4on<rinKKlO'cOa04NC'H4 

K-  •«M*-«^rDK»»-*r'K.»p<M  r>  i\j  4>ON.tfy«HMU'^4U'nj4»'JN^N 


»— « ro 

I 4-> 

•-D  =j 
Q. 

UJ  E 
—I  o 

CO  o 

<c 


U.  UJ  UJ  UJ  UI 

CD  (D  ^ *V  oo  O'  ifl  O O 4 |T  W (\JC»  M «c  KJ  o W O 4 nj  lO  J Pi  Pj  O U)  «D  O 

*“•  O Ip  IP  t O *-•  PJ  <T*  <D  Ifi  W KJ  (T<  K lO  K IP  K*  «)  *■«  IP  O'  o N M IP  K-  O'  O O' 

CT*  O'  P''  K-  »-*  CO  rl(OiO  KO*  C K r«4K  4 O -•  « in  K1  M O'  O (O^  M O' vl  ^ 

K.  • Kj  ^ -4  WK  ^ MK  «OW  r (VI  4vP  K ^CVJCO^K'lPfr^lP  PJ  .4  Pj  PJ  C\J 


t0S!?^r^N^2r2'SDHO','H'OHNIfiNOH,,,onH'ONN«« 

UJ  r)NiO>fKOOuJKK^'frt.'KNOKNK'rK^KOrt"f"fn< 

(SJ  NNt>Of*)rir3nN^NKU)N-cy;^^ONK^CN«NNlOnH 
r?  ’"CMiTu')v*Mv1<‘.oocnino-cONuinojiOu')«>OK4ifiK 
ro  -<  -•  -o  -o  k-c>o«-oio-h«<i  ^ (si  >o  -c  in  ^ •«  h k 


3oaoco-.-toooooooooo^ 

jooooooooooaooooooo 

• f + * f ♦ i •♦■♦•■♦■♦♦♦'♦'••■■♦•♦I 
JLULUUJUJUJlJUUJlUUJlUUJUJUJUJUiUiUiUI 
XiG-OCn-ONj',  iTK-iaJM-O-OPrtcOU) 
■cor,  cOM^O'u‘)K>rcC.cC-HCSirt«4>in&*0 
)Oi,i.’,3NOOOK'C-Ov:Cf'4NH^ON 


o-O'T^- 

u*)  C rfl  CSI  f 

5-  3 « f 

r,  zj  « <.  ( 


IOCOOOOOOOOOOOOOOOC 


HrtHHKHNNNONKKK»c-0'OKOK'rHH04NNNOO-rH 
o.’DO'Ok-o^O'J-p,  or^OfKCsiOOP'PjfKK-'OO^r^r^^moooK' 
r)no-CMin4onKKn>iT-cor/n>f4<irKKr,  OKOOrt'O 
KrJMfJ^^iflKO'ON-cOO/i-O'flnr,  •TK'T'OHnuT^NK'O'OO 

rK  |K  fN  rK  fK  PK  PKfK  rK  fK  K K ps. 

M MM  ««  H H «t  H H H H 


i-hhN00--»c0h0-h--0. 
>00000000000000000( 
I I I I ♦ ♦ I I I ♦ I ♦ I I I I ♦ 

1 til  I1J  .11  Ll  1 II?  ^ ’ll  fU  HI  .11  111  .11  .11  1 4 1 II-'  III  , 

fiC'0»‘i»iO»’K(CK>P'C'*KxicC«rO< 
icniN,i1>rOOK«0"(KC'OKN£r: 
«*innlKN03nO'D«lG-'r  -•  fK  — • . 
‘KN(bK*c*MilP»«#Hr3|,/KK^| 


JOOOOOOOOOOOOOOOOOt 


N-fOON4^)ON440N-f'OON4-OON«/.OON4-OON4'00 

KKKOOOG^«««KKNK',(nr,n444-ri1‘lir»k‘-ilT'<!X-04K 

-^-^«rxr>iP^<\t(st<KirvtP4rstP4rs.f>jPtrs,ft(sif't(s.fsiKJtNtC^rjrsjcsirjM(S4{M 

>r<r>r>f<r>r^'4-4>i<r'P^'fv>r-4«r>4<p'/>f>r>r«r>r'r<f>r<r'p'r 


•Scale  factor  CUMSF2  equals  168.89.  The  RKS  error  for  POLYN  equals  .0154 
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APPENDIX  K 


Polynomial  Curve  Fitting  Subroutine 


Polynomial  Curve  Fitting  Subroutine  (PLSCF)(Ref  15) 

Identification: 

PLSCF  Polynomial  Least  Squares  Curve  Fit 
CDC  6600  FORTRAN  Subroutine 

Purpose.  PLSCF  will  compute  the  best  fitting  polynomial  p(x)  to  a 
discrete  set  of  data  points  (x.j,  y^,  i=l,  NP  such  that 

NP  2 

£ [y.r  - p(x.)]  = minimum 

i=l  1 1 

The  coefficients  of  p(x)  will  be  returned  for  polynomials  of  degree  <6. 

For  higher  degree  polynimal  fits  and,  if  requested,  for  polynomials  of 
degree  <6,  the  Chebyshev  polynomial  multipliers  are  returned  (see  remarks). 
Control : 

Dimension  X(NP) , Y(NP) ,W(NP) ,C(NC) ,W0RK(N+1 ) (N+6)/2) ) where 
N=NDEG  , and  NC=N+1  if  NDEG>0,NC=(N+1 ) (N+2)/2  if  NDEG<0 

CALL  PLSCF  ( X , Y , W , NP , NDEG , NMAX , C , I N , XD , X0 , WORK, I ER) 

where 

X(I)  = Input  array  of  distinct  independent  variables. 

Y(I)  = Input  array  of  dependent  variables. 

W(I)  = Input  weights  corresponding  to  each  data  point.  If  W(l)  0, 
all  W(I)>are  taken  as  +1  and  need  not  be  defined. 

NP  = Number  of  given  data  points. 

NDEG  = Degree  of  least  square  polynomial  requested.  If  NDEG  is 

negative,  polynomials  of  all  degrees  0 to  NMAX  will  be  computed. 
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NMAX  = Degree  of  calculated  least  square  polynomial  (0  to  |NDEGj). 
This  is  the  maximum  degree  polynomial  fit  for  which  no  loss 
of  significance  was  indicated. 

C(I)  = Output  coefficient  array  as  follows: 

p(x)  = C(l)  + C(2)  x + C(3)  x2  + ...  + C(N  + 1)  x2 
If  NDEG  is  negative,  the  coefficients  are  returned  as  follows: 
PQ(x)  = C(l) 

P^x)  = C(2)  + C( 3)  x 
P2(s)  = C(4)  + C(5)  + C(6)  x2 

IN  Input  Flag 

= 0 and  NDEG  <6,  the  coefficients  of  p(x)  are  returned  in  C. 
t 0 or  NDEG  6,  the  Chebyshev  polynomial  multipliers  are 
returned  in  C as  follows: 

P(x)  = C(l)  TQ(t)  + C(2)  T] ( t)  + ...  + C(N+1 ) TN(t) 

where 

t = x . XD  + X)  and  T. (t)  - Chebyshev  polynomials 

If  NDEG  is  negative,  the  multipliers  of  each  degree  are 
returned  in  C as  follows: 

P0(x)  = C(l)  T0(t) 

P-,  ( x ) = C( 2)  T0(t)  + C( 3)  T-|  ( t) 

P2(x)  - C(4)  T0(t)  + C(5)  T^t)  + C(6)  T2(t) 

XD,X0  = Output  multiplicative  constant  XD  and  additive  constant  X0 
for  lineary  transformation  of  argument  range. 

WORK(I)  = Working  storage  area  of  ((N+l)  (N+6)/2)) 

IER  Output  Flag 

= 0,  no  errors 
= 1,  errors  in  dimension 
* 2,  coinciding  arguments 
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Methods. 
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- Determining  least  squares  polynomial  fits  to  a set  of  discrete 
data  is  generally  meaningful  only  for  low  degree  polynomials.  Determining 
the  polynomial  fit  for  higher  degrees  in  a straightforward  manner  involves 
solving  a linear  system  of  equations  AC  = B,  where  A is  a positive 
definite  matrix  which  is  frequently  ill-conditioned. 

- This  routine  computes  the  polynomial  p(x)  in  the  form  p(x)  = c-| 

Tq( t)  + C2  T^(t)  + ...  + cm  where  T.  are  the  Chebyshev  polynomials 

and  t is  a linear  transformation  of  the  original  data,  (t  = (2x-xmax 

+ xmin^^xmax"xmin^  = x • XD  + x°)>  This  method  results  in  a remarkable 
improvement  of  the  normal  equations. 

Remarks.  With  this  subroutine,  one  can  obtain  polynomial  fits  for 
degrees  up  to  6.  In  addition,  the  Chebyshev  coefficients  may  be  obtained 
for  any  degree.  If  one  wishes  to  use  a high  degree  fit,  one  can  evaluate 
the  Chebyshev  expansion  for  a specified  argument  x using  subroutine  CHEB 
with  argument  t = x.XD+XO  and  the  calculated  coefficient  vector  of  the 
Chebyshev  expansion. 

The  Chebyshev  expansion  of  the  polynomial  p(x)  gives  a much  better 
indication  of  the  accuracy  of  the  approximation  than  the  coefficient 
vector  of  the  polynomial  itself.  If  the  specified  degree  of  the  poly- 
nomial is  too  high,  the  last  terms  of  the  Chebyshev  expansion  will  be 
uniformly  small  compared  to  the  coefficients  in  front.  The  degree  may 
be  reduced  by  the  number  of  small  trailing  coefficients  without  unduly 
enlarging  the  overall  error.  An  upper  bound  for  the  error  introduced1 
by  neglecting  the  last  terms  of  the  Chebyshev  expansion  is  given  by  the 
sum  of  the  absolute  values  of  these  terms. 
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Storage.  PLSCF  uses  672g  words  of  storage  and  no  COMMON.  In 
addition  to  the  input  and  output  vector  arrays  X,Y,  W and  C,  a working 
storage  area  W0RK((N+l)(N=6)/2)),  N= | NDEG | , must  be  supplied. 


TABLE  K-I 


? 

Constants  for  Polynomial  Ax  + Bx  + C 


Value  of  A 

Value  of  B 

Value  of  C 

Interval  for  x 

Decimal 

(Octal)* 

Decimal 
(Octal )* 

0.00  to  0.04 

3765E-04 

(177776) 

. 1 435E-01 
(000726) 

. 5743E+00 
(044601) 

0.04  to  0.10 

-.3263E-03 

(177765) 

.4974E-01 

(003136) 

.7961E+00 

(062746) 

0.10  to  0.20 

-.2890E-02 

(177641) 

. 1033E+00 
(006470) 

. 5080E+00 
(040405) 

0.20  to  0.30 

-.6618E-02 

(177477) 

.1403E+00 

(010766) 

.41 52E+00 
(032444) 

0.30  to  0.40 

-.11 34E-01 
(177214) 

.1715E+00 

(012763) 

. 3635E+00 
(027205) 

0.40  to  0.60 

- . 1 973E-01 
(176571) 

.2114E+00 

(015416) 

. 31 59E+00 
(024157) 

0.60  to  0.80 

- . 3397E-01 
(175646) 

.2588E+00 

(020441) 

.2762E+00 

(021533) 

0.80  to  1.0 

.5085E-01 

(174575) 

. 3010E+00 
(023205) 

.2499E+00 

(017774) 

♦These  octal  equivalent  values  for  A,  B,  and  C were  used  in  sub- 
routine POLYN.  I i i 
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